Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,069,743 C→T A33T (GCA→ACA)  rutG ← pyrimidine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,7430CT97.1% 89.1 / NA 35A33T (GCA→ACA) rutGpyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/9);  minor base G (0/1);  total (25/10)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTC  >  NC_000913/1069496‑1069980
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
ccGCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGt                                                                                                                                                                                                                                            >  1:371712/1‑251 (MQ=255)
  gCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAACTGGAAAGATTGGGATCCAGCCCCATCAGCATCGGCCTTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCt                                                                                                                                                                                                                                          >  2:87722/1‑251 (MQ=255)
      aGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGc                                                                                                                                                                                                                                      >  1:415251/1‑251 (MQ=255)
               cTGATATTCGGGTAAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACCAAAGCCGAGCTGGAGCCGAGATAAATGGGGCCGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCACATTACTGTCGCCTGAGCAAACGGCa                                                                                                                                                                                                                              >  2:176100/1‑250 (MQ=255)
                     ttCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTc                                                                                                                                                                                                                       >  1:414192/1‑251 (MQ=255)
                              aTGCCCTGACCGTTAAAGCCCGTCGCGGCGATAACAACCCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAACGAACAACAGTAACGTGCCGATCCACGCCATTCAAATGGAAAGAATGGGATACAGCCCCATCCGACTCGGCATTAACCCCGTCGCGCAAAACATCGCCACCGCGTGTTGAACGACCATTACTGTAGTCTGGACCAACGGCAATCGTTCATCCGgcgc                                                                                                                                                                                                              >  2:314206/1‑251 (MQ=255)
                                                           gATCACCACGCCGACAAACGCCGCGCTGGAGCCGAGATAACTGGAGACGCGCCCGCCGGTCATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACACCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGAtt                                                                                                                                                                                 >  2:61443/1‑251 (MQ=255)
                                                             tCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGACGCtttt                                                                                                                                                                               >  2:280465/1‑251 (MQ=255)
                                                              caccacGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAACAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGAtttt                                                                                                                                                                               >  1:315413/1‑250 (MQ=255)
                                                                 cacGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACt                                                                                                                                                                           >  2:235504/1‑251 (MQ=255)
                                                                                          ccGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGc                                                                                                                                                  >  1:25673/1‑251 (MQ=255)
                                                                                                    tGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGt                                                                                                                                        >  1:62248/1‑251 (MQ=255)
                                                                                                      gggccgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGCCATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCCATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGACAGTGAGGAAACCCGAACATTGCCATCAGCTGtct                                                                                                                                      >  2:133997/5‑251 (MQ=255)
                                                                                                             gcccgccGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTa                                                                                                                                 >  1:470177/1‑249 (MQ=255)
                                                                                                                          aaaGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGc                                                                                                                  >  1:326938/1‑251 (MQ=255)
                                                                                                                                      tAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATc                                                                                                        >  1:425531/1‑249 (MQ=255)
                                                                                                                                                 tCCCCGACATTAAAATGCACAGATTGGGATCCAGCCCCATCAGCATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGACCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGCAAACCGAACATTGCCATCAGCTGTCTCCTTACGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATacc                                                                                           >  2:198717/1‑249 (MQ=255)
                                                                                                                                                 cccccGACATAAAAATGAAAAGATGGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCCCCAAACATCGCCACCGGGTGTTGACCGCCCATTACTGGCGTCTGAGCAAACGGCAATCGTTCATCCGGGGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCAGCAGCTGTCTCCTTAATGAGGTTAACAAGCAGTGCGTATCAGCGCGTGATAAc                                                                                           <  2:470177/250‑1 (MQ=255)
                                                                                                                                                   cccGACATTAAAATGTAAAGTGTGGGATCCAGCCCCATCAGAATGGGCATTAACTCCGTCCAGCCAAACATCGCCACCGCGTGTAGACCGCCCATTACTGTCGTCGGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTATGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTg                                                                                         <  2:262650/251‑1 (MQ=255)
                                                                                                                                                      gACATTCAAATGGAAAGCTTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCACTCGTTCATCCGGCGCGACCACGCCGCTCTCTGAAGAGGTCGCTTTTAACTGCCAGCGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACCAGCAGGGCGCATCAGCGCGTGATATCTGCgc                                                                                      >  2:177517/1‑250 (MQ=255)
                                                                                                                                                                   aaaGATTGGCATCCACCCCCATCAGAATCGGCATGAACCCCGCCGCGCCAAACATCCCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCa                                                                         <  2:19980/251‑1 (MQ=255)
                                                                                                                                                                      gATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGcc                                                                      >  1:14891/1‑251 (MQ=255)
                                                                                                                                                                        ttGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCg                                                                    >  1:17794/1‑251 (MQ=255)
                                                                                                                                                                           ggATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAACCATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACCTTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGAGCGCATCAGCGCGTGGTAACTGCGATCGAACCACACCAGCCCGt                                                                   >  1:409621/1‑249 (MQ=255)
                                                                                                                                                                           ggATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAg                                                                 <  1:235504/251‑1 (MQ=255)
                                                                                                                                                                            gacccAGCCCCAACAGAATCGGCATTAACACCGTCGCGCCAAACACCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgg                                                                <  1:173798/248‑1 (MQ=255)
                                                                                                                                                                             aTCCAGCCCCATCAGAATCGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCacac                                                                           >  2:344912/1‑239 (MQ=255)
                                                                                                                                                                             aTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCacac                                                                           <  1:344912/239‑1 (MQ=255)
                                                                                                                                                                                         cAGAATCGGCTGTAACACCGCCGCGCAAAACATCGCCACCGCGTGTTGAACGCCCTTTCCTGTCGTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACg                                                   <  2:497767/251‑1 (MQ=255)
                                                                                                                                                                                         cAGAATCGGCATTAACACCGTCGCGCCAAACAGCGCCCCCGCGTGTTGAACGCCCATTCCTGTCGTCTGAGCAAACGGCAAGCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGAc                                                    <  2:414192/250‑1 (MQ=255)
                                                                                                                                                                                                ggCATTCACACCGTCGCGCCAAACATCGCCACCGAGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGGGTGGTGTGACGATGAATc                                            >  1:127507/1‑251 (MQ=255)
                                                                                                                                                                                                     aaaCACCGTCGCCCCAAACACCGCCCCCGCGTTTTGACCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCTCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                       <  2:415251/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gtgtTGACCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTCCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGcc         <  1:453946/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gtgtTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGcc         >  1:173415/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTc  >  1:473623/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
CCGCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTC  >  NC_000913/1069496‑1069980

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