Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,092,703 A→G F41L (TTC→CTC)  frvA ← putative PTS enzyme IIA component FrvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7030AG98.7% 259.2 / NA 75F41L (TTC→CTC) frvAputative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/43);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATG  >  NC_000913/4092472‑4092949
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
acaGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCCGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGGTTGAGTCGGGGAt                                                                                                                                                                                                                                      >  2:700254/1‑250 (MQ=255)
  tgtgtGCCGCTGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGTCGTCACCGGCACCAAATCCAGTAGAGGGCATTTTTTCGAGCAGCAGAAGGGTTTGCAGGGGCTGGTGTGAGTCGGTGGTaa                                                                                                                                                                                                                                    >  2:726620/2‑250 (MQ=255)
   gtttgCCGATGAATTTGACCTTAGCCTTATCGCCGCTCTGCGGCACGCCGAGGCAGATCCAGCAATTGACGGCTTCGCCATCGCTTGCTTCCCAGTCAGTAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTGTAACGCAGGCGCGTTTACCGGGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGGGCAGTTATTCGCGCAGGAGCAGGGGTTGCAGGAGCTTGTGTGAGGCGTTGATaaaa                                                                                                                                                                                                                                  <  2:425325/248‑1 (MQ=255)
         cGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGttt                                                                                                                                                                                                                            <  1:16144/251‑1 (MQ=255)
           atgatTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACTCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGATTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTg                                                                                                                                                                                                                          >  1:418210/1‑251 (MQ=255)
                tttGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACg                                                                                                                                                                                                                     >  1:477870/1‑251 (MQ=255)
                    acctgtgccTCTCCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCGTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACTGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTAt                                                                                                                                                                                                                 <  1:41144/244‑1 (MQ=255)
                     ccTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTAtt                                                                                                                                                                                                                <  1:363554/251‑1 (MQ=255)
                             tctTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTAACGTGCGGCACGGCGACACCGGAACCACATCCAGTCGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGCGCTGGTGTGAGTCGGTGCTAAACACGTTTTGTAACGCTATTGTCGCCa                                                                                                                                                                                                         >  1:394547/1‑250 (MQ=255)
                                ccgccgcTTTGCGGCACGCCGAGGCAGGTCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAAGAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTg                                                                                                                                                                                                     <  2:431663/250‑1 (MQ=255)
                                ccgccgcTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACt                                                                                                                                                                                                      <  2:760101/249‑1 (MQ=255)
                                 cgccccgttccGGCCCGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAAGAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTTCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGATTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGc                                                                                                                                                                                                    <  1:601172/241‑1 (MQ=255)
                                  gccgcTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCt                                                                                                                                                                                                   <  1:653015/251‑1 (MQ=255)
                                         gtcGGCACGCCGAGGCAGATCCAACAATTGACGTCTGCGCCATCGCTGGCTGTTCAGTCAATAGCCTGCGCTTTTCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGATCCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAAc                                                                                                                                                                                            <  2:70605/249‑1 (MQ=255)
                                           cGGCACGCTGAGGCGGATCCAGCTATTGACGTCTTAGCCAGCGCGGGCGTTTCAGTCAATAGCCTTCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGTAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACgg                                                                                                                                                                                          <  2:448682/251‑1 (MQ=255)
                                             gctcGCCGAGGCCGATCCAGCAAGTGTCGTCGTCGCCATCGCTGGCTTGCCAGGCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGATTAACGCAGTCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTGTTCGCGCAGCAGCAGTGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACgtaa                                                                                                                                                                                        <  2:155349/248‑4 (MQ=255)
                                                 gCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATa                                                                                                                                                                                      >  1:24089/1‑249 (MQ=255)
                                                  ccGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGGTAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAgc                                                                                                                                                                                   >  2:344675/1‑251 (MQ=255)
                                                   cGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCAAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAgcg                                                                                                                                                                                  >  2:553968/1‑251 (MQ=255)
                                                    gAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCGTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGACCGGCGCCACCGGAACCACATCCAGAAGAGTGCCTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCAAACTGCTTCAGAACGGAAtcaacgc                                                                                                                                                                                 >  2:56942/1‑245 (MQ=255)
                                                      ggCAGATCCAGCAATTGACGTCTTCGCCAGCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCa                                                                                                                                                                               <  1:379855/251‑1 (MQ=255)
                                                      ggCAGATCCAGCAATGGACGTCTTCGTCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCAGTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAgcgaca                                                                                                                                                                               <  1:379854/251‑4 (MQ=255)
                                                               aGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGCGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                                       >  1:789199/1‑250 (MQ=255)
                                                                     tGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTc                                                                                                                                                                >  2:269896/1‑251 (MQ=255)
                                                                               gCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATa                                                                                                                                                       >  1:796985/1‑250 (MQ=255)
                                                                                  agcGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGc                                                                                                                                                   <  1:167798/249‑1 (MQ=255)
                                                                                  aTCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGACTAAGCACAATTGCCCTGAATATTCAGGTCAATACatc                                                                                                                                                   >  2:200079/1‑249 (MQ=255)
                                                                                     gCTGGCATTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGGTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTg                                                                                                                                                <  2:210380/251‑1 (MQ=255)
                                                                                        ggCTTTACAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAGCGGTTGTTTAACGCAGGCGCTTTTCACGTGAGGCACGGCGACACCGGTACCAAAGCCAGTAGATGGCATTTTTTCGCGCAGCAGCGGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGTTAGTGCGGCCATCTGGTTGAGAACAGAATAAGCGCCATTGTCCTGAATATTCTGGTCAATGCAGCTTGc                                                                                                                                               <  2:741578/249‑1 (MQ=255)
                                                                                         gCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGt                                                                                                                                            >  1:260052/1‑251 (MQ=255)
                                                                                                   cAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTTGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGc                                                                                                                                   <  1:543665/250‑1 (MQ=255)
                                                                                                   cAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcc                                                                                                                                  <  2:239216/251‑1 (MQ=255)
                                                                                                    aaTATCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCa                                                                                                                                 <  2:523539/251‑1 (MQ=255)
                                                                                                     aTAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGACCTGCCAt                                                                                                                                >  2:388135/1‑251 (MQ=255)
                                                                                                       aGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTa                                                                                                                              <  1:168112/251‑1 (MQ=255)
                                                                                                       aGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTGTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAAGAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAtt                                                                                                                               <  1:387706/250‑1 (MQ=255)
                                                                                                            gcgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                                         <  2:336564/251‑1 (MQ=255)
                                                                                                            gcgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGGCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATc                                                                                                                           <  2:477870/249‑1 (MQ=255)
                                                                                                            gcgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCACTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCTATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                                          <  1:824824/250‑1 (MQ=255)
                                                                                                             cgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGa                                                                                                                        <  1:600119/251‑1 (MQ=255)
                                                                                                                  tGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTcc                                                                                                                    <  1:269896/250‑1 (MQ=255)
                                                                                                                      ggCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaa                                                                                                               <  1:465010/251‑1 (MQ=255)
                                                                                                                         gAATAATACGAACGGTTGTTGAACGCTGGCGCGTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGGTCAGTTTGTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGGCGGTGATAGAGCCGATGTGTAACGCTAGTGGCGCCAACTGCTGCAGAACGGAATAAGGGCCGGTGCCCTGAATATTCGGGGCAATACAGGTGGGAGGAAGAGGTGCGAGTAATCGGACCCCATaaaaga                                                                                                            <  2:384698/251‑3 (MQ=255)
                                                                                                                           ataataCGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatata                                                                                                          <  1:344675/251‑1 (MQ=255)
                                                                                                                              ataCGAACGGTTGTTTAAAGCAGGCTCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATGTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATAttt                                                                                                       <  2:137639/251‑1 (MQ=255)
                                                                                                                               taCGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTCGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATACGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTTCAGTAAGAGCTGCCATTAATCGGACTCCATAAACTATCtttt                                                                                                      >  2:158061/1‑251 (MQ=255)
                                                                                                                                aCGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCCTAAAATATAtttt                                                                                                      >  2:471572/1‑250 (MQ=255)
                                                                                                                                   aaCGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATc                                                                                                  <  1:146942/251‑1 (MQ=255)
                                                                                                                                     cGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAAATGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGATTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCa                                                                                                 >  2:51633/1‑250 (MQ=255)
                                                                                                                                       gttgttTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAtt                                                                                              >  1:215221/1‑251 (MQ=255)
                                                                                                                                            ttacagCAGGCGCTTGTACCGTGCGTTACGGCGTCCCCGGAACTAAATCCTGTAGTGTGCATTTTTTCTCGCAGCAGCCTGGTTTGCCGGATCGTGTGTGAGCCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGTCATTGCCCGGAATATTCAGGTCAATACAGCTTGCATTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTaa                                                                                          <  2:282435/245‑1 (MQ=255)
                                                                                                                                                          tttACCGTGCGGCACGGCTACACCGGAACCAAATCCAGGAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTGGCAGGAGCTGGTGTGAGGCGGTGATAAAACCGATGTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCATGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATc                                                                                                  <  2:65109/228‑1 (MQ=255)
                                                                                                                                                          tttACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCATGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATc                                                                                                  >  1:65109/1‑228 (MQ=255)
                                                                                                                                                            tACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAg                                                                          <  1:489405/250‑1 (MQ=255)
                                                                                                                                                             aCCGTGCGGCTCGGCGACACCTGAACCTAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCTCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtg                                                                        <  2:579227/251‑1 (MQ=255)
                                                                                                                                                             aCCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTc                                                                                                                     >  2:373132/1‑206 (MQ=255)
                                                                                                                                                             aCCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTc                                                                                                                     <  1:373132/206‑1 (MQ=255)
                                                                                                                                                                gTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtaatg                                                                     >  2:650178/1‑251 (MQ=255)
                                                                                                                                                                         ggCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCCGTAAGAGCTGCCATTCATCGGACTCCATAAAATCTATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGat                                                            >  1:521778/1‑251 (MQ=255)
                                                                                                                                                                           cGACGCCGGAACCAAATTCAGTTGAGTGCATTTGTTCGCGCGTCAGCAGGGTTTGCAGGAGTTGGTGTGTGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCACCTGCTTCAGAACGTAAGAAGCGCCATTGCTCTGAATATTCAGGTGAATACAGCTTGCAGTAAGAGCTGTCATTAATCGGACTCCATAAAATATATGTTAATCAATTGTAATGGTGGATCGATTAAGTGATGAGAGGCGatat                                                          <  2:834823/251‑1 (MQ=255)
                                                                                                                                                                                       cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGACTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                                         >  2:845410/1‑176 (MQ=255)
                                                                                                                                                                                       cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                                         <  1:845410/176‑1 (MQ=255)
                                                                                                                                                                                         aaTCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                                             >  1:647667/1‑250 (MQ=255)
                                                                                                                                                                                                tAGAGTGGATTTTTTCTCGCAGCATACGGGTATGCAGGAACTGGAGTGAGTCGGTGATAAAACCGTTTTGTAACGATATTGTCTCCAAATGCTTCAGAACGGAATAAGCGCCATTGCCATGACTCTTCAGATCAATACAGCTGGCAGTAAGAGCTGCCATTAATCGGACTCCATCAAATATATTTTAATCAATTGTATTGGAGGAGCGCTTCGGTGCTAAGAGGAGATCTGTCTGTAAACAATGAAAAttt                                     >  2:81693/1‑251 (MQ=255)
                                                                                                                                                                                                     tGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTACACAATGACAATTTGAGAc                                >  2:508640/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc                                <  1:650178/250‑1 (MQ=255)
                                                                                                                                                                                                                gcgcAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAgg                     >  2:685611/1‑251 (MQ=255)
                                                                                                                                                                                                                  gcagcagcaGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGTCGATATGTCGGTAATCAATGACAATTTGAGACAGAGTGAAAGGTc                   <  1:246507/251‑1 (MQ=255)
                                                                                                                                                                                                                    agcagcagGGTTTGCAGGAGCTGGAGTGCGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGTCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGACGGAAAACAATGACAATTTGAGCCAGAGTGAAAGGTc                   >  1:332936/1‑249 (MQ=255)
                                                                                                                                                                                                                      cagcagGGTTTGCAGGGGGTGGTTTGAGTCGGTGATAAAACCGTGTTGTAACGCTATTGTCGCCAACCGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATGCAGGTCAATACAGCTTGCATTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAGTGTAATGGTTGAGCGATTAAGTGATTAGAGGCGATATGTCGGTAAACAATGGCAATTTGAGCCAGATTCAAAGTTCAGAt               <  2:539673/251‑1 (MQ=255)
                                                                                                                                                                                                                      cagcagGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGACCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGAt               <  1:553968/251‑1 (MQ=255)
                                                                                                                                                                                                                          agGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCCTTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGc           >  2:235285/1‑251 (MQ=255)
                                                                                                                                                                                                                                tGCAGGAGCTGGGGTGAGTAGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAAAGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGCGCGATTAAGTGATGAGAGGAGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGt      >  2:9023/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAAt   <  1:200079/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGGCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATg  <  1:543816/251‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATG  >  NC_000913/4092472‑4092949

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