Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,579,657 A→C 23.7% V62G (GTG→GGG)  symE ← toxic peptide regulated by antisense sRNA symR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,6570AC23.7% 90.5 / 3.7 54V62G (GTG→GGG) symEtoxic peptide regulated by antisense sRNA symR
Reads supporting (aligned to +/- strand):  ref base A (36/3);  new base C (5/9);  total (41/13)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  NC_000913/4579424‑4579889
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
gtaattttttgaggggggggAAATTTTTTTGTGAGTCGTGGGAAGTCCCCCGGGGGTTAACCGGGAGTAACCCGTTTTACGCTACTTTCTGTTTACCGGCAATCCCCCCAAAAACCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCGGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTg                                                                                                                                                                                                                                        <  1:358765/224‑1 (MQ=255)
gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGCCAATACAGCCTTCCATCACTTTg                                                                                                                                                                                                                                        >  2:610974/1‑239 (MQ=255)
gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAAGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGGAGGGACAAAAACGCCTTCCACAACtttt                                                                                                                                                                                                                                        >  1:730179/1‑238 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCCGTTCGCTCTCCTCGGCGGAGGGTGGGTGGGCGGGGGGGCCAATACAGCCCTCCATCACTTTGAAATCTAccc                                                                                                                                                                                                                             >  1:618761/1‑238 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCATTTGACAttaacc                                                                                                                                                                                                                              >  2:180033/1‑233 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGGTGGGCGGGTAGGACAATAAAGCCTTCCATCACCTTGGCATCTAcc                                                                                                                                                                                                                              >  1:69914/1‑238 (MQ=255)
             gAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCCTTTGACATCTCccgccg                                                                                                                                                                                                                          >  2:374625/1‑239 (MQ=255)
               gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACAACAACCGCCGt                                                                                                                                                                                                                         >  2:874211/1‑238 (MQ=255)
               gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACACTTTGCACACCTGGCGTAGCGACTGCCTCAGTTCGCTCCCCTCGGGGGCGGGGGGTTGGGCGGTAAGGCCAACACAGCCTTCCACCCACTTGAAATCTAcccccgg                                                                                                                                                                                                                         >  1:526232/1‑237 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGGGGGGGGGTGGGGGGGAGGGCCAATAAAGCCTTCCCCCCCTTTGACATCACCCGCCGTg                                                                                                                                                                                                                        >  1:467750/1‑238 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGGCAATACAGCCTTCCATCACTTTGCCATCTACCGCCGTg                                                                                                                                                                                                                        >  2:212738/1‑238 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTg                                                                                                                                                                                                                        >  2:710451/1‑238 (MQ=255)
                ggggAAAGTTTCTACTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTg                                                                                                                                                                                                                        >  2:750574/1‑238 (MQ=255)
                          tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGCCAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCaaa                                                                                                                                                                                                             >  2:570379/1‑239 (MQ=255)
                                  gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGGGGCAAAACcggcg                                                                                                                                                                                                      >  2:279026/1‑238 (MQ=255)
                                  gTCAGTGGAAGGCTACCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCATTAAACGCCTGCACCTGCTTTTGTTTACGCGCCGCCAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGCGGTTGGGCGGCGCGTCCAATACAGCCTTCCTTCACTTTGACATCTACCGCCGTGCCAGTAGCAAAACcggcg                                                                                                                                                                                                      >  2:125741/1‑238 (MQ=255)
                                                cccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCCCCGTGCCAGTGGCAAAACCGGCGGCTTCCAACCCCTg                                                                                                                                                                                        >  2:460628/1‑238 (MQ=255)
                                                 ccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCCTCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCATCCCCTGa                                                                                                                                                                                       >  2:731156/1‑238 (MQ=255)
                                                            ccGGGAGTAACTCCGTTTACGCGTCTTTTGTTTTACCGCAAATCACCCAAAAAAACCCCTGCCCCTGTTTTTGTTTCCGCCCCGGCAGTTGGCACACTGGGGGGAGCGACTGAACCAGTCCGCTCTCCCCGGCGGCGGGTGGTGGGGCGGTGAGGAAAATCCAGCCTCCCATCACTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTg                                                                                                                                                                           <  1:523816/239‑1 (MQ=255)
                                                                 aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCAAGCCACTGACCTTTCAGGGTGATGGc                                                                                                                                                                      >  2:233168/1‑239 (MQ=255)
                                                                 aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGCGTGGGCGGGACGGACAATCCAGCCTTCCATCACTTTGTCATCTACCGCCGTGCCCGGGGCAAAACCGACGGCTTCCACCCAATGACCCTTCAGGGTAATgac                                                                                                                                                                      >  1:176042/1‑237 (MQ=255)
                                                                        cctttttccctttttttttttctccGTCAATCCTCCAAAAAAACGCCCGCCCCTTTTTTTGGTTCCGCGCCGACTGTTGGCCCCCCGGGCGTAGCGTCAGAATCAGTTCGCTCCCCTCGGCGGCGGGGGGTGGGGGGGGGAGGAAAATACAGCCTTCCATCCCTTTGCCATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATa                                                                                                                                                               <  1:460628/216‑1 (MQ=255)
                                                                          tGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCACCCACTTACCTTTCCGGGGGATgcgggggaaaag                                                                                                                                                             >  1:802403/1‑228 (MQ=255)
                                                                                      ttCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCCCTGACCTTTCAGGGGGATGGCGGGAATACGGCTGTAATCCgg                                                                                                                                                >  2:445083/1‑239 (MQ=255)
                                                                                              aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGGATACGGCTGGAATCCGGGTGGCgg                                                                                                                                         >  2:452424/1‑237 (MQ=255)
                                                                                              aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGAGGGCGGGAAAACGCGCTGAAAACCGGGGAgga                                                                                                                                          >  1:422850/1‑236 (MQ=255)
                                                                                              aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCACCCACTGCCCTTTCCGGGGGCGGGCGGGGAAACGGCCGTAATCCGGGGGGCGa                                                                                                                                         >  1:880184/1‑238 (MQ=255)
                                                                                              aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGATTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACAGCATCAGTTCGCCCTCCTCGGCGTCGGGGGGTTGGGCGGTGAGGGCAATGCAGCATTCCATCCCTTTGACATCTACCGCCGTGCCAGTGGCGAAACAGGAGGCTTCCAGCGCATGACCTTTCCGGGGGAGGGCGGGATAACGGCTGGAATCCgggggggca                                                                                                                                         >  1:662247/1‑232 (MQ=255)
                                                                                                   cAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGCGGGGAATACGGCTGTAAACCCGGGAGCGAccccca                                                                                                                                   >  1:332963/1‑234 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGCCCTTTCAGGGTGATGGGGGGAATACGGCCGTAATCCGGGTAGCGACTCGCaaaaaccgggggggcaag                                                                                                                  >  1:747405/1‑225 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGGTGGCGGGGATAAGGCTGTAAACCGGGCAGCCACCCGCCAAAACGACGGGGaaaaa                                                                                                                 >  1:102816/1‑235 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGGGaaccccaaaaccaccgggggcctt                                                                                                                  >  2:715111/1‑216 (MQ=255)
                                                                                                                         gCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGGAATCCGGGTAGCGACTCGCATAACCAACGGTGACATGACgg                                                                                                             >  2:186802/1‑239 (MQ=255)
                                                                                                                         gCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGAGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGGATACGGCTGGAATCCCGGGAGGGACCCGCCTaaaccgcgggggaatcaggg                                                                                                              >  1:602589/1‑221 (MQ=255)
                                                                                                                          cACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCAAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGAACCCCTTAACCGACGGTGACATGACGGt                                                                                                            >  2:125141/1‑239 (MQ=255)
                                                                                                                           cccggCTTTTGTTTACGGCCCGAAAGTTTGCCACCCGGGCGAGGGCAAGGCATCAGTTCGCTCTCCGCGGGGGGGGGTGGTTGGGGGGTGAGGAAAATACAGCTTTCCATCACTTTGACATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGtt                                                                                                           <  2:571765/235‑1 (MQ=255)
                                                                                                                              tGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGCCGGTTAtt                                                                                                        >  2:416819/1‑239 (MQ=255)
                                                                                                                                    tGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCgg                                                                                                  >  2:130179/1‑239 (MQ=255)
                                                                                                                                                       tGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATGCAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGACTACGGCTGTAATCCGGGTAACGACTCGCATAAACGAAGGGGGACGAGAGGTTTTTTGCCGGGGGAACTTCCGGTTCGAAc                                                                               >  1:336741/1‑239 (MQ=255)
                                                                                                                                                               tGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACTTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGGGACATGACGGTTATTTGCCGGGGGGACTTCTGCTTCGAACGGtttggc                                                                       >  2:509635/1‑235 (MQ=255)
                                                                                                                                                                       gagcggCATCAGTTCCGTCCCCCGGGCGGGGGGTGGTTGGGCGGTGAGGAAAATCCAGCCTCCCATCCCTTGGACATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAAt                                                               <  2:102816/234‑1 (MQ=255)
                                                                                                                                                                            gcagtaGTTCGCTCGCGTCGGCGGGGGGTGGTGGGGCGGTGATGACAATCCAGCCTTCCATCCTGTTGACATCTACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCgtg                                                          <  2:440662/234‑1 (MQ=255)
                                                                                                                                                                                  gTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCACGTGGCAAAACCGGCGGCTTCCAGCCCCTGCCCTTTCAGGGTGATGGCGGGAATACGGCCGTAATCCGGGTCGCGACTCCCCTAACCCCCCGTGCCATACCAGTTATTTTCCGGGGGGACGTCTGCGTCGCACGGGTGTGCAACAGAAAGCgtgtgcggg                                                    >  1:430592/1‑234 (MQ=255)
                                                                                                                                                                                            ctggggggggggTTGTTGGGGGGTAGGGCAAATCCACCCTTCCACCCCTTGGACTGCTCCCGCCGGGGCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTa                                          <  2:430592/230‑1 (MQ=255)
                                                                                                                                                                                                 gcgggggTTGTGGGGGGGTGGGACAAATCCCGCCTTCCACCCTTTGCACTTCACCCGCCGGGCAAGGGGCAAAACCGGCGGCTGCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATtctc                                     <  1:710451/235‑1 (MQ=255)
                                                                                                                                                                                                  ggcggGTGGTGGGGGGGTGAGGACATTACAGCCTTCCACCCTTTGGACATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTGCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATtctc                                     <  2:764659/238‑1 (MQ=255)
                                                                                                                                                                                                      ggtgttTGGGCGGTGAGGACAATCAACCCTCCCATCCCTTTGCCATCTACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAgg                                 <  2:400553/238‑1 (MQ=255)
                                                                                                                                                                                                      ggtggtTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTTTGCAATAGCATGCGTGTCAGTCATAACTGCTATTTTCCAGGa                                >  2:464561/1‑239 (MQ=255)
                                                                                                                                                                                                                ggggaaGGAAAATCCAGCCTCCCATCCTTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgatt                       <  2:513305/233‑1 (MQ=255)
                                                                                                                                                                                                                 ggggtaggaaaatcacgCCTTCCATCTCTTTGCCATGTGCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgatt                       <  2:618761/222‑1 (MQ=255)
                                                                                                                                                                                                                     aGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCCTAACCGACGGGGACATGCCGGTTATTTGCCGGGGGGACTTCTGCTTTGAAAGGTTGGGCAATAGAATGCGTGTCAGTCATAACAGCCATTCCCCAGGAAAAgggattgggatt                 >  1:695099/1‑239 (MQ=255)
                                                                                                                                                                                                                         aaaaatcAGCCTTCCACCCATTGGAACTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAg               <  2:832193/232‑1 (MQ=255)
                                                                                                                                                                                                                                  ccTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGgt     >  2:487995/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGGGACATGACGGTTTTTTGCCGGGGAGACTTTTGCTTCGAACGGTTGTGCAAAAGAATACGGGTCCGGCATAAATTCCATTATCCAGGAAAAGTGATTGAGAATAACGAAGCgggggga  >  1:294727/1‑235 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
GTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  NC_000913/4579424‑4579889

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