Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,288 C→T D141N (GAC→AAC)  fsaB ← fructose‑6‑phosphate aldolase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,2880CT98.7% 268.4 / NA 74D141N (GAC→AAC) fsaBfructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/47);  minor base A (1/0);  total (27/47)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTG  >  NC_000913/4139047‑4139534
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
tAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCCGGGGCGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGGGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGAACGAATGCCGGTTCCGCCCTg                                                                                                                                                                                                                                               >  1:497730/1‑251 (MQ=255)
tAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCCTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTg                                                                                                                                                                                                                                               >  2:359903/1‑251 (MQ=255)
       gAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCt                                                                                                                                                                                                                                        <  1:615380/251‑1 (MQ=255)
        agtagtGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTa                                                                                                                                                                                                                                       <  2:146264/251‑1 (MQ=255)
         gtagtGCTAAATGCGGCATACCAGTCGTGTTCGAACGTCTCTATAGCTGACTCTACCGCAGGGGTGATGAGCATTTGTTGCGCTACAACTAAGGGAAGGGTGATGGAGTCACATCCGGCCAGTAAACAGTCCAGCGCCTGACGCGTCGTTTTAAAGATGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTAGGCATCTAc                                                                                                                                                                                                                                      <  2:240519/251‑1 (MQ=255)
              gCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCTCCCTGGGCATCTACGCGGt                                                                                                                                                                                                                                 >  2:266032/1‑251 (MQ=255)
               ccAAATGCGGCATTCCCGTCGTGTTCGCATTTCTCTATAGCGGACGCTACCGCGGGGGTGTTGAGCATTTGTTTCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGtt                                                                                                                                                                                                                                <  2:63728/251‑1 (MQ=255)
                 aaaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTaa                                                                                                                                                                                                                              <  1:316855/251‑1 (MQ=255)
                 aaaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTa                                                                                                                                                                                                                               >  1:253185/1‑250 (MQ=255)
                   aTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAAc                                                                                                                                                                                                                             >  1:399321/1‑250 (MQ=255)
                     gCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACATTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACAt                                                                                                                                                                                                                           <  1:337928/250‑1 (MQ=255)
                     gCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACAt                                                                                                                                                                                                                           <  1:449537/250‑1 (MQ=255)
                        gCATTCCATTCGTGTTCGAACTCCTCGATACCTGACGCTTCTGCAGGGGTGTTGAGCATTTGGTGCGCTAAATCTAAGGGAAGGGGGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTTGATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACgg                                                                                                                                                                                                                       <  2:333728/251‑1 (MQ=255)
                         catgccAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCGAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGa                                                                                                                                                                                                                      <  1:125538/247‑1 (MQ=255)
                         cATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGa                                                                                                                                                                                                                      <  1:331233/251‑1 (MQ=255)
                         cATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGa                                                                                                                                                                                                                      <  1:209033/251‑1 (MQ=255)
                              cAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAAc                                                                                                                                                                                                                 >  1:583302/1‑251 (MQ=255)
                                   gtgtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtt                                                                                                                                                                                                            <  2:455517/251‑1 (MQ=255)
                                     gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                          <  1:597854/251‑1 (MQ=255)
                                                               aCCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCtaataa                                                                                                                                                                                   <  1:41538/248‑1 (MQ=255)
                                                                ccGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCt                                                                                                                                                                               >  2:378922/1‑251 (MQ=255)
                                                                 cGCAGGGGTGTTGAGCATTTGTTGCGCTCCATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGTAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCtt                                                                                                                                                                              <  1:90883/251‑1 (MQ=255)
                                                                 cGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCtt                                                                                                                                                                              >  2:166587/1‑251 (MQ=255)
                                                                 cGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCtt                                                                                                                                                                              <  1:35393/251‑1 (MQ=255)
                                                                          gtTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTgcggcg                                                                                                                                                                      <  1:172290/250‑1 (MQ=255)
                                                                             gAGCATTTGTTGCCCTACTGGTAGGGGAAGGGGGATGGATTCACATCCGGCCAGTAAACAGTCCAGCGCCTGGCGCGTCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTTATAACCCTTGTGCGGCGCtat                                                                                                                                                                  <  2:454367/251‑1 (MQ=255)
                                                                                     tttgCGCCAAATCTAAGGGCAGGGTGATGGTTTCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGGTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagc                                                                                                                                                           <  2:587181/249‑1 (MQ=255)
                                                                                        gcgcTACCTCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAttt                                                                                                                                                                                                           >  1:298141/1‑199 (MQ=255)
                                                                                        gcgcTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAttt                                                                                                                                                                                                           <  2:298141/199‑1 (MQ=255)
                                                                                         cgcTACATCTAAGGGCAGGGGGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGc                                                                                                                                                      >  2:567583/1‑251 (MQ=255)
                                                                                             aCATCGAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCa                                                                                                                                                    <  1:331114/249‑1 (MQ=255)
                                                                                               aTCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGACCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACACTTGTGCGGCGCTATATACAGCAGTGCCACgtg                                                                                                                                                >  2:381936/1‑251 (MQ=255)
                                                                                                   aaGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtag                                                                                                                                             <  1:429785/250‑1 (MQ=255)
                                                                                                            gTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCctc                                                                                                                                   <  1:376941/251‑1 (MQ=255)
                                                                                                               atggctttacaTCCTGCTAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCAGTTCTAACAGCGTTTGCAGCTCCTGAGCCGTACGAATGCCGTTTCCTCTCTGGGCATCTACGCGGTTAACATACGGAGCAGCGTCTGTTGCCCCTGCCAGTTCGGCTAATAACCCTTGTGCGGCGCTAGATACAGCAGTGCCAAGTGTAGTAATACCCTCttt                                                                                                                                <  2:569772/243‑1 (MQ=255)
                                                                                                                 ggATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCttttt                                                                                                                              <  1:359903/251‑1 (MQ=255)
                                                                                                                     tctcaTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCa                                                                                                                           <  2:322907/247‑1 (MQ=255)
                                                                                                                           ccTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGTGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtttt                                                                                                                    <  2:253185/251‑1 (MQ=255)
                                                                                                                              gCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTa                                                                                                                   <  1:385718/249‑1 (MQ=255)
                                                                                                                              gCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAAt                                                                                                                 <  2:242766/251‑1 (MQ=255)
                                                                                                                                cAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtg                                                                                                               <  1:167578/251‑1 (MQ=255)
                                                                                                                                  gTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCCTGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgct                                                                                                             <  2:289718/251‑1 (MQ=255)
                                                                                                                                  gTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgct                                                                                                             <  1:191096/251‑1 (MQ=255)
                                                                                                                                    aaaCAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTGACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCGTTTGTCAGTATGTTAATtgctg                                                                                                            <  2:159412/250‑1 (MQ=255)
                                                                                                                                           ccatggcCTGACGCCGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAGCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTGCAGTATTTTAATTGCTGCCAGAcc                                                                                                     <  2:452863/245‑1 (MQ=255)
                                                                                                                                                gcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCgg                                                                                                >  1:390912/1‑250 (MQ=255)
                                                                                                                                                         ggcgTGTTAAAGCGGGCTGCCAGCACCATGCTTTCTGGCGCGTGAATTTATAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAAAGCAGTGCCAAGTGTAGAAATACCCTATTTTTTCAGCATTTTAATTGCTGCCAGACCTTCGGAAGTCACCgg                                                                                      >  2:382232/1‑251 (MQ=255)
                                                                                                                                                              tttAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGAttt                                                                                 >  2:208383/1‑251 (MQ=255)
                                                                                                                                                                        gctgCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAAt                                                                       >  1:32893/1‑251 (MQ=255)
                                                                                                                                                                             cAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCg                                                                  <  2:243670/251‑1 (MQ=255)
                                                                                                                                                                                 aCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                              >  2:274725/1‑251 (MQ=255)
                                                                                                                                                                                 aCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                              <  2:276690/251‑1 (MQ=255)
                                                                                                                                                                                      gCTTTATGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGa                                                                                                       >  1:425179/1‑205 (MQ=255)
                                                                                                                                                                                         ttCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTc                                                        <  1:266032/249‑1 (MQ=255)
                                                                                                                                                                                             gtcgcgTTCATTTCTAAAGGCGTTTCCCGCTCCTGCACCGTAATAATTCCTTTGCCGCCCTTGGCATCTGCGCGGGTATCATACGGAGCAACGTATTTTTCCCCTGCCAATGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTATAACAACAATTCCCGGAATAGCGTCTCGCa                                                   <  2:592517/248‑1 (MQ=255)
                                                                                                                                                                                             ggCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCGTCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCa                                                   <  1:519490/250‑1 (MQ=255)
                                                                                                                                                                                                          cTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttct                                     >  1:630616/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCAGCGTTTGTAGAACCTGAACAGTACGAAAGACGTATCCGCCCAGACTATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTAAGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTCTTCAGAATTTTAATTGCTGCCAGACCTTCGGCAGTCACCGGGCTTTTCACCACAATACCCGGACTAGCGTCACACAAGCGCTTCGCttcttc                                   >  2:540495/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCttcttc                                   <  1:274725/251‑1 (MQ=255)
                                                                                                                                                                                                            aaCAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCttcttc                                   >  1:589093/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTACGCCCTGGGCATATAAGCGGTTAACATACGGAGCAACGAATTTTGCCCCTGCCAGTGCGGCTAATACCCCTTGTGCGGCGCTACATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCACGCTCTTCGGAAGTCACCGGGATTTTCACCACAATACCCAGCATCGCGTCGCGAAGACGCTCTGCttcttc                                   >  2:629817/1‑251 (MQ=255)
                                                                                                                                                                                                               aGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTccac                                <  1:79876/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccgttccAGCTCCTGATCCTTACGAATGCCGTTGCCGCCCTGGGCATCTACGCTGTTAACATACGGAGCAACGTCTTTTGCCCCTGCCAGTGCGGCTAATTACCCTTGTGCGTCGCTATATACATCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTTCCAGTCCTTCGGAAGTCACCGTGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTGCTTccacca                              <  1:71032/245‑1 (MQ=255)
                                                                                                                                                                                                                      gCAGCGCCTGAAACGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGCGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATccc                          <  2:395266/250‑1 (MQ=255)
                                                                                                                                                                                                                        aGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGGCGCGCAGGCGCTTCGCTTCTTCCACCATccctgc                       <  2:603724/251‑4 (MQ=255)
                                                                                                                                                                                                                          cTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGACTAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTgcgc                    >  2:110660/1‑251 (MQ=255)
                                                                                                                                                                                                                             cTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTc                 >  2:236284/1‑251 (MQ=255)
                                                                                                                                                                                                                               gAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTcgc               >  1:524273/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTc          >  2:97896/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTACTATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCGTCGGAAGTCTCCGGTATTTTCTCCACAATACCCGGAATAGCGTCGCGCAGGCGATTCGCTGCTTCCACCATCCCTGCTCGTCGCGGCTCAt        <  2:588150/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGt     >  1:31651/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCt   <  1:208383/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCt   <  1:148783/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCTGCGCGTCGCGGCTCATGGTCTg  >  1:392962/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTG  >  NC_000913/4139047‑4139534

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