Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,269,897 T→G 8.8% noncoding (40/66 nt) rdlB → sRNA antisense regulator affects LdrB translation; proposed addiction module in LDR‑B repeat, with toxic peptide LdrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,269,8970TG 8.8% 308.2 / 15.5 111noncoding (40/66 nt)rdlBsRNA antisense regulator affects LdrB translation; proposed addiction module in LDR‑B repeat, with toxic peptide LdrB
Reads supporting (aligned to +/- strand):  ref base T (51/50);  new base G (6/4);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

GTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGA  >  NC_000913/1269762‑1270036
                                                                                                                                       |                                                                                                                                           
gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAcc                                                                                                                                         >  2:1653742/1‑140 (MQ=37)
     gccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCa                                                                                                                                    <  2:1659786/140‑1 (MQ=34)
          gTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                                                                                               <  2:1254712/140‑1 (MQ=34)
           tCGTGCCAGAAAGTCATGGGAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGc                                                                                                                              >  1:910732/1‑140 (MQ=21)
           tCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGc                                                                                                                              >  1:910772/1‑140 (MQ=34)
            cGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCTTATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTTTTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                             <  1:1561867/140‑1 (MQ=21)
            cGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                             <  1:1065138/140‑1 (MQ=34)
            cGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                             <  1:752120/140‑1 (MQ=34)
            cGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                             <  1:1561907/140‑1 (MQ=34)
            cGTGCCAGAAACTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                             <  1:395535/140‑1 (MQ=17)
              tGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggg                                                                                                                           >  2:386362/1‑140 (MQ=34)
              tGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggg                                                                                                                           >  1:154309/1‑140 (MQ=34)
                      aatCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtct                                                                                                                   <  1:1609547/138‑1 (MQ=11)
                                 aCTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCaa                                                                                                        <  1:134238/140‑1 (MQ=32)
                                  cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                                                                       <  2:311974/140‑1 (MQ=32)
                                     cgcgAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                                                                                                    <  1:45655/140‑1 (MQ=32)
                                      gcgAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                   <  1:338525/140‑1 (MQ=32)
                                      gcgAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                   <  1:1099938/140‑1 (MQ=32)
                                      gcgAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                   <  1:1783687/140‑1 (MQ=32)
                                            gTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacc                                                                                             <  1:906600/140‑1 (MQ=32)
                                               atatAGCCGCCTTGGTGTAATGACAACATTTGGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGGTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCCCCAgg                                                                                          >  2:213132/1‑140 (MQ=255)
                                                    gccgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATaa                                                                                     >  2:436559/1‑140 (MQ=32)
                                                     ccgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATaaa                                                                                    <  1:182712/140‑1 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCACCCAATGCCACCAGGGATAAAGcc                                                                                 >  2:1807362/1‑140 (MQ=17)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:699850/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  2:1807505/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:699911/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:293753/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:1123684/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  2:32037/1‑140 (MQ=32)
                                                        ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:1832894/1‑140 (MQ=32)
                                                         cTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGccc                                                                                <  1:1149573/140‑1 (MQ=32)
                                                            gttgtAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCg                                                                             <  1:783154/140‑1 (MQ=33)
                                                                  aTGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACAt                                                                       >  1:1832067/1‑140 (MQ=37)
                                                                      caacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcg                                                                   <  2:1859637/140‑1 (MQ=37)
                                                                       aacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCAttcccc                                                                  >  1:267945/1‑136 (MQ=255)
                                                                       aacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcgc                                                                  >  1:213982/1‑140 (MQ=37)
                                                                       aacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcgc                                                                  >  2:6358/1‑140 (MQ=37)
                                                                       aacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCCAACATTTCGcc                                                                 >  1:213778/1‑140 (MQ=2)
                                                                         caTTTTGCGGCTATTCTTGAAGTGTATGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCt                                                                <  2:553246/140‑1 (MQ=37)
                                                                              tGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCg                                                           >  1:1151183/1‑140 (MQ=37)
                                                                                cGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGa                                                         <  1:1203250/140‑1 (MQ=37)
                                                                                 ggCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                                                                                                     >  2:674801/1‑79 (MQ=34)
                                                                                 ggCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                                                                                                     <  1:674801/79‑1 (MQ=34)
                                                                                    tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCCCAACATTCCGCCACACCGGATAAt                                                     >  2:896363/1‑140 (MQ=255)
                                                                                     aTTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCTACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAtt                                                    <  1:381764/140‑1 (MQ=255)
                                                                                     aTTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAtt                                                    <  1:535466/140‑1 (MQ=39)
                                                                                            aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                             <  1:268240/140‑1 (MQ=255)
                                                                                            aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                             <  1:1706743/140‑1 (MQ=255)
                                                                                            aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATACCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                             <  2:968756/140‑1 (MQ=255)
                                                                                              gtgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTgg                                           <  1:1832672/140‑1 (MQ=255)
                                                                                               tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:1584122/140‑1 (MQ=255)
                                                                                               tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:257049/140‑1 (MQ=255)
                                                                                               tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:257039/140‑1 (MQ=255)
                                                                                               tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:1584134/140‑1 (MQ=255)
                                                                                                         cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGccccc                                                                              >  2:1514799/1‑94 (MQ=33)
                                                                                                           aGAATAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGAGCGGGGGTTTTCTCTTTCCAGCAACCAATGACACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATACCGGCTTGGTGAGGATAATACg                              >  2:1674507/1‑140 (MQ=11)
                                                                                                                aGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCa                         >  1:1699741/1‑140 (MQ=40)
                                                                                                                  cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAt                       <  2:1123684/140‑1 (MQ=40)
                                                                                                                   ccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAtt                      <  2:172096/140‑1 (MQ=40)
                                                                                                                   ccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAtt                      <  1:439513/140‑1 (MQ=40)
                                                                                                                   ccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAtt                      <  1:438862/140‑1 (MQ=40)
                                                                                                                      cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                               >  2:964664/1‑112 (MQ=39)
                                                                                                                      cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                               <  1:964664/112‑1 (MQ=39)
                                                                                                                      cGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAGTCCGGCTTGGTGTGGATACTACGTCTCAATTCAt                   >  2:1843373/1‑140 (MQ=17)
                                                                                                                       gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAttcttc                  >  1:1346260/1‑140 (MQ=255)
                                                                                                                       gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATc                  >  1:845915/1‑140 (MQ=40)
                                                                                                                       gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATc                  >  1:403705/1‑140 (MQ=40)
                                                                                                                         tCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATctt                <  1:1081172/140‑1 (MQ=40)
                                                                                                                         tCTGTTGTCAGGATTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATctt                <  1:1078231/140‑1 (MQ=255)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAGAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:192555/1‑140 (MQ=255)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGa                                     >  2:871086/1‑116 (MQ=39)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGa                                     <  1:871086/116‑1 (MQ=39)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCATTTCATcttcac             >  2:922057/1‑140 (MQ=255)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:795428/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  1:1695957/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  1:1015144/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  1:419528/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  1:1796941/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:194500/1‑140 (MQ=40)
                                                                                                                            gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTACGGGTTGGGGTGGGTACTACGTCTCAATTCTTcttcac             >  1:1206906/1‑140 (MQ=14)
                                                                                                                            gttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCaa                        >  1:1678577/1‑129 (MQ=25)
                                                                                                                            gttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCaa                        <  2:1678577/129‑1 (MQ=25)
                                                                                                                             ttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcact            >  2:1497870/1‑140 (MQ=40)
                                                                                                                             ttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcact            >  2:1496264/1‑140 (MQ=40)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:1662598/110‑1 (MQ=17)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          >  2:971785/1‑110 (MQ=17)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:1318421/110‑1 (MQ=17)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          >  2:1662598/1‑110 (MQ=17)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          <  1:971785/110‑1 (MQ=17)
                                                                                                                             ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                          >  2:1318421/1‑110 (MQ=17)
                                                                                                                              tgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcactt           <  1:1640826/140‑1 (MQ=40)
                                                                                                                                  aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                                                              <  2:1651102/85‑1 (MQ=33)
                                                                                                                                  aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                                                              >  1:1651102/1‑85 (MQ=33)
                                                                                                                                  aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATc       >  1:1860295/1‑140 (MQ=40)
                                                                                                                                  aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATc       >  1:1860565/1‑140 (MQ=40)
                                                                                                                                  aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATc       >  1:1860228/1‑140 (MQ=40)
                                                                                                                                       ttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa  >  1:526497/1‑140 (MQ=38)
                                                                                                                                       ttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa  >  1:1468003/1‑140 (MQ=38)
                                                                                                                                       ttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa  <  2:112323/140‑1 (MQ=38)
                                                                                                                                       |                                                                                                                                           
GTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGA  >  NC_000913/1269762‑1270036

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: