Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,390 G→T 42.1% R123L (CGC→CTC)  nemR → transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,3900GT42.1% 55.6 / 66.8 66R123L (CGC→CTC) nemRtranscriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base G (22/16);  new base T (12/16);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTCATCTGAAA  >  NC_000913/1726104‑1726671
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
gggATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCa                                                                                                                                                                                                                                                                                         >  1:646901/1‑289 (MQ=255)
    tGGGGCTAAGCGAATTACTAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGcgcg                                                                                                                                                                                                                                                                                    <  2:551689/289‑1 (MQ=255)
      gggCTAAGCGAATTACTAAAAACCGCGTAAGTGCCGAAATGGTCCTTATATCACTACTTTCGCTCTAAAGTAGCGTTTGGCGTTGCAATGCTTGAGCGTAATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCGTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGcgcg                                                                                                                                                                                                                                                                                    <  2:506153/288‑1 (MQ=255)
            aGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAAAACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCCGAAGATATGCTCAGCGCTATGGATa                                                                                                                                                                                                                                                                             >  1:506153/1‑289 (MQ=255)
                 aTTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGt                                                                                                                                                                                                                                                                        >  2:203270/1‑289 (MQ=255)
                               cTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGAGGATAAAGGTGCTCGCGGCGTGATc                                                                                                                                                                                                                                                         >  2:318483/1‑289 (MQ=255)
                               cTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGa                                                                                                                                                                                                                                                           >  2:364473/1‑288 (MQ=255)
                               cTGAAGTGCCGAAAGGGTCCTTCGATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGa                                                                                                                                                                                                                                                           >  1:186690/1‑288 (MQ=255)
                                tGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCGGCAATCCGGCAAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGACGTGTGCGCTCTGTCAGAAGATATGCGCAGCAGGATGGATAAAGGTGCTGGGGGCGTAATc                                                                                                                                                                                                                                                         >  2:604046/1‑289 (MQ=255)
                                      gCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTg                                                                                                                                                                                                                                                   <  1:74203/289‑1 (MQ=255)
                                                       aTCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGa                                                                                                                                                                                                                                  >  2:360015/1‑289 (MQ=255)
                                                                ttCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGAGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGt                                                                                                                                                                                                                        <  1:318483/289‑1 (MQ=255)
                                                                ttCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCATATTACCAGCAAACACTGGACCAGTTTTGCCAACATGGAACCATCACTGGTTGCCTGCCCGTAAACCTCTCTGCCGCAGTGTGCGATCCGTCAGTACGTATGCGCAGCCCGATTGATAAAGGTGCTCGCGGCTTGACCGCCCTTCCCTCTCTCGCCTGGGAAAATGGcct                                                                                                                                                                                                                         >  2:402888/1‑288 (MQ=255)
                                                                  cGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTg                                                                                                                                                                                                                       >  1:135707/1‑289 (MQ=255)
                                                                       tAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAAc                                                                                                                                                                                                                  <  1:312594/289‑1 (MQ=255)
                                                                                 tttGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATgg                                                                                                                                                                                                                            <  2:773685/269‑1 (MQ=255)
                                                                                 tttGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATgg                                                                                                                                                                                                                            >  1:773685/1‑269 (MQ=255)
                                                                                   tGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATAACACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTACCAACATGGAACCATCAGAGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTACGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCACGCGGCGTGATCGCCCTGCTGTCTCCGGCGCTGGAAAATGGCCTGGAGACCCATTGTTTaa                                                                                                                                                                                                        >  2:302544/1‑287 (MQ=255)
                                                                                       gTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTGACCtttt                                                                                                                                                                                                  <  1:17153/289‑1 (MQ=255)
                                                                                         tGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCCACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGCAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGAGCTCGCGGCGTGACCGCCCTGCTCTCTCAGGCGCTGGAAACGGGCCGTGAGAACCATTGTTTAACCTTTtgt                                                                                                                                                                                                >  2:520664/1‑289 (MQ=255)
                                                                                                    aGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCCATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGATTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGACACTGGCCGTGAGAACCATTGTTTAACCGTTTGTGGCGAAcc                                                                                                                                                                                        >  2:604384/1‑286 (MQ=255)
                                                                                                                 cgcATATCACCACCTCATTAATGGATTGCTGCAATCCGTCGAAGGCACCTACCGCGACCGCATACTGGCTTATTACCAGCACACACCGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACATTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACaa                                                                                                                                                                        <  2:760336/289‑1 (MQ=255)
                                                                                                                   cATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGcac                                                                                                                                                                          >  1:624109/1‑285 (MQ=255)
                                                                                                                   cATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGcac                                                                                                                                                                          <  2:624109/285‑1 (MQ=255)
                                                                                                                   cATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGt                                                                                                                                                                      <  1:285900/289‑1 (MQ=255)
                                                                                                                   cATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGt                                                                                                                                                                      >  1:457762/1‑289 (MQ=255)
                                                                                                                    atatCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTCTTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGTTGCTCGCGGCGTGATCGCCCTGCTCCCTCAGGCGCTGTAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCCCaacgg                                                                                                                                                                     >  2:447143/1‑286 (MQ=255)
                                                                                                                    atatCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGt                                                                                                                                                                      <  2:635578/288‑1 (MQ=255)
                                                                                                                         aCCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTa                                                                                                                                                                <  1:325647/289‑1 (MQ=255)
                                                                                                                           cAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACg                                                                                                                                                              <  1:364473/289‑1 (MQ=255)
                                                                                                                             gCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAATCTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGcgcg                                                                                                                                                                                                                                                                                    <  2:95900/169‑1 (MQ=255)
                                                                                                                             gCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGcgcg                                                                                                                                                                                                                                                                                    >  1:95900/1‑169 (MQ=255)
                                                                                                                                  tgactgaGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGAGAACCGCTGCAACAGGCACAAGTGCTTTACGCACtttg                                                                                                                                                       >  2:662233/1‑289 (MQ=255)
                                                                                                                                     gttaGTTGCTGAAATCCGGCGAGGGATGCTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCAAACATGGAACCGTCAGTCGTTGCCTTACCGTAAAACCCTCTGCCGAAGTGTTCGATCTGTAAGAAGATATGCTCAGCGCGATGGATAAAGGTTCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCt                                                                                                                                                    <  2:414967/286‑1 (MQ=255)
                                                                                                                                        aGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTg                                                                                                                                                  >  2:292870/1‑288 (MQ=255)
                                                                                                                                                         gAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCa                                                                                                                                <  2:135707/289‑1 (MQ=255)
                                                                                                                                                         gAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCa                                                                                                                                <  2:186690/289‑1 (MQ=255)
                                                                                                                                                                       cgACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAAGCTGCAGGCCAAAATTTCGCGCAGt                                                                                                                  >  1:536944/1‑289 (MQ=255)
                                                                                                                                                                         aCCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGtt                                                                                                                 >  1:468947/1‑288 (MQ=255)
                                                                                                                                                                                        attaCCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACaa                                                                                                                                                                        >  2:513285/1‑218 (MQ=255)
                                                                                                                                                                                        attaCCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACaa                                                                                                                                                                        <  1:513285/218‑1 (MQ=255)
                                                                                                                                                                                              aGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTg                                                                                                                                       >  2:66090/1‑245 (MQ=255)
                                                                                                                                                                                              aGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTg                                                                                                                                       <  1:66090/245‑1 (MQ=255)
                                                                                                                                                                                               gCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTg                                                                                          >  1:122858/1‑289 (MQ=255)
                                                                                                                                                                                                 aaaCACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCGCTGCCGCAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGAGTTATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGc                                                                                        <  1:41547/289‑1 (MQ=255)
                                                                                                                                                                                                    cacTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCa                                                                                     >  1:125390/1‑289 (MQ=255)
                                                                                                                                                                                                        gAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTa                                                                                 <  1:203270/289‑1 (MQ=255)
                                                                                                                                                                                                          aCCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCCTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaa                                                                               <  1:604384/289‑1 (MQ=255)
                                                                                                                                                                                                                tttGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaaaa                                                                            <  1:360015/286‑1 (MQ=255)
                                                                                                                                                                                                                 ttGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACatt                                                                        >  1:815914/1‑289 (MQ=255)
                                                                                                                                                                                                                              acgatcgttGGTAGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGGGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTg                                                           <  2:125390/281‑1 (MQ=255)
                                                                                                                                                                                                                                  tCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTAATGCGACGCCTGCCGt                                                       >  1:782316/1‑289 (MQ=255)
                                                                                                                                                                                                                                   cAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaa                                                                              <  1:852652/265‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaa                                                                              >  2:852652/1‑265 (MQ=255)
                                                                                                                                                                                                                                       ggTTGCCTGCAAGTAAAACTCCCTGCCGAATTGTGCGATCTGTCAGAAGATATGCTCAGCTCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAg                                                   <  1:534139/288‑1 (MQ=255)
                                                                                                                                                                                                                                           gCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGc                                              >  2:854885/1‑289 (MQ=255)
                                                                                                                                                                                                                                            ccTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACAATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTGGCGCGAATCTGCAGGCCAAACTTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAACAATTATTGCGACGCCTGCCGTTTAGCAGGcct                                            >  2:279749/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                  gAAGTGTGCGATCTGTCAGAAGATATGCTGAGCGCGATGGATAAAGCTGCTGGCGGCGTGATCGCCCTGCTCTCTCAGGAGCTGGAAAATGGCCGTGAGAACACTTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGTTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCTGGCATTTTTTATCACCAGACGACCgg                       <  1:207419/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCggg                      <  1:259500/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGTGTGCGATCTGTCAGAAGATATGCTCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCggg                      <  2:709896/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tgCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGACTCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAACCGCGCTGGCCCATTTAAAAAACATTATTGCGACGCCTTCCGGTTAGCAGGCATTTTTTCTCACCAGACGACCGGGAGCCt                 >  2:33150/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                               tGTCAGAAGATATGGGCAGAGCGATGGATGAAGGTGCTCGCCGCCTGATCGCCCTGCTATCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGTCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCTAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTc          <  2:637273/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  cAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCCCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGCCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGCCTGCCCGGTGCCTCTCTGTCa         >  2:215029/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                      agaTATGCGCAGCGCGATGGATAAAGGTGATCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTGGAGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTCATCTGaaa  <  1:279749/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           tGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTa                                                                                 >  1:216136/1‑206 (MQ=255)
                                                                                                                                                                                                                                                                                           tGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTa                                                                                 <  2:216136/206‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
GGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTACCAGCAAACACTGAACCAGTTTTGCCAACATGGAACCATCAGTGGTTGCCTGACAGTAAAACTCTCTGCCGAAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTCATCTGAAA  >  NC_000913/1726104‑1726671

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: