Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,579,656 A→C 24.2% 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,6560AC24.2% 76.6 / 2.6 48V62G (GTG→GGG) symEtoxic peptide regulated by antisense sRNA symR
Reads supporting (aligned to +/- strand):  ref base A (33/2);  new base C (5/8);  total (38/10)
Fisher's exact test for biased strand distribution p-value = 1.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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

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: