Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,455:1 +G 100% intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G93.8% 204.5 / 1.5 64intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (32/26);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.

TCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCG  >  NC_000913/3560179‑3560700
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
tggccggcGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGGGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGTGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                                           <  1:82216/271‑2 (MQ=255)
gcggcacaGGTGAGAAGGATGGTAAAATGTTCTTTTACCATGAACTGGTTAGCAACGTTGAGATTTTTGGTGACAATGCGCAAATTGTTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGTCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                                           <  1:152069/274‑2 (MQ=255)
    cACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTTCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTTCCCCGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATCATGGcac                                                                                                                                                                                                                                                          >  4:299094/1‑273 (MQ=255)
    cACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTAGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGCGCCAGGGCGTGTTAACCGACCGGGAAGGCAGCGCGCGACCGCCCGGATGGcgc                                                                                                                                                                                                                                                          >  4:62654/1‑273 (MQ=255)
       cGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                                   >  1:146555/1‑278 (MQ=255)
               atgatgCGAAAATCTTCTTTTACCATCAACGGGTTAGCAACGTTGAGATTGTTGGTGAAAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGc                                                                                                                                                                                                                                          <  1:42654/279‑1 (MQ=255)
                   tgCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                                           >  4:204718/1‑274 (MQ=255)
                       aaaaTCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                                                      >  2:339637/1‑275 (MQ=255)
                            cttcttTTACCATCAACGGTTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                                                                                                                                  <  2:62654/274‑1 (MQ=255)
                                     ccATCAACGGGTTAGAAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                                                                                                                                        <  2:342384/275‑1 (MQ=255)
                                      cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTTCCACGGCGTGTTAGCCGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGACACGGTTTTGCCCCGCCAGCCCATCGAGGt                                                                                                                                                                                                                       >  4:151276/1‑275 (MQ=255)
                                       agtaaCGGGTTGGGAACGTTGGGATTGTTGGTGACAAGGCGCAAAGTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCAGTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                                      <  2:81831/272‑1 (MQ=255)
                                         cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTACCCGAACTGGACGGCCGCGCTGCCCCGCCTTGATGGCCACGGATCAGGTTTTGCGCCGCCAGCCCATTGGGGTcgc                                                                                                                                                                                                                    >  4:180784/1‑275 (MQ=255)
                                                ttAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                                                                                                                              >  2:129608/1‑274 (MQ=255)
                                                        gTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGTAAGGCAGCGCCGCACCCCCATGATGTCAGAGGTAGCGGTTTTGCTCCGCCCGCCCATTTAGGTCgcggcggcacggcgcgg                                                                                                                                                                                                     >  4:89916/1‑264 (MQ=255)
                                                              aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                                                                                                                                                                              <  2:263498/275‑1 (MQ=255)
                                                                  ttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                                                                                                                          >  2:537981/1‑275 (MQ=255)
                                                                   tgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTAGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCCGGTTTTGCTCCGCCAGCTCATTGAGGTCTCGGCGAATAGTCTACGGGCTGACGgagcac                                                                                                                                                                                      >  1:119089/1‑275 (MQ=255)
                                                                    ggTGACAATGCGCAAATTGCTGTGATTGATCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                                                                                                                                        <  2:304668/275‑1 (MQ=255)
                                                                    ggTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgac                                                                                                                                                                                         >  2:471193/1‑273 (MQ=255)
                                                                          aaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCCCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCCCTAc                                                                                                                                                                              >  1:15281/1‑279 (MQ=255)
                                                                                 aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                           <  2:151276/275‑1 (MQ=255)
                                                                                 aaaTTGCTGGGACTTGGCAGTGGGTGGGTTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGGGCCCTTGGGGATTTGGCCCGCCACTTGGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                           <  4:537981/275‑1 (MQ=255)
                                                                                  agttGATGTGAGTGAGCAGGGGGGGCGCAACCGCTTACAGCGTGGGTGCGAGATCGATAAACAGCGGCGGGCCATGGGGGATTGGGTACGCCACTTTGCGGGCGATTGGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                           <  4:129608/272‑1 (MQ=255)
                                                                                       cTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCTGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                      <  2:474056/274‑1 (MQ=255)
                                                                                       cTGGGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                      <  2:270491/274‑1 (MQ=255)
                                                                                           gATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                                                                                                                 <  2:180784/275‑1 (MQ=255)
                                                                                               gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGGGCCGAAATCGCTAAACCGTGTCGAGCCATTGGGGATTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                             <  2:286717/275‑1 (MQ=255)
                                                                                                  ccgggcgtgcgCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                          <  2:89916/271‑1 (MQ=255)
                                                                                                    gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                            >  2:119148/1‑239 (MQ=255)
                                                                                                    gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                         >  2:127582/1‑274 (MQ=255)
                                                                                                     tgcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGGTCAGGTTTTGCTCCGCCCGCCCCTTGCGGTCGCGGCGAATATTCTGCGGGCTGACGGCGAAATGCTCTACCCGCTCTTCGGTTCTGACATAACCctg                                                                                                                                                   >  1:60308/1‑279 (MQ=255)
                                                                                                      gcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCCCAGGAACAGGTTTTGCTCGACCACCACCGTGAGGTGTGGGCGAATAGCGCTGCGGCGGACCGGAAACTGCCCGACCAACTCTTGGGATCTTACATAAcc                                                                                                                                                      >  4:286717/1‑275 (MQ=255)
                                                                                                        cggggttcccGCTTCGGGCGTGGTGCCTGTATGGGTAACACGGGTCGAGCCTTGGGGGAGTGTCTCGCCCCCTTTGCGGGCGTTGCGCGATTTTCCTTCGGACGGGGTGGCCTTGCGATAGTGGTACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCAAGCTAATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCGCTCTTCGGTACTGACATAACCct                                                                                                                                                    <  4:481890/267‑1 (MQ=255)
                                                                                                         tggccTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                    <  2:299094/270‑1 (MQ=255)
                                                                                                         ggcgtTACCGCTTCCGGGGTGGGGCGGATATGGATAAAAAGGGTCGAGCCATTGGGGATTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                    <  4:339637/273‑1 (MQ=255)
                                                                                                           cgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATACCCctgct                                                                                                                                                 >  2:137718/1‑275 (MQ=255)
                                                                                                               aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCCCATTGAGGTCGGGGCGAATAGTCTGCTGGCTGATGGAGTAATGCTCTACCAGCTTTTCGGTGCCGACATAACCctgctttgt                                                                                                                                             >  2:502193/1‑271 (MQ=255)
                                                                                                                   cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGGTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                    >  1:18552/1‑264 (MQ=255)
                                                                                                                   cTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                                                                                                                          >  2:110573/1‑274 (MQ=255)
                                                                                                                      ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAAATGGCAGCACNCGCCCCGCCGCCAGGATGGCGCAGGAACAGGCTTTGCTCCGCCAGCTCATGGAGGTCGCGGGCATTAGTAGGCGNGATGACGACGAGATGCTCGACACGCGCCTCNNTGATGGAATGAGCCATCTGTTTAACCCGTTCNa                                                                                                                                  >  3:65997/1‑279 (MQ=255)
                                                                                                                             ggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAACAGGTTTTGCTCCGCCAGCGCATTGa                                                                                                                                                                                                                          >  2:80766/1‑185 (MQ=255)
                                                                                                                                         gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                                                                                                                           >  2:428293/1‑251 (MQ=255)
                                                                                                                                           atagcacgtcgtggCCTTTGTGGATTTGCTCCGCCAGTGGGCGGGTGATTCTCGCTTTTTCTTCGGTTTGGGTTGCCTTGCAGTCGTGCCAAGGCGTGTGAGCCGAAATGGAAGGGTGCGCCGCACGGCAATGATGGCGCAGAATCATGTTTTTCTCCGCCAGCTCATTGAGGTCGCGGAGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTAGGTAATTACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                                                                                                                  <  4:23284/261‑1 (MQ=255)
                                                                                                                                            aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                                 <  2:204718/274‑1 (MQ=255)
                                                                                                                                             aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                       <  4:428293/251‑1 (MQ=255)
                                                                                                                                                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGGTCTACCAACCCTTCGGTAATGAAATAACCCAGCTGTTTAACCAGTTAGATAAATCCGTTTGGACATTGtgttt                                                                                                         >  4:509352/1‑275 (MQ=255)
                                                                                                                                                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                    >  1:172431/1‑233 (MQ=255)
                                                                                                                                                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCACCAGATCAGGTTTTGCTACGCCAGCTACTTGAGTTCGGTGCGAAGAGTCTGCGGGCTGCCGGCGCAACGCTCTACCACCACTTCGTGACTGCCCTAAACCTCCTGTTTACCCAGTTCCATAATACCCTTGTACCATTCtgttt                                                                                                         >  4:173433/1‑275 (MQ=255)
                                                                                                                                                    tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCAACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttt                                                                                                        >  2:492705/1‑273 (MQ=255)
                                                                                                                                                             tGGGTATTTGCTCCGCCGCTTTGCGGGCGTTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCTATCCTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                             <  1:65997/276‑1 (MQ=255)
                                                                                                                                                                    ttccTCCGCCCCTTGGCGGGCTAGGGTCGCTTTTGCTTCGGTCTGGGTGACCTTTCGATAGGTCCACGGAGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCCCAGGTTCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                                       <  4:502193/272‑1 (MQ=255)
                                                                                                                                                                                  gcgtggCGATGCGCTCTTTGTCTCGCGTGAGGCGGGGCGTGCGGTCGTGCCCCGGGGTGTTAACCGAACTGGAAGCCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCATCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                         <  4:349271/271‑1 (MQ=255)
                                                                                                                                                                                       gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCAGCCGGAGAAATGCTCTACCAGCCATTCGGTACTGCCCTAACCATGCTGTTTAACCAGGTTGATAATACCGTTTCTACGTGTGTTTTGTTTCATTCTCACATCCCTGGAATTTttttcgtttgc                                                                    >  4:334213/1‑273 (MQ=255)
                                                                                                                                                                                                 cgttttCTTCGGGCGGGGTGGCCTGGCGATCGGGCCACGGGGTGTTAACCGACATGGAAGGCAGCGCCGCACCGCAAGGATGGGGCAGGATCAGGTTTGGCTCCGCCAGCTCATTTAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAg                                                          <  4:137718/273‑1 (MQ=255)
                                                                                                                                                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTTACATAACCCTGCTGTTTAACCAGTTCGATAATACAGTAGTTACGTTGTGTTTGTTTCATTTATAAACCCCTGGAATTATTTTCGTTTTCACGCATTGGGCGAATGcaaca                                               >  4:427661/1‑275 (MQ=255)
                                                                                                                                                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                               >  2:259650/1‑274 (MQ=255)
                                                                                                                                                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAAGTTGTGTTTCTTTCATTTATAAATCCCGCGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                               >  2:552215/1‑274 (MQ=255)
                                                                                                                                                                                                                 gggggCCTTGCGGTCGGGCCACGGGGTGTTAACCGAATTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTGTGCGGGCTGACGGTGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                         <  4:425227/272‑1 (MQ=255)
                                                                                                                                                                                                                 ggTGGCCTGTCGATCGTGTCACGTCTTGTCAACCGAACTGGAAGGCAGTGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGATCCGCCAGCTCATTGAGGTAGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGATCTACCAGCTCTTCGGTACTGACATAACCCGGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                         <  2:30420/275‑1 (MQ=255)
                                                                                                                                                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTGTTCGGTACTGACATAACCCTGATGGTTAACCAGCTCGATAATACCGTTGTGACGTTGTGGTTGTTTCATGTATAAATTCCTGAAAGTAGTTTCGTTTTCGCGAATTAAGCGAATCAAAAAAAGcataccat                                   >  4:501505/1‑268 (MQ=255)
                                                                                                                                                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                             >  2:425227/1‑265 (MQ=255)
                                                                                                                                                                                                                        ttGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTCTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTa                                  <  2:333281/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTAAAAATCCCTGGAATTATTTTCGTTTTCGACCATTGACCAATTCCAAAAAACCCACCCCTAAACCCACGCCCACACCGGCAAtctgt        >  4:338360/1‑274 (MQ=255)
                                                                                                                                                                                                                                                     aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAAATGCTCTACCAGCTCTTAGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTATGTTTCTTTTATAAATACCTGGAATTATTTTCGTTTTCGCGCATTGAGCGATTCAACAAAAGCCATAGCAAAACCCACGGACAACCCGACGATGTGTGCTCCg  >  1:151841/1‑278 (MQ=255)
                                                                                                                                                                                                                                                         gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTTTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGAATTTAGCGAATCAACAAAAGCCATAGCTAAACCCACGCCCAACCCGGGCATGTGGGCCCCg  >  4:580304/1‑274 (MQ=255)
                                                                                                                                                                                                                                                          aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                                          >  2:104400/1‑107 (MQ=255)
                                                                                                                                                                                                                                                              cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                      <  4:104400/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
TCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCG  >  NC_000913/3560179‑3560700

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