Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,343 T→C T70A (ACG→GCG)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,3430TC97.7% 135.4 / ‑4.2 44T70A (ACG→GCG) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (19/24);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

ACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCT  >  NC_000913/66102‑66581
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAAc                                                                                                                                                                                                                                        >  1:385792/1‑251 (MQ=255)
                               caccGTTGATTGCTGCGTCGTTCATTTGGCGGGGGCAGGGAATGGTTCCGTAGAAATAGTCGGCGGTGGG‑GGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGCCGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGAc                                                                                                                                                                                                         <  2:170752/248‑1 (MQ=255)
                                accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCGCCAATGGACGGGAATGCCTGGTAGAGCAGTCGGTGAGGTGGCGTATCGGGGGAGGGCTTTTTAGCACCTTCAACCCCTT‑CACCGGTTTCGGTGCTAaccaccaccc                                                                                                                                                                                                       >  2:90375/2‑250 (MQ=255)
                                accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGAGCAGGCGGTGAGTTGGCGGATCGGGGGGGGGCTTTTTCGCACCTTCAACCACTT‑CCCCGGTTTCGATGCTAACCACGAcc                                                                                                                                                                                                        >  2:420630/2‑250 (MQ=255)
                                 ccGTTGGTTTCTGCGTCAGTCATTTTGCGGGTGCAGGGGATGATGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGTCCGCCTGCGCCCAGATAGAGGCGTGGCGCGCGGGTGTGTGGACCATGCCGCCAATGGACGGGAATGCCTGATAGAGGAGTCGGTGAGGTGGCCTATCGGAGGAGGGGTTTTTCGCACCTTCAACCACTTTCACGGGTTT‑GATGcgaaacccggccc                                                                                                                                                                                                       >  2:80916/1‑238 (MQ=255)
                                 ccGTTGATTTCTGGGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCTTAGAAAGAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAAGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCat                                                                                                                                                                                                      <  1:289008/251‑1 (MQ=255)
                                                      aTTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGCTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCt                                                                                                                                                                                 >  1:157203/1‑251 (MQ=255)
                                                                  cAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGcc                                                                                                                                                                     <  1:124492/251‑1 (MQ=255)
                                                                             gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                                                                                                          <  1:198964/251‑1 (MQ=255)
                                                                             gCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGa                                                                                                                                                          <  1:413474/251‑1 (MQ=255)
                                                                                   gAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATTGACGCCGGAAGGTTTGATcaca                                                                                                                                                     >  1:133066/1‑250 (MQ=255)
                                                                                    aaaTAGTCGGCT‑TGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACaaa                                                                                                                                                   <  1:379824/251‑1 (MQ=255)
                                                                                               g‑tgggtggTGCCGGTTGTTGGATGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGTCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCCCCTTCAACCACTT‑CACCGGTTTCTATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCt                                                                                                                                        <  1:423751/251‑1 (MQ=255)
                                                                                                        tgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAAGCCGCGCTCGCGATc                                                                                                                                 >  1:293713/1‑250 (MQ=255)
                                                                                                                   gaaaccgcttggcccGCTTGCGCCCGGATAGTGGCGTGGCGTGAGTGTGTGTGCACAATGCCGCCAATGGACGGTAATGCCTGATAGAGCAGTCGGTGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCAGTTTCAACCACTT‑CACCGTGTTCGATGCTAACCACGTCCATAGCGTCATCGCTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGAc                                                                                                                     <  1:427298/239‑1 (MQ=255)
                                                                                                                           ctggcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGcccc                                                                                                             >  2:160563/5‑251 (MQ=255)
                                                                                                                                  gCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAg                                                                                                      <  1:41979/251‑1 (MQ=255)
                                                                                                                                       cgcCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGa                                                                                                 <  1:160563/251‑1 (MQ=255)
                                                                                                                                           cAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATAGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTTCCCCATGTGAGCGTGAcc                                                                                               >  1:88820/1‑249 (MQ=255)
                                                                                                                                              atagtggcgtggcGCGAGTGTTGGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGtt                                                                                          <  1:438626/248‑1 (MQ=255)
                                                                                                                                                      gtggcgCGAGGGGGTGTGCACAAGGCCGACAATGGACGGGAATGACTGATAGAGCAGTCGGTGAGTTGGCGTATAGGAGGAGGGCTTTTGCGCACCTTCAACCACTT‑CACCGGTTGCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGTCCCATGTGAGCGTGACCAGGTTATGTTTTg                                                                                  <  1:185370/251‑1 (MQ=255)
                                                                                                                                                        ggcgCGATTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTgg                                                                                 <  1:114107/250‑1 (MQ=255)
                                                                                                                                                          cgcgAGTGTGTGTGCCCAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTTAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAg                                                                              <  1:2956/251‑1 (MQ=255)
                                                                                                                                                          cgagAGTGTGTGTGCACAATGCCGCCAATGGAAGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAGAGACGCAGCGCTCGCGTTCAACGGCGCTGACGTTGCCCCATGGGAGCGTGACAAGGTTATGTTTTGGc                                                                                >  2:180700/4‑249 (MQ=255)
                                                                                                                                                                 ttttgtgctcactGCCGCCAATGGTCGGGAACGCCTGATAATGCAGGCGGTTATTTGGCGTATCGGGGGAGGGCTTTGTCGCACCTTCAACCTCTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATAGACGCCGGAAGGGTTGATCACAAAGACGCCGCGCTCGCGATCAACTGCGCTGCCGTTGCCCCATGTGAGCGTGACCAGGTTATTTTTCGGCAGCGCCAg                                                                        <  2:251870/238‑1 (MQ=255)
                                                                                                                                                                        cacaATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGcct                                                                >  1:255199/1‑249 (MQ=255)
                                                                                                                                                                         acaATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGcctc                                                               >  1:110375/1‑248 (MQ=255)
                                                                                                                                                                              gccgccAATGGACGGGAATGCCTGATAGAGCAGTAGGTGAGTTTGCGGATCGGAGGAGGGCGTTTTCGCACCTTCAACCACTT‑CACCGGGTGCGAGGCTAACCGCGACCATATCGTCATCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACTAGGTTATGTTTTGGCAGCGCCATGTGGGCTTCTAAtt                                                          <  2:426025/251‑2 (MQ=255)
                                                                                                                                                                                  ccAAAGGACGGGAATGCATGATTGAGCAGTCTGTGAGGTGGCGTATCGGGTGTGGGCGTTTTCGCACCTTCAACCACGT‑CACCGGGTTCGATGCTAACCACGACCATCGCGTCAGCGGTCATGATTCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGTGCCGACGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCTCCAGGTTGGCCTCTAATatct                                                       <  2:235404/250‑4 (MQ=255)
                                                                                                                                                                                    aaTGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCGCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGc                                                    >  1:109344/1‑251 (MQ=255)
                                                                                                                                                                                      tGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGt                                                  >  1:463043/1‑251 (MQ=255)
                                                                                                                                                                                         acgggAATGCCTGATAGAGCCGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGCCCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTCATACCTGGCGTTTg                                               >  2:445770/3‑251 (MQ=255)
                                                                                                                                                                                         acgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                               <  2:310626/249‑1 (MQ=255)
                                                                                                                                                                                           gggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                               >  1:282904/1‑249 (MQ=255)
                                                                                                                                                                                                     gATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGGCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCGTGTGAGCGTGACCAGGTTAGGTTGTGGCAGCGCCAGGTTGGCTTCTAATACCTGGTGTTTGAGATCTTCTa                                     <  1:407734/249‑1 (MQ=255)
                                                                                                                                                                                                     gATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCATGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTaa                                    <  2:418595/250‑1 (MQ=255)
                                                                                                                                                                                                        agagCAGTCGGTGAGTTGGCGTATCGGAGTAGGGCTCTTTCGCACCTTCAACCACTT‑CCCCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGCCGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCATCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACAt                                 <  2:129173/250‑1 (MQ=255)
                                                                                                                                                                                                          agcagtCGGTGAGGTGGCGTATCGGAGGCGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGACAGCGGTCATGATGCTGTACTCGACGCCGGAAGGTTTGATAACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTGATAACTGACGTTTGAGATCTTCTAAAATGtt                              >  2:399782/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggTGAGTTGGCGTAGCGGAGGAGTGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCt                       <  2:199514/251‑1 (MQ=255)
                                                                                                                                                                                                                     aGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAATGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATCCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGt                   >  1:244836/1‑251 (MQ=255)
                                                                                                                                                                                                                             tatCGGAGGCGGGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTACCCACGACAATATCGTCAGCGGTCATGCTGCTGTAATCGTCGCAGTAAGGTTTGATCACAAAGTAGACGCGCTCGCGATCAACGGCACTGACGTTGCCCCATGCGAGCGCGACCAGGTTATATTTTGGCAGCGCCAGGTTGACCGCTAATACCTTGGTTTTGGGATCTTCTACAATGTTAACTACTGCGTGACGGagt           >  2:334826/3‑249 (MQ=255)
                                                                                                                                                                                                                                  gaggagGTCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATAGCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTGGCCCCATGTGAGCGTGACCAGTTTATGTTTTGGCATCGCCAGGTTGGCCTCTAATGCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCtt      <  1:440028/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gaggGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCCCGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGCCGTTGCCCCATGTGAGCGAGACCAGGTTATGTCTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTATCTCCTTCTTGCCGGTTGCGCTTTGc   >  2:282039/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aggGCTTTTTCGCACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCt  <  1:400719/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
ACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCG‑TGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTT‑CACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCT  >  NC_000913/66102‑66581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑