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

GTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  NC_000913/4579424‑4579889
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
gtaattttttgaggggggggAAATTTTTTTGTGAGTCGTGGGAAGTCCCCCGGGGGTTAACCGGGAGTAACCCGTTTTACGCTACTTTCTGTTTACCGGCAATCCCCCCAAA‑AACCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCGGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTg                                                                                                                                                                                                                                        <  1:358765/224‑1 (MQ=255)
gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGCC‑AATACAGCCTTCCATCACTTTg                                                                                                                                                                                                                                        >  2:610974/1‑239 (MQ=255)
gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAAGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGGAGGGAC‑AAAAACGCCTTCCACAACtttt                                                                                                                                                                                                                                        >  1:730179/1‑238 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCCGTTCGCTCTCCTCGGCGGAGGGTGGGTGGGCGGGGGGGCC‑AATACAGCCCTCCATCACTTTG‑AAATCTAccc                                                                                                                                                                                                                             >  1:618761/1‑238 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCCATTTG‑ACAttaacc                                                                                                                                                                                                                              >  2:180033/1‑233 (MQ=255)
          aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGGTGGGCGGGTAGGAC‑AATAAAGCCTTCCATCACCTTG‑GCATCTAcc                                                                                                                                                                                                                              >  1:69914/1‑238 (MQ=255)
             gAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCCCTTTG‑ACATCTCccgccg                                                                                                                                                                                                                          >  2:374625/1‑239 (MQ=255)
               gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACAACAACCGCCGt                                                                                                                                                                                                                         >  2:874211/1‑238 (MQ=255)
               gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACACTTTGCACACCTGGCGTAGCGACTGCCTCAGTTCGCTCCCCTCGGGGGCGGGGGGTTGGGCGGTAAGGCC‑AACACAGCCTTCCACCCACTTG‑AAATCTAcccccgg                                                                                                                                                                                                                         >  1:526232/1‑237 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGGGGGGGGGTGGGGGGGAGGGCC‑AATAAAGCCTTCCCCCCCTTTG‑ACATCACCCGCCGTg                                                                                                                                                                                                                        >  1:467750/1‑238 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGGC‑AATACAGCCTTCCATCACTTTG‑CCATCTACCGCCGTg                                                                                                                                                                                                                        >  2:212738/1‑238 (MQ=255)
                ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTg                                                                                                                                                                                                                        >  2:710451/1‑238 (MQ=255)
                ggggAAAGTTTCTACTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTg                                                                                                                                                                                                                        >  2:750574/1‑238 (MQ=255)
                          tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGCC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCaaa                                                                                                                                                                                                             >  2:570379/1‑239 (MQ=255)
                                  gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGGGGCAAAACcggcg                                                                                                                                                                                                      >  2:279026/1‑238 (MQ=255)
                                  gTCAGTGGAAGGCTACCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCATT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGCCAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGCGGTTGGGCGGCGCGTCC‑AATACAGCCTTCCTTCACTTTG‑ACATCTACCGCCGTGCCAGTAGCAAAACcggcg                                                                                                                                                                                                      >  2:125741/1‑238 (MQ=255)
                                                cccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCCCCGTGCCAGTGGCAAAACCGGCGGCTTCCAACCCCTg                                                                                                                                                                                        >  2:460628/1‑238 (MQ=255)
                                                 ccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCCTCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCATCCCCTGa                                                                                                                                                                                       >  2:731156/1‑238 (MQ=255)
                                                            ccGGGAGTAACTCCGTTTACGCGTCTTTTGTTTTACCGCAAATCACCCAAAA‑AAACCCCTGCCCCTGTTTTTGTTTCCGCCCCGGCAGTTGGCACACTGGGGGGAGCGACTGAACCAGTCCGCTCTCCCCGGCGGCGGGTGGTGGGGCGGTGAGGAA‑AATCCAGCCTCCCATCACTTTG‑CCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTg                                                                                                                                                                           <  1:523816/239‑1 (MQ=255)
                                                                 aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCAAGCCACTGACCTTTCAGGGTGATGGc                                                                                                                                                                      >  2:233168/1‑239 (MQ=255)
                                                                 aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGCGTGGGCGGGACGGAC‑AATCCAGCCTTCCATCACTTTG‑TCATCTACCGCCGTGCCCGGGGCAAAACCGACGGCTTCCACCCAATGACCCTTCAGGGTAATgac                                                                                                                                                                      >  1:176042/1‑237 (MQ=255)
                                                                        cctttttccctttttttttttctccGTCAATC‑CTCCAAAAAAACGCCCGCCCCTTTTTTTGGTTCCGCGCCGACTGTTGGCCCCCCGGGCGTAGCGTCAGAATCAGTTCGCTCCCCTCGGCGGCGGGGGGTGGGGGGGGGAGGAA‑AATACAGCCTTCCATCCCTTTG‑CCATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATa                                                                                                                                                               <  1:460628/216‑1 (MQ=255)
                                                                          tGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑CCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCACCCACTTACCTTTCCGGGGGATgcgggggaaaag                                                                                                                                                             >  1:802403/1‑228 (MQ=255)
                                                                                      ttCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCCCTGACCTTTCAGGGGGATGGCGGGAATACG‑GCTGTAATCCgg                                                                                                                                                >  2:445083/1‑239 (MQ=255)
                                                                                              aCCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGGATACG‑GCTGGAATCCGGGTGGCgg                                                                                                                                         >  2:452424/1‑237 (MQ=255)
                                                                                              aCCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGAGGGCGGGAAAACGCGCTGAAAACCGGGGAgga                                                                                                                                          >  1:422850/1‑236 (MQ=255)
                                                                                              aCCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCACCCACTGCCCTTTCCGGGGGCGGGCGGGGAAACG‑GCCGTAATCCGGGGGGCGa                                                                                                                                         >  1:880184/1‑238 (MQ=255)
                                                                                              aCCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGATTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACAGCATCAGTTCGCCCTCCTCGGCGTCGGGGGGTTGGGCGGTGAGGGC‑AATGCAGCATTCCATCCCTTTG‑ACATCTACCGCCGTGCCAGTGGCGAAACAGGAGGCTTCCAGCGCATGACCTTTCCGGGGGAGGGCGGGATAACG‑GCTGGAATCCgggggggca                                                                                                                                         >  1:662247/1‑232 (MQ=255)
                                                                                                   cAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGCGGGGAATACG‑GCTGTAAACCCGGGAGCGAccccca                                                                                                                                   >  1:332963/1‑234 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGCCCTTTCAGGGTGATGGGGGGAATACG‑GCCGTAATCCGGGTAGCGACTCGCaaaaaccgggggggcaag                                                                                                                  >  1:747405/1‑225 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGGTGGCGGGGATAAG‑GCTGTAAACCGGGCAGCCACCCGCCAAAACG‑ACGGGGaaaaa                                                                                                                 >  1:102816/1‑235 (MQ=255)
                                                                                                                     gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGGGaaccccaaaaccaccgggggcctt                                                                                                                  >  2:715111/1‑216 (MQ=255)
                                                                                                                         gCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGGAATCCGGGTAGCGACTCGCATAACCA‑ACGGTGACATGACgg                                                                                                             >  2:186802/1‑239 (MQ=255)
                                                                                                                         gCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGAGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGGATACG‑GCTGGAATCCCGGGAGGGACCCGCCTaaaccgcgggggaatcaggg                                                                                                              >  1:602589/1‑221 (MQ=255)
                                                                                                                          cACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCAAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGAACCCCTTAACCG‑ACGGTGACATGACGGt                                                                                                            >  2:125141/1‑239 (MQ=255)
                                                                                                                           cccggCTTTTGTTTACGGCCCGAAAGTTTGCCACCCGGGCGAGGGCAAGGCATCAGTTCGCTCTCCGCGGGGGGGGGTGGTTGGGGGGTGAGGAA‑AATACAGCTTTCCATCACTTTG‑ACATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGtt                                                                                                           <  2:571765/235‑1 (MQ=255)
                                                                                                                              tGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGCCGGTTAtt                                                                                                        >  2:416819/1‑239 (MQ=255)
                                                                                                                                    tGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCgg                                                                                                  >  2:130179/1‑239 (MQ=255)
                                                                                                                                                       tGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATGCAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGACTACG‑GCTGTAATCCGGGTAACGACTCGCATAAACGAAGGGGGACGAGA‑GGTTTTTTGCCGGGGGAACTTCCGGTTCGAAc                                                                               >  1:336741/1‑239 (MQ=255)
                                                                                                                                                               tGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACTTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGGGACATGACGGTTATTTGCCGGGGGGACTTCTGCTTCGAACGGtttggc                                                                       >  2:509635/1‑235 (MQ=255)
                                                                                                                                                                       gagcggCATCAGTTCCGTCCCCCGGGCGGGGGGTGGTTGGGCGGTGAGGAA‑AATCCAGCCTCCCATCCCTTGG‑ACATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAAt                                                               <  2:102816/234‑1 (MQ=255)
                                                                                                                                                                            gcagtaGTTCGCTCGCGTCGGCGGGGGGTGGTGGGGCGGTGATGAC‑AATCCAGCCTTCCATCCTGTTG‑ACATCTACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCgtg                                                          <  2:440662/234‑1 (MQ=255)
                                                                                                                                                                                  gTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCACGTGGCAAAACCGGCGGCTTCCAGCCCCTGCCCTTTCAGGGTGATGGCGGGAATACG‑GCCGTAATCCGGGTCGCGACTCCCCTAACCC‑CCCGTGCCATACCAGTTATTTTCCGGGGGGACGTCTGCGTCGCACGGGTGTGCAACAGAAAGCgtgtgcggg                                                    >  1:430592/1‑234 (MQ=255)
                                                                                                                                                                                            ctggggggggggTTGTTGGGGGGTAGGGCA‑AATCCACCCTTCCACCCCTTGG‑ACTGCTCCCGCCGGGGCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTa                                          <  2:430592/230‑1 (MQ=255)
                                                                                                                                                                                                 gcgggggTTGTGGGGGGGTG‑GGACAAATCCCGCCTTCCA‑CCCTTTGCACTTCACCCGCCGGGCAAGGGGCAAAACCGGCGGCTGCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATtctc                                     <  1:710451/235‑1 (MQ=255)
                                                                                                                                                                                                  ggcggGTGGTGGGGGGGTGAGGAC‑ATTACAGCCTTCCACCCTTTGG‑ACATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTGCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATtctc                                     <  2:764659/238‑1 (MQ=255)
                                                                                                                                                                                                      ggtgttTGGGCGGTGAGGAC‑AATCAACCCTCCCATCCCTTTG‑CCATCTACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAgg                                 <  2:400553/238‑1 (MQ=255)
                                                                                                                                                                                                      ggtggtTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTTTGCAATAGCATGCGTGTCAGTCATAACTGCTATTTTCCAGGa                                >  2:464561/1‑239 (MQ=255)
                                                                                                                                                                                                                ggggaaGGAA‑AATCCAGCCTCCCATCCTTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgatt                       <  2:513305/233‑1 (MQ=255)
                                                                                                                                                                                                                 ggggtaggaaaatcacgCCTTCCATCTCTTTG‑CCATGTGCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgatt                       <  2:618761/222‑1 (MQ=255)
                                                                                                                                                                                                                     aGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCCTAACCG‑ACGGGGACATGCCGGTTATTTGCCGGGGGGACTTCTGCTTTGAAAGGTTGGGCAATAGAATGCGTGTCAGTCATAACAGCCATTCCCCAGGAAAAgggattgggatt                 >  1:695099/1‑239 (MQ=255)
                                                                                                                                                                                                                         aaaaatcAGCCTTCCACCCATTGG‑AACTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAg               <  2:832193/232‑1 (MQ=255)
                                                                                                                                                                                                                                  ccTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGgt     >  2:487995/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGGGACATGACGGTTTTTTGCCGGGGAGACTTTTGCTTCGAACGGTTGTGCAAAAGAATACGGGTCCGGCATAAATTCCATTATCCAGGAAAAGTGATTGAGAATAACGAAGCgggggga  >  1:294727/1‑235 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
GTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAAT‑AAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAC‑AATACAGCCTTCCATCACTTTG‑ACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACG‑GCTGTAATCCGGGTAGCGACTCGCATAACCG‑ACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  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: ‑