Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,185,1990TG29.6% 134.8 / 21.8 73T35P (ACC→CCC) thiMhydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (15/35);  new base G (21/0);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATAAAAGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCG  >  NC_000913/2184949‑2185448
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
gATAAAAGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGt                                                                                                                                                                                                                                                                      <  2:295395/251‑1 (MQ=255)
              aCAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCAcatcat                                                                                                                                                                                                                                                       <  2:543780/251‑1 (MQ=255)
                   aaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAAGGTTTCACCACATCATTGGTc                                                                                                                                                                                                                                                 >  1:185101/1‑251 (MQ=255)
                      ggcggcGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                                                                                                                                                                                                                             >  1:443086/1‑251 (MQ=255)
                               aaTCGAGCGCACCCGCCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGcaca                                                                                                                                                                                                                                    <  2:443086/251‑1 (MQ=255)
                                          cccaccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCCCCACATCATTGGTCATGCAGTGCACAAGAGGGGAATg                                                                                                                                                                                                                         <  2:80758/251‑1 (MQ=255)
                                          cccaccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAAGGTTGGACCACACCATTGGTCAAGCAGAGCCCAAGAGGGGGATg                                                                                                                                                                                                                         >  2:134434/1‑251 (MQ=255)
                                           ccaccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCAATGGTCAAGCAGTGGACAAGAGGGGGATGt                                                                                                                                                                                                                       >  1:137305/1‑251 (MQ=255)
                                             accGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGTACCACATCAGTGGTCAGGGAATGCACAagagagggaggttg                                                                                                                                                                                                                      >  1:121251/1‑241 (MQ=255)
                                               cGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGt                                                                                                                                                                                                                   <  2:238685/251‑1 (MQ=255)
                                                      ggATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGaaaa                                                                                                                                                                                                             <  1:108839/249‑1 (MQ=255)
                                                        ntCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGGTCTGCAGTGCACAAGAGGGGAATGTGGGGGAAAAAggg                                                                                                                                                                                                         >  2:241641/2‑250 (MQ=255)
                                                              cGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGGGCACAAGAGGGGAATGGTGGTGAGAAAGGGGTaa                                                                                                                                                                                                    >  2:171208/1‑249 (MQ=255)
                                                                 ccAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGGAAAAAGGTGTAACGCGt                                                                                                                                                                                               <  1:206937/251‑1 (MQ=255)
                                                                   aGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgc                                                                                                                                                                                             <  1:510828/251‑1 (MQ=255)
                                                                            gagaGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgc                                                                                                                                                                                    <  2:285153/251‑1 (MQ=255)
                                                                                     ttgcttgcTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACCTCCATGGTCATGCAGGGCCCAAAGGGGGAATGGTTGGGGAAAAAGGGTAACGCGTGCGCAGCTTTCGCTggatcc                                                                                                                                                                           >  2:39968/1‑246 (MQ=255)
                                                                                           gctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTTCACCCACTCATTTGTCCATCAGGGCCCAAAAGGGGGATGTTGGGGAAAAAAGGGTAAACCGGGCGCAGATTGCGCTGAAACCAGCAgg                                                                                                                                                                   >  2:222800/1‑251 (MQ=255)
                                                                                               aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCACATCCTTTGGCATGCAGGGCACAAGAGGGGGATGTTGGGGGAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAACCCGCAgggcgg                                                                                                                                                               >  1:63894/1‑247 (MQ=255)
                                                                                                  ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACtt                                                                                                                                                           <  1:273464/251‑1 (MQ=255)
                                                                                                      gcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGGACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                                                                                                                                       >  1:224764/1‑251 (MQ=255)
                                                                                                      gcaCGCATCGCCTGAGCGCGTGGCTGCGGCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCATCAGGGCGTCTTGCAt                                                                                                                                                       <  2:495282/251‑1 (MQ=255)
                                                                                                           cATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGttt                                                                                                                                                  <  1:23211/251‑1 (MQ=255)
                                                                                                            aTCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGATTTCACAACACCAATTGTCATGCCGGGCCCAAGAGGGGGAAGGTGGTGAAAAAAGTGTAACGCGTGCGCCGATTGCGCTGGAccccgccgcgcgccctgcatattttt                                                                                                                                                   >  2:211295/1‑228 (MQ=255)
                                                                                                              cGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAGCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCt                                                                                                                                               <  2:527786/251‑1 (MQ=255)
                                                                                                               gCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTc                                                                                                                                              <  1:376877/251‑1 (MQ=255)
                                                                                                                 cTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGGCGTGCACAAAAGGGGGATGTTGGTGAAAAAAGGGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTTCTCCt                                                                                                                                            >  1:502778/1‑251 (MQ=255)
                                                                                                                  tGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTTAAGGGTTGGAACACACCCTTGGTCATGGAGGGCACAAGGGGGGGGTGTTGGGGGAAAAAGGGTCAACCGTGCGCAGATTGCGCGGAACCCAGCAGGGAGACTTGctcagatttcttctg                                                                                                                                           >  2:528881/1‑237 (MQ=255)
                                                                                                                      cgcgTGGCTGCGTCAGTTTGCAAACGTTAATCAACAAAGCACTGTCGATAGCCGCAAACTGTCTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCTGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCATATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGATCCTGCCaa                                                                                                                                       <  2:352032/251‑2 (MQ=255)
                                                                                                                                cGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCCTTTGTCCTGCAGTGCACAAAAGGGGGATGTTGGTGGAAAAAGGGTAACGCGTGCGCAGATTTCGCTGGAACCCGCCGGGCGACTTGCATAGTTTGCTCCTGCCCTAAAGTgg                                                                                                                                >  2:395418/1‑247 (MQ=255)
                                                                                                                                   cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGaagaag                                                                                                                           <  1:134434/250‑1 (MQ=255)
                                                                                                                                   aggtgtGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGc                                                                                                                          <  2:16747/249‑1 (MQ=255)
                                                                                                                                         gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATACCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAc                                                                                                                   <  2:445371/251‑1 (MQ=255)
                                                                                                                                         gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAc                                                                                                                   <  1:172900/251‑1 (MQ=255)
                                                                                                                                           cAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCCCATCATTGGGCATGCAGTGCACCAGAGGGGGGTGGTGGTGGAAAAGGGGTAACGCGAGCGCAGATTTCGCTGAACCCCGCCAGGCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAGAAGCaaagaact                                                                                                                  >  1:490011/1‑245 (MQ=255)
                                                                                                                                            aaCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACAACAATGGTCATGCCGTGCACAAGAGGGGAATGTTGGTGGAAAAGGGGTAACGCGTGCGCCAATTTCGCTGGACCCAGCCGGTCGACTTGGATAGTTTTCTCCTGCCATAAAGTGCAGAAGCAATGACCtg                                                                                                                >  1:290500/1‑251 (MQ=255)
                                                                                                                                             aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTTAAGGGTTGGACCCCCCCATTTGTCCTGCAGGGCACAAGAGGGGGAGGTTGGTGGAAAAGGGGTAACGCGTGCGGAGATTTGGCCGAACCCAGCAGGGCGACCTGCCTAGTTTGCTTCTGCCCTAACGGGaagaagaaaaggacttg                                                                                                                >  2:431142/1‑240 (MQ=255)
                                                                                                                                                     caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGt                                                                                                       <  1:524745/251‑1 (MQ=255)
                                                                                                                                                             cACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTTGTCATGCAGTGGACAAGAGGGGGATGTTGGGGAAAAAGGTGTAACGCGGGCGCCGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTccc                                                                                               >  1:144937/1‑251 (MQ=255)
                                                                                                                                                              aCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGACTTcccc                                                                                              >  1:343350/1‑250 (MQ=255)
                                                                                                                                                              aCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGTTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGTGGAAAAAGGGTAACGCGTGCGCAGGTTTCGCTGAACCCCGCAGGTCGACCTGCATAGTTTGCTCCTGCCATAACGGGAAGAAGCAAAGACCTGGTGGTCCGTGACTTcccc                                                                                              >  2:135558/1‑250 (MQ=255)
                                                                                                                                                                gggCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTAc                                                                                            <  1:467/250‑1 (MQ=255)
                                                                                                                                                                   cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCt                                                                                         <  1:39968/251‑1 (MQ=255)
                                                                                                                                                                     aTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                                                                                                                                 <  2:290344/179‑1 (MQ=255)
                                                                                                                                                                     aTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAgggc                                                                                                                                                                 >  1:290344/1‑177 (MQ=255)
                                                                                                                                                                            aaaaCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATc                                                                                <  2:284234/250‑1 (MQ=255)
                                                                                                                                                                               actgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAg                                                                             <  1:96076/251‑1 (MQ=255)
                                                                                                                                                                                  gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAATGGTTTCCCCCCCACATTGGTCATGAAGGGCAAAAGAGGGGAAAGTTGGTGGAAAAGGTGGAAAGCCGGCGCAGATTGCGCTGCAACCACCAGGGCGACGTGCCTAGTTTTGTCCTGCCAAAAAGTGAAAAACAAAACACCTGGGGGTCCCTTGCATCCCTACGCTGGCATCTTCCAGAtt                                                                          >  2:549834/1‑250 (MQ=255)
                                                                                                                                                                                   actgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCa                                                                        <  1:431142/251‑1 (MQ=255)
                                                                                                                                                                                       gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTg                                                                    >  1:253691/1‑251 (MQ=255)
                                                                                                                                                                                             tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATAcn                                                               <  2:253691/250‑2 (MQ=255)
                                                                                                                                                                                               ggTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCAACACATCACTGGCCAAGCACCGCACAAGAGGGGAATGATGGTGGAAAAGGGGGAACTCGTGCGCCGGCTTCGCTGCAACCAGCAGGGCGGCCTTGATAGTTTGCTCCTGCCATAATGTGAAGAAGCAATGACCTGGCGGTCCGCGAATTCCCTACGCTGGCAATATCCCGGTCATGTGAAACgggg                                                            >  2:100169/1‑250 (MQ=255)
                                                                                                                                                                                                   tCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATAATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTAtt                                                         <  2:486964/250‑1 (MQ=255)
                                                                                                                                                                                                     gATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAAACCCGCAGGTCGGCTTGCCTAGTTTGCTCCTGCCATAACGGGAAGAAGCAATGAACTGGTGGGCCGGGGCTTCCCTACGCTGGCATTATTCAGATCCAGGGGTACGGGGATTtct                                                      >  2:345891/1‑251 (MQ=255)
                                                                                                                                                                                                        aaCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCTCATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAg                                                   <  1:282780/251‑1 (MQ=255)
                                                                                                                                                                                                          ccATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACCTCATTGGGCATGCAGGGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGAATGCGCTGAAACCCGCAGGGCGACTTGCCTAATTTTCTCCTGCCCTAACGTGAAAAAGCAAATACCTGGGGGTCCGTGACTTCCCTACCCCGGCATTTTCCAGATCAGGTGATACGGGGAtttttccgcc                                                 >  2:508240/1‑244 (MQ=255)
                                                                                                                                                                                                              cGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCCTGCCGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGGAAAGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGCTTTCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTAATCAGATCAGGTGATACGGGTATTTCTCAGccctcc                                             >  1:191641/1‑247 (MQ=255)
                                                                                                                                                                                                                     gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCCACACATCATTGGTCATGCAGTGCACAAGAGGGGAATTTTGGTGGAAAAAGTGTAACGCGTGGGCAGAGTTAGCTGAACCCCACAGGGCGGCGTGGACAGTTTTCTCCTGCCCTAAAGGGAAAAAGCAATGAACCAGGGGTCCGTGAATTCCCTACCCTGGCATTATCCAGATCAGGGGGTACGGGGAATTCTCAGCCTTCCCGCAGa                                       >  2:156657/1‑250 (MQ=255)
                                                                                                                                                                                                                     gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGGAAAAAGGTGTAACGCGCGCGCAGAATGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGACTTCCCTACGCTGGCATTATCCAAATCAGGGGATACGGGGATTTCTCAGCCTTCCCGCAGa                                       >  1:247917/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAg                                     <  1:222800/251‑1 (MQ=255)
                                                                                                                                                                                                                       tGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCCCCACATCATTGGGCATGCAGCGCACAAGAGGGGGATGTTGGGGGAAAAGGTGTAACGCGTGGGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGGGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGGTACGGGGAATTCTCAGCCTTCACGCACAAgg                                    >  2:318166/1‑251 (MQ=255)
                                                                                                                                                                                                                         cACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGc                                  <  2:224764/251‑1 (MQ=255)
                                                                                                                                                                                                                          aCCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAAGTCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAAAAGCAATGACCTGGGGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGGGATACGGGGATTTCTCAGCCTTCACGCCGAAgggg                                  >  1:481580/1‑249 (MQ=255)
                                                                                                                                                                                                                           ccGAGCGCCAGCAAGGTATTGGCGGTAAAGGATTGCACCCCATCCTTGGTCCTGCAGTGCACAAAAGGGGAATGTTGGTGAAAAAGGAGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGGTTGCTCCTGCCATAAAGTGAAGAAGCAATGACCTGGGGGGCCGTGACTTCCCCACGCTGGCAATATCCAGATCAGGGGATACGGGGATTTCTCAGCCTTCACGCAGAAggggac                                >  1:33977/1‑248 (MQ=255)
                                                                                                                                                                                                                            cGAGCGCCAGCAAGGTATTGGCGGCAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAACAAAGACCTGGGGGGCCGTGACTTTCCTACGCTGGCCTTATCCAGATCAGGGGATACGGGGATTTTTCAGCCTTCACGCAGAAGGGCAcc                               >  1:358350/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gcCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAg                          <  2:185101/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGTATTGGCGGTAAAGGTTTTCACCAAATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAAGGGTAACGCGTGCGCCGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGGTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGGGGTAAGGGGATTTTTCCGCCTTCCAGCCGAAGGGCACCCCGAGGCGTTTg                   >  2:25860/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggTATTGGCGGTAAAGGGTTGCACCAAATCATTTGTCATGCAGGGGCCACGAGGGGGATGTTGGGGGAAAAGGTGTAACGCGGGCGAAGAATGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTTCTCCTGCCATAACGTGAAGAAGCAAAGACCTGGGGGTCCGTGACTTCCCTACGCTGGCCTTAACCCGAACAGGGGATACGGGGATTTCTCCGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTgg                  >  2:449317/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGtt                <  2:481580/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTg               <  2:358350/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATaa      >  1:108907/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        aGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGTTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAAt    <  1:463342/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATc   >  1:60088/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           gTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCg  <  2:144937/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
GATAAAAGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCG  >  NC_000913/2184949‑2185448

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