Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,269,897 T→G 8.9% 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.9% 203.1 / 6.1 74noncoding (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 (47/20);  new base G (4/3);  total (51/23)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGA  >  NC_000913/1269759‑1270036
                                                                                                                                          |                                                                                                                                           
tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                                                                                            <  1:1425957/140‑1 (MQ=37)
tCGGTGCAGCCAGGTCGTGCCAGAAACTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTTGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                                                                                            <  2:936793/140‑1 (MQ=12)
   gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTGTCAAGATTAGCCCCCGTTCTGTTGTCAGGTGTTAcc                                                                                                                                         >  2:1309850/1‑140 (MQ=255)
         ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                                                                                                                                   >  1:701242/1‑140 (MQ=34)
         ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                                                                                                                                   >  1:701428/1‑140 (MQ=34)
         ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                                                                                                                                   >  1:701434/1‑140 (MQ=34)
         ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                                                                                                                                   >  1:185016/1‑140 (MQ=34)
              tCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGc                                                                                                                              >  2:1294013/1‑140 (MQ=34)
                 tGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggg                                                                                                                           >  1:1320446/1‑140 (MQ=34)
                                  aaccTGCGCGAGCGACATATAGCCGCCTTGTTGTAACGACAACAATTTGCGGCTATTCTTGAAGACTCCGGATTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                                                                                          <  2:298364/137‑1 (MQ=11)
                                  aaaCTGCGCGAGCGTCATATAGCCGCCGTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGGGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                                                                                          <  2:898235/140‑1 (MQ=17)
                                     cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                                                                       >  2:699151/1‑140 (MQ=32)
                                     cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                                                                       >  2:483627/1‑140 (MQ=32)
                                     cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                                                                       >  2:699471/1‑140 (MQ=32)
                                     cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                                                                       >  2:699420/1‑140 (MQ=32)
                                      tGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTCTTGTTACAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                                                                      <  1:1230118/140‑1 (MQ=12)
                                        cgcgAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                                                                                                    >  1:834385/1‑140 (MQ=32)
                                              cGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccac                                                                                              <  1:1350202/140‑1 (MQ=32)
                                               gTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacc                                                                                             <  1:1627453/140‑1 (MQ=32)
                                               gTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacc                                                                                             <  1:905092/140‑1 (MQ=32)
                                                tCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacca                                                                                            <  1:1291204/140‑1 (MQ=32)
                                                          gccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGc                                                                                  <  2:185016/140‑1 (MQ=32)
                                                           ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  2:1299/1‑140 (MQ=32)
                                                           ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcc                                                                                 >  1:445818/1‑140 (MQ=32)
                                                             ttgttgTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcccc                                                                               >  2:1374946/1‑140 (MQ=32)
                                                             ttgttgTAATGACAACATTATGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGcccc                                                                               >  2:853208/1‑140 (MQ=17)
                                                                ttgaAATGACAACATTTTGCGGCTATTCTTGAAGTTTCTGGTTTCAATATAAGCCCACGTTCAGTAGACAGGTAGTACCTCTCAACGAGCGGGGGTTTTCTCTATCCAGCAACCAAAGCCACCAGGGATAAAGCCCCCGc                                                                            >  1:1674038/1‑140 (MQ=11)
                                                                          aacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCCCACCATTccgc                                                                  >  2:1021462/1‑140 (MQ=255)
                                                                                      cTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATaa                                                      >  1:237545/1‑140 (MQ=37)
                                                                                       tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAt                                                     >  1:524233/1‑140 (MQ=38)
                                                                                       tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAt                                                     >  2:1174256/1‑140 (MQ=38)
                                                                                       tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAt                                                     >  1:1052897/1‑140 (MQ=38)
                                                                                       tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAt                                                     >  1:1052905/1‑140 (MQ=38)
                                                                                       tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCAACACCGGATact                                                     >  1:1422767/1‑138 (MQ=255)
                                                                                        aTTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCATCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAtt                                                    <  1:1595563/140‑1 (MQ=255)
                                                                                             tGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                               <  1:705809/140‑1 (MQ=255)
                                                                                             tGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                               <  2:321079/140‑1 (MQ=255)
                                                                                               aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                             <  1:1299/140‑1 (MQ=255)
                                                                                                aGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGGTTTACCTCTCAACGTGCGGGGGGTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTg                                            >  2:477238/1‑140 (MQ=255)
                                                                                                   gtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAgg                                                                                          >  2:844370/1‑91 (MQ=33)
                                                                                                   gtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAgg                                                                                          <  1:844370/91‑1 (MQ=33)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCATTGCGCCTCACCGGATAATTCCGGCGTGgtgt                                        >  1:207415/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:1248684/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:1248701/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:1248186/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:498041/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  2:433963/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:834947/1‑140 (MQ=255)
                                                                                                    tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:987560/1‑140 (MQ=255)
                                                                                                    tCTGGTATCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCATTGCGCCTCACCGGATAATTCCGGCTTGgtgt                                        >  1:836873/1‑140 (MQ=255)
                                                                                                     cTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGtgtg                                       <  1:232047/140‑1 (MQ=255)
                                                                                                      tGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCGTGGTGTgg                                      >  2:205512/1‑140 (MQ=255)
                                                                                                      tGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGAAACATTGCGCCAAACCGGAAAATTCCGGCTTGGTGTgg                                      >  2:205473/1‑140 (MQ=255)
                                                                                                       tggttCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGa                                     <  2:204681/137‑1 (MQ=1)
                                                                                                            cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCACCGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGTTAATTCCGGCTTGGTGTGGAtacta                                <  2:588563/140‑1 (MQ=255)
                                                                                                            cAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGAtacta                                <  2:1374818/140‑1 (MQ=17)
                                                                                                               gATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGt                             <  1:338483/140‑1 (MQ=38)
                                                                                                                  tAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGTCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGtctc                          >  2:123844/1‑140 (MQ=255)
                                                                                                                  tAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGtctc                          >  2:488949/1‑140 (MQ=40)
                                                                                                                     cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAAt                       <  2:621542/140‑1 (MQ=40)
                                                                                                                         cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCAt                   >  2:1246833/1‑140 (MQ=40)
                                                                                                                          gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACTTTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATc                  <  1:1447996/140‑1 (MQ=255)
                                                                                                                          gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATc                  <  1:1502332/140‑1 (MQ=40)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCg                                                                             <  2:451012/76‑1 (MQ=32)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACTTTGCGACTAACCGGAAAATTCCGGGTTGGTGTGGATAATACGACTAAATTCATcttcag             >  2:855047/1‑139 (MQ=11)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:640874/1‑140 (MQ=40)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:402338/1‑140 (MQ=40)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:738199/1‑140 (MQ=40)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:908681/1‑140 (MQ=40)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACAACGTCTCAATACATctacac             >  2:308545/1‑140 (MQ=255)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGTGCGGGGATTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  1:387422/1‑140 (MQ=255)
                                                                                                                               gttgtCAGGTTTTACCTCTCAACGAGCGGGGGTTTTCTCAATCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCTGATAATACAGGCTTGTTGTGGATACAACGTCTCAATTCATcttcaa             >  2:1572015/1‑139 (MQ=11)
                                                                                                                               gttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:791584/1‑140 (MQ=21)
                                                                                                                               gttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAAATCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac             >  2:319632/1‑140 (MQ=17)
                                                                                                                                ttgtCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCATTGCGCCTCACCGGATAAATCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcact            >  2:1648228/1‑140 (MQ=12)
                                                                                                                                   tCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcacttca         <  1:319632/140‑1 (MQ=21)
                                                                                                                                        ttttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCATTGCGCCTCACCGGATAATACCGGCTTGGTGTGGATACTACGACTCAATTCAACTACACTACAACCCt    >  2:1639251/1‑140 (MQ=11)
                                                                                                                                          ttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa  >  1:1615409/1‑140 (MQ=38)
                                                                                                                                          ttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGAAAAATGGCGCAACACCGCATAATTCCGGGGTGTGGTGTATAATACTCCACAAAACATCTTCCCTAAAACCCTGa  >  1:1615479/1‑140 (MQ=2)
                                                                                                                                          |                                                                                                                                           
TCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGA  >  NC_000913/1269759‑1270036

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