Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,450,008 T→C 7.7% I442V (ATC→GTC)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,450,0080T→C 7.7% 52.3 ‑0.7 ‑0.2274I442V (ATC→GTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (118/135):  total (126/148)
Fisher's exact test strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.83e-01

Only 100 of 287 total aligned reads displayed.

TCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAA  >  NC_007779/3449908‑3450108
                                                                                                    |                                                                                                    
tCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGAt                                                                                                      <  1:4161643/101‑1 (MQ=255)
tCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCACCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCACCTGGATAGCTTTACCTTCGAt                                                                                                      >  1:8954488/1‑101 (MQ=255)
   aGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAg                                                                                                   <  1:8816440/101‑1 (MQ=255)
     cGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTa                                                                                                 <  1:6930281/101‑1 (MQ=255)
     cGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTa                                                                                                 <  1:11671716/101‑1 (MQ=255)
          gCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGt                                                                                            >  1:11889890/1‑101 (MQ=255)
           ccgTAAGTGAACAGCCATCTGGTCCCCATCGAAGTCGGCGTTATATGCCGCAGAAACCAGCGGGTGCAGCTGGATAGCTTTCCCTTCGATCAGTCCCGGtt                                                                                           <  1:42144/99‑1 (MQ=255)
           cGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGtt                                                                                           >  1:6256353/1‑101 (MQ=255)
           cGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGtt                                                                                           <  1:6971223/101‑1 (MQ=255)
                 gTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATg                                                                                     >  1:408055/1‑101 (MQ=255)
                   gAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGcc                                                                                   >  1:7639572/1‑101 (MQ=255)
                    aaCAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCt                                                                                  >  1:12810551/1‑101 (MQ=255)
                    aaCAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCt                                                                                  <  1:5395223/101‑1 (MQ=255)
                    aaCAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCt                                                                                  <  1:3812329/101‑1 (MQ=255)
                          cATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATAc                                                                            <  1:1107121/101‑1 (MQ=255)
                          cATCTCGTCACCATCGAAGTCGGCGTTATATGCCGCCCAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGCTCAACTGCCTGGATAc                                                                            >  1:6526797/1‑101 (MQ=255)
                           aTCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATAcc                                                                           >  1:7143513/1‑101 (MQ=255)
                             cTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTTGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCa                                                                         >  1:13935089/1‑101 (MQ=255)
                             cTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCa                                                                         <  1:10391151/101‑1 (MQ=255)
                              tGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAg                                                                        <  1:6457003/101‑1 (MQ=255)
                              tGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAg                                                                        >  1:11264019/1‑101 (MQ=255)
                                 tCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACg                                                                     <  1:1509505/101‑1 (MQ=255)
                                 tCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACg                                                                     <  1:13543349/101‑1 (MQ=255)
                                 tCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCCGGATACCCAGACg                                                                     <  1:4824469/101‑1 (MQ=255)
                                   aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGt                                                                   <  1:8086315/101‑1 (MQ=255)
                                   aCCATCGAAGTCGGCGTTATATGCCGCACAAACCAACGGGTGCCGCTGGATAGCATTACCTTCGATCATTACCGGGTTAAATGCCTGGATACCCAGACGGt                                                                   >  1:12206913/1‑101 (MQ=255)
                                    ccATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTg                                                                  <  1:7564272/101‑1 (MQ=255)
                                     cATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGc                                                                 <  1:9830563/101‑1 (MQ=255)
                                      aTCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCa                                                                >  1:4759431/1‑101 (MQ=255)
                                      aTCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCa                                                                >  1:3088556/1‑101 (MQ=255)
                                        cGAAGTCGGCGTTATATGTCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCaga                                                              <  1:332844/101‑1 (MQ=255)
                                         gAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCagag                                                             >  1:13318953/1‑101 (MQ=255)
                                           aGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTc                                                           >  1:9328535/1‑101 (MQ=255)
                                           aGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTc                                                           <  1:4281865/101‑1 (MQ=255)
                                           aGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTc                                                           <  1:2018154/101‑1 (MQ=255)
                                           aGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTc                                                           <  1:1833564/101‑1 (MQ=255)
                                             tCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCgg                                                         <  1:271124/101‑1 (MQ=255)
                                              cGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGt                                                        <  1:5025907/101‑1 (MQ=255)
                                                gCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCACAGACGGTGc                                                      >  1:6030743/1‑101 (MQ=255)
                                                 cGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCa                                                     <  1:38694/101‑1 (MQ=255)
                                                  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCAc                                                    >  1:14123212/1‑101 (MQ=255)
                                                  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCAc                                                    <  1:12594675/101‑1 (MQ=255)
                                                  gTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCAc                                                    <  1:2790058/101‑1 (MQ=255)
                                                    tataTGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACgg                                                  >  1:1037534/1‑101 (MQ=255)
                                                     atatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGt                                                 >  1:7338217/1‑101 (MQ=255)
                                                     atatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGt                                                 <  1:3642637/101‑1 (MQ=255)
                                                      tatGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGtt                                                <  1:223605/101‑1 (MQ=255)
                                                       atGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTc                                               <  1:5457237/101‑1 (MQ=255)
                                                       atGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTc                                               <  1:14328357/101‑1 (MQ=255)
                                                       atGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTc                                               >  1:13606056/1‑101 (MQ=255)
                                                       atGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTc                                               <  1:8622934/101‑1 (MQ=255)
                                                        tGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTca                                              <  1:12003423/101‑1 (MQ=255)
                                                          cctcacatgccAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTcagc                                            <  1:8145173/92‑1 (MQ=255)
                                                          ccGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTcagc                                            <  1:11405257/101‑1 (MQ=255)
                                                             cacaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGt                                         <  1:1975341/101‑1 (MQ=255)
                                                              gctatccAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                        >  1:989531/6‑101 (MQ=255)
                                                              acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                        >  1:1650133/1‑101 (MQ=255)
                                                              acaAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTa                                        <  1:5882029/101‑1 (MQ=255)
                                                               caAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAc                                       <  1:8461577/101‑1 (MQ=255)
                                                                aaaCCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGCTGCAGAGTCGGAGCACGGTTCAGCAGTAcc                                      >  1:2900058/1‑101 (MQ=255)
                                                                aaaCAAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTAcc                                      >  1:3907486/1‑101 (MQ=255)
                                                                 tgccAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCCAACGCCTGAATACCCAGGCGGGGCAGCGTAGGTGCACGGTTCAGCAGGACCg                                     >  1:10777983/3‑101 (MQ=255)
                                                                 aaCCAGCGGGCGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGCACCg                                     >  1:6895639/1‑101 (MQ=255)
                                                                    cAGTGGGTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGgt                                  <  1:10601924/101‑1 (MQ=255)
                                                                    cAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgt                                  >  1:6143545/1‑101 (MQ=255)
                                                                     aGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtg                                 >  1:2988051/1‑101 (MQ=255)
                                                                      gtgggTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGgtgt                                <  1:7796063/99‑1 (MQ=255)
                                                                      gtgggTGCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGGGCAGCGTAGGTGCACGGTTCAGCAGGACCGGgtgt                                >  1:2104192/3‑101 (MQ=255)
                                                                      gCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGgtgt                                <  1:1907612/101‑1 (MQ=255)
                                                                       cGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGtt                               <  1:14056222/101‑1 (MQ=255)
                                                                         ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg                             >  1:10254218/1‑101 (MQ=255)
                                                                         ggTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTcg                             <  1:4573416/101‑1 (MQ=255)
                                                                            gCAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcggta                          <  1:13054887/101‑4 (MQ=255)
                                                                             cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa                         <  1:9633913/101‑1 (MQ=255)
                                                                             cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa                         >  1:13178364/1‑101 (MQ=255)
                                                                             cAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGa                         <  1:7509390/101‑1 (MQ=255)
                                                                             cAGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcggtaa                         <  1:11752288/101‑5 (MQ=255)
                                                                              aGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGAt                        >  1:11490872/1‑101 (MQ=255)
                                                                              aGCTGAATTGCCTTACCCTCAACCAGGACCGGCTCGAACGCCTGAATACCCAGGCGGTGCAGCGTAGGTGCACGGTTCAGCAGGACCGGGTGTCcggtaat                        <  1:3739033/101‑6 (MQ=255)
                                                                                cTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATaa                      >  1:1234390/1‑101 (MQ=255)
                                                                                cTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATaa                      >  1:12517571/1‑101 (MQ=255)
                                                                                  ggATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACt                    <  1:10983965/101‑1 (MQ=255)
                                                                                  ggATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACt                    >  1:12700765/1‑101 (MQ=255)
                                                                                   gATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACtt                   <  1:261335/101‑1 (MQ=255)
                                                                                    aTAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCGGCAGTACCGGGTGTTCGCGGATAACTTc                  <  1:3425725/101‑1 (MQ=255)
                                                                                     tAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCg                 <  1:11803030/101‑1 (MQ=255)
                                                                                     tAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCg                 <  1:2160146/101‑1 (MQ=255)
                                                                                        cTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTcc              <  1:8298210/101‑1 (MQ=255)
                                                                                        cTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTcc              <  1:13740292/101‑1 (MQ=255)
                                                                                        cTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGCGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTcc              >  1:470065/1‑101 (MQ=255)
                                                                                         tttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCa             <  1:8008047/101‑1 (MQ=255)
                                                                                          ttACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAg            <  1:6390128/101‑1 (MQ=255)
                                                                                             ccTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGat         <  1:552567/101‑1 (MQ=255)
                                                                                             ccTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGat         <  1:10687138/101‑1 (MQ=255)
                                                                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat       <  1:1730591/101‑1 (MQ=255)
                                                                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat       <  1:9149309/101‑1 (MQ=255)
                                                                                               ttCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGatat       >  1:13465664/1‑101 (MQ=255)
                                                                                                tCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATc      <  1:5088850/101‑1 (MQ=255)
                                                                                                   aTCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCa   >  1:12044490/1‑101 (MQ=255)
                                                                                                    tCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCaa  >  1:4565050/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
TCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGAAGTCGGCGTTATATGCCGCACAAACCAGCGGGTGCAGCTGGATAGCTTTACCTTCGATCAGTACCGGTTCAAATGCCTGGATACCCAGACGGTGCAGAGTCGGTGCACGGTTCAGCAGTACCGGGTGTTCGCGGATAACTTCGTCCAGGATATCCCAA  >  NC_007779/3449908‑3450108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.