Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,1580GC22.2% 216.7 / 15.6 90G304G (GGC→GGGaraGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (17/53);  new base C (17/3);  total (34/56)
Fisher's exact test for biased strand distribution p-value = 1.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGTGGTGATCGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTA  >  NC_000913/1983920‑1984395
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
gaaTGTGGTGATCGGCTTTGTTTTTTTCCCACCCGTTGTTGATTACACAACCGCCGAGCCCATGTTTTCGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCg                                                                                                                                                                                                                                                    <  1:1675495/239‑1 (MQ=255)
           tcgtctttGTTTTCTTCCCACCCTTGGTTGATTCCACAGCCGCCGAGCACATGTTTACGTCTGGCACTGATGTGGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGcc                                                                                                                                                                                                                                       <  2:828418/233‑1 (MQ=255)
           tCGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGACAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGGTGGGCGTTCTAACCCCGACCCGCGGGGGGTTGCGTGCCCCCGAACacg                                                                                                                                                                                                                                          >  2:1680035/1‑235 (MQ=255)
            cgtcATTGTTTTCTTCCCACGCGGTTTGGTGTACAAAACCCGCGAGCACAGGTTTTCGTCTGGATCTGATGGTGATATTGGCGCGAACGGGGGGCACGGGAATAATGCCTTCCGCTTTGCGAGCTTCCGGGCAGTGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCt                                                                                                                                                                                                                                      <  2:2108907/236‑1 (MQ=255)
            cGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGGATTTCGTCCCGCCAAAAATGc                                                                                                                                                                                                                                        >  1:1339240/1‑237 (MQ=255)
             gactttGTTTTCTTCCCAACCGTTTTTGATTCCACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCt                                                                                                                                                                                                                                      <  2:2026862/234‑1 (MQ=255)
                    ttttCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGGTTTGCGTCCCGCCAAACATg                                                                                                                                                                                                                                          >  1:1492344/1‑228 (MQ=255)
                         ttcCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGGCGATATAAACCTGAACGGCGGTGGATTTGGGTCGCGCCAAACAGCCCttctatgatttc                                                                                                                                                                                                                           >  2:198407/1‑228 (MQ=255)
                         ttcCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGGTGGTCGGTATAAACCTGACCGGCGGTGGTTTTCGTCCCCCCAAACAAGCCTTTCTTTAATTct                                                                                                                                                                                                                          >  1:908020/1‑237 (MQ=255)
                          tcCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGGTTTGCGTCCCCCCCAACAATCCTTTCCTTAAtttcct                                                                                                                                                                                                                        >  1:2106386/1‑235 (MQ=255)
                              aaCCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTTCGTCCCGCCCAAACTGCCTTTCATTAAATCGCTAcg                                                                                                                                                                                                                     >  1:237485/1‑238 (MQ=255)
                              aaCCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTTCCTTCCCCCAAAAAAACCTTTCCTTAATTCCCTAcc                                                                                                                                                                                                                     >  1:1595374/1‑237 (MQ=255)
                               cccGTTTTTGATTACCCACCCGCCGAGCACATGTTTACGTCTGGCCCTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  <  1:1888597/238‑1 (MQ=255)
                                    tgtggATTCAACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgcc                                                                                                                                                                                                             <  2:237485/239‑1 (MQ=255)
                                        gATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTa                                                                                                                                                                                                         <  2:890474/239‑1 (MQ=255)
                                          ttACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGGTTTTCGGCCCCCCaaa                                                                                                                                                                                                                                              >  2:689746/1‑202 (MQ=255)
                                           tACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGGCCCCCCCAACAATCCTTTTATTAATCCGCCACcgccccccccaccca                                                                                                                                                                                                      >  1:711959/1‑234 (MQ=255)
                                              acaACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTGCGTCCCCCCACACATGCCTTTTTTTAATTCTCTACCCCCCGCCCCCCCCAGa                                                                                                                                                                                                   >  1:873088/1‑238 (MQ=255)
                                                        aGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGGTTTTGGTCCCGCCAAACATTCCTTTCATTAATTCGCTACGCCCCCCCCCCACCCGACCCAACAAccc                                                                                                                                                                                        >  1:159805/1‑239 (MQ=255)
                                                          cacaTGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGGCGGTATAAACCTGACCGGCGGGCGTTTTGCTCACGCCCCAAACTGCCTTCATTAATGCGTCACCGCCCTCCCCTACCCCGACCAAACACCCca                                                                                                                                                                                       >  2:765504/1‑239 (MQ=255)
                                                          caagtGTTTGCGTCGGGCACGTGTATTCATATTGGCGCAAACGCAATGCACGGGAATACTGCCTTCCGCTTTGCGACTTCCGGGGCAGAGCCTCATGCCTGGGGCAATGGCGTGGCTCGGTTTACTAATCTCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCCTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCca                                                                                                                                                                                       <  2:845963/234‑1 (MQ=255)
                                                                tttACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCaaa                                                                                                                                                                                              <  1:801984/225‑1 (MQ=255)
                                                                tttACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCaaa                                                                                                                                                                                              <  1:1796642/225‑1 (MQ=255)
                                                                  tACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTTCGTTCCCCCCAACAAGCCTTTTCTTTATTTGCCACGCCCCCCCCCCCCCAAACCAAACAACCCCACAAtttt                                                                                                                                                                               >  1:1340344/1‑237 (MQ=255)
                                                                   aCGTCTGGCACGGATGTTGATATTGTCGCTACCGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTcac                                                                                                                                                                             <  1:1895462/239‑1 (MQ=255)
                                                                       cTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGGCCCCCCAAAAATGCCTTTCATTAAATTGCTACGCCCCGCCCCCACCCAACCAAAAAACCCAACAATTTcaccaca                                                                                                                                                                         >  1:558641/1‑238 (MQ=255)
                                                                        tGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCGTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTg                                                                                                                                                                        <  1:765504/239‑1 (MQ=255)
                                                                          gCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCCTCCCCCCCAACAAGCCTTTCCTTAATTTCCCACCCCCCCCCCCCACCACACCCAACAACCCAACAAATTTACCACcgcg                                                                                                                                                                      >  1:2091137/1‑239 (MQ=255)
                                                                          gCAATGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCACAAATGCCCTTTCTTTAATCTCTACGCACCTCCCCTACCCGACACAACAGCCccaca                                                                                                                                                                                    >  2:801984/1‑225 (MQ=255)
                                                                                aTGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGGCGATATAAACCTGACCGGCGGGGTATTTCTGTCCCCCCAACAAGACCTTTCTTTAATCGCGCCGCCCCCCCCCCTCCCCACCCAAAAACCCCAACAATTTACCCACTGTGAAcc                                                                                                                                                                 >  2:333967/1‑238 (MQ=255)
                                                                                   ttGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGcc                                                                                                                                                              <  1:1761327/237‑1 (MQ=255)
                                                                                         ttGTCGCGAACGGAGTGGACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCt                                                                                                                                                                                                                                      <  1:1749276/162‑1 (MQ=255)
                                                                                           gTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTGCGTGCCGCCCAAAAAGGCCTTTATTAAATTGCCACGCCCCGCCCCTACCAGACCAAAAAGCCcaaaa                                                                                                                                                                                    >  1:223845/1‑206 (MQ=255)
                                                                                             cgcgAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTTCGTTCCGCCAAACAAGCCCTTTATTAATTCTCCACGCCCCCCCCCCCCCCAACCCAAAAAGCCCCACAATTTTCCCCCTGCAAACCCCCACACTTATTg                                                                                                                                                    >  2:1761327/1‑238 (MQ=255)
                                                                                             cgcgAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGGCGATATAAACCTGACCGGCGGGGGTTTTGGTCCCCCCACACATGCCTTTCATTAATATGCTCAGACCCGCCCCTACCCGCACACACCACCCAACCAATTTACCACTGCGCACCCCCCGACTTATTgta                                                                                                                                                  >  2:630778/1‑237 (MQ=255)
                                                                                              gcgAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCg                                                                                                                                                 <  1:198407/239‑1 (MQ=255)
                                                                                                gAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGGGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGACCCGCCCCTACCAGACCAAACAGCCca                                                                                                                                                                                       <  1:964399/200‑1 (MQ=255)
                                                                                                   cGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGCCGATATAAACCTGTCCGGCGGGGTGTTGTTGCCCCCCCAACATGCCctttacttaat                                                                                                                                                                                                                            >  2:1749276/1‑154 (MQ=255)
                                                                                                        tGCACGGGAATAATGCCTTCCGCTTTGCGACCTTCCGGGCAGAGCTTCCTGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGTTCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGcc                                                                                                                                                              <  2:786290/216‑1 (MQ=255)
                                                                                                         gacggggAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCGGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCt                                                                                                                                      <  1:93951/235‑1 (MQ=255)
                                                                                                         gCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCt                                                                                                                                      <  1:2063919/239‑1 (MQ=255)
                                                                                                               gAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGATTATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAAt                                                                                                                                                                                                                             <  1:1225604/149‑1 (MQ=255)
                                                                                                               gAATAATGCCTTCCGATTTGCGTGCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGc                                                                                                                                 <  2:774131/238‑1 (MQ=255)
                                                                                                                  taatGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGATTTGCGTCCCCCCCAAAATGCCTTTTCTTAATTCGCTACGCCCCGCCCCCACCCAACCAAAAAACCCAAAAAATTTCCCACTGCGAACCGCCCGACTTATTg                                                                                                                                                    >  1:786290/1‑216 (MQ=255)
                                                                                                                   aatGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTca                                                                                                                            <  1:214955/239‑1 (MQ=255)
                                                                                                                   aagtccTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCCCCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTca                                                                                                                            <  1:1680035/235‑1 (MQ=255)
                                                                                                                    atGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCt                                                                                                                                                                                                                                      <  1:1625234/135‑1 (MQ=255)
                                                                                                                    atGCCTTCCGCTGTGCGATCTTCCGGGCAGAGTATCATGCCTGCGGCATTGGCGTGGCTCGGTTTACGAAAATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTcac                                                                                                                           <  2:908020/239‑1 (MQ=255)
                                                                                                                    aatgcTTCCGCTTTTGCGACTTGCCGGGAGAAGCATCGTGCCTGCGGAAATGGGTTGGCTCGGTTTACGAATTTCGCTCGCCTGTTGGTGGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCCCCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTca                                                                                                                            <  2:1147177/235‑1 (MQ=255)
                                                                                                                     tGCCTCCCGTTTGCCTATCTTCCGGGCAGAGCATCCTGCCTGCGGCAAGTGCGTGGCTAGGTTTGCGAATATCGATCGGCTTTTGGTCGATATAAACCTGCCCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTcaca                                                                                                                          <  2:1339240/239‑1 (MQ=255)
                                                                                                                        cTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCAGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGc                                                                                                                        <  1:1441434/238‑1 (MQ=255)
                                                                                                                         ttCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  >  1:670302/1‑149 (MQ=255)
                                                                                                                         ttCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  <  2:670302/149‑1 (MQ=255)
                                                                                                                         ttCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCAGCGGGGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCCCCCCTCCCCAACCAAACAACCCCAAAATTTTCCCACTCCGAACCGCCCGCCTTATTGGCTGACGCACGCCTGGTGCTTTCACACCa                                                                                                                       >  2:1441434/1‑238 (MQ=255)
                                                                                                                                       tCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGgc                                                                                                          <  1:1594947/238‑1 (MQ=255)
                                                                                                                                         ttCCGGGCAGAGCATCATGGCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaac                                                                                                                                                                                     <  1:2063949/162‑1 (MQ=255)
                                                                                                                                          tCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAACCCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGctc                                                                                                      <  1:918182/239‑1 (MQ=255)
                                                                                                                                            cGGGCAGCGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGCCCAAACAGCCCAACAATTTCCCCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGctcct                                                                                                    <  1:1522102/239‑1 (MQ=255)
                                                                                                                                            cGGGCAGAGCATCATGCTTGCGGCATTGGCGGGGCTCGGTTTACGAATATCGACCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGctcct                                                                                                    <  1:2063423/239‑1 (MQ=255)
                                                                                                                                            cGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGctcct                                                                                                    <  1:1022103/239‑1 (MQ=255)
                                                                                                                                                 cAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaac                                                                                                                                                                                     <  1:1745503/155‑1 (MQ=255)
                                                                                                                                                 cAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTcac                                                                                                                                                                             <  1:1966018/163‑1 (MQ=255)
                                                                                                                                                    agCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAAc                                                                                                                                                                                                                                             <  1:656409/99‑1 (MQ=255)
                                                                                                                                                       atcatGCCTGCGGCAAGGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAgg                                                                                                           <  1:39120/222‑1 (MQ=255)
                                                                                                                                                              cTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGGTTTTCGTCCCCCCCAACAAGCCcttta                                                                                                                                                                                                                                  >  2:656409/1‑97 (MQ=255)
                                                                                                                                                                  ggCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGAAATAAACCTGACCGGCGGTTATTTGCGTCCCGCAAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGACCCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGc                                                                                <  1:1366536/238‑1 (MQ=255)
                                                                                                                                                                         gCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGcc                                                                                                                                                              <  1:1161428/153‑1 (MQ=255)
                                                                                                                                                                                    tttACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCaaa                                                                                                                                                                                              <  1:644994/111‑1 (MQ=255)
                                                                                                                                                                                    tttACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGcc                                                                                                                                                              <  1:1682511/142‑1 (MQ=255)
                                                                                                                                                                                    tttACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGc                                                                                                                                                  <  1:1582929/154‑1 (MQ=255)
                                                                                                                                                                                           atatCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGc                                                                                                                        <  1:564721/173‑1 (MQ=255)
                                                                                                                                                                                              tcgatcgGCTGTTGGTCGATATAAACCTGACCGGCGGGGGTTTGGGTCCCGCCAAAAATGCCCTTCTTTAATTCGCTACGCCCCGCCCCCACCAGACCCAACAGGCccaca                                                                                                                                                                                    >  2:644994/1‑111 (MQ=255)
                                                                                                                                                                                              tcgatcgGCTGTTGGTCGATATAAACCTGACCGGCGGGGGTTTGCGTTCCCCCCAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCCACCAGACCCAAAAACCCCACAAATTTACCACTGCGCAACCCCCAAACTtttt                                                                                                                                                     >  2:1682511/1‑142 (MQ=255)
                                                                                                                                                                                                    gcgtGTTGCTCCATATACCCCTGCCCGTCGGTGATTTGCTCCCCCACAAACATGCCTTTCTTTAGTCCGCACCGCCCCGCCCCGCCCAGACCAAACGGCCCAACAACTTACCCACTGCGAACCGCCAGACTTATTGGCGTCGGCCCGCCTGGTGCTTTCACAGCATCAAGCCGTAGGCCCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCa                                              <  2:1662801/235‑1 (MQ=255)
                                                                                                                                                                                                     gCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  <  2:1363618/75‑1 (MQ=255)
                                                                                                                                                                                                     gCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  >  1:1363618/1‑75 (MQ=255)
                                                                                                                                                                                                      cttttgtccGATATAACCCTGCCGGGCGTTGATTTGCGCCCCCCAAACCATGCTTTCCATTAATTCGCTCCGCCCCGCCCCTCCCAGACCAAACACCCCAACAATTTCACCCCTGCGACCCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATc                                            <  2:1340344/230‑1 (MQ=255)
                                                                                                                                                                                                                 tataACCCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCCCTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGc                                                                                <  1:1587011/191‑1 (MQ=255)
                                                                                                                                                                                                                 tataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCt                                                                                                                                                                                                                        <  1:1031428/58‑1 (MQ=255)
                                                                                                                                                                                                                 tataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCaaa                                                                                                                                                                                              <  1:939708/82‑1 (MQ=255)
                                                                                                                                                                                                                 tataAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAAc                                                                                                                                                                  <  1:1921945/109‑1 (MQ=255)
                                                                                                                                                                                                                     aaCCTGCCCGCCGGGTATTTGCGTCCCGCAAAACATCCCTTTAATTAATTCGCCACGCCCCGCCCCGACCAGACCAAACAGCCCAACAAGTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGcgt                            <  2:2106386/239‑1 (MQ=255)
                                                                                                                                                                                                                           aCCGGCGGTGGTTTTCGGCCCGCCAAAAATGCCTTTCATTAAATCGCTACGCCCCGCCCCTACCCGACCAAACAGCCCAACAATTTCACCACCGCGAACCGCCAGACTTATTGGCGTACGCCCGCCTGGGGCCTTCACAGAATCCAGACGGAGGCGCCCCCCCCCATCACTAACCGGGTTCCACCCGTAGa                                                                      >  2:1587011/1‑191 (MQ=255)
                                                                                                                                                                                                                           aCCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  >  1:91880/1‑53 (MQ=255)
                                                                                                                                                                                                                           aCCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                                                                                                                                                                                  <  2:91880/53‑1 (MQ=255)
                                                                                                                                                                                                                           aCCGGCGGTGATTTGCGTCCCCCCCAACAATCCTTTCATTAATTCGCTACGCCCCCCCCCTCCCAGACCCCACAACCcaaca                                                                                                                                                                                    >  2:939708/1‑82 (MQ=255)
                                                                                                                                                                                                                             cgggggTGTTTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCa                    <  1:333967/239‑1 (MQ=255)
                                                                                                                                                                                                                                    tGATTGGCGCCCCCCCAACCATGCCTTTCATTAATTCGCTCCGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACctg               <  1:630778/238‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAgg                                                                                                           >  2:1363475/1‑134 (MQ=255)
                                                                                                                                                                                                                                                  cGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAgg                                                                                                           <  1:1363475/134‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTa  >  1:682049/1‑239 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
GAATGTGGTGATCGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCAGAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTACGCACGCCTGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGGTTGCCAGCCGTAGATATCGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTA  >  NC_000913/1983920‑1984395

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