Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,985,1010TG40.3% 113.9 / 14.1 63T275P (ACC→CCC) unknown16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base T (8/29);  new base G (25/0);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCG  >  1/1984852‑1985340
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
ttAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGgtg                                                                                                                                                                                                                                                      <  1:362149/251‑1 (MQ=255)
           gTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTTAATGCCgg                                                                                                                                                                                                                                            >  1:177612/1‑249 (MQ=255)
                              gCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGGTAAAGGCGTCCTGAAAAAGC‑GGA‑TT‑ggga                                                                                                                                                                                                                     >  1:620353/1‑249 (MQ=255)
                                    cTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATc                                                                                                                                                                                                               <  1:595143/251‑1 (MQ=255)
                                          cAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGCCATGAAAAAGC‑ggggtgggaa                                                                                                                                                                                                                      >  1:95336/1‑232 (MQ=255)
                                            gCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGGCCTGAAAAAGC‑GGG‑TT‑GGGAATGATCAGAACGaa                                                                                                                                                                                                       >  1:302061/1‑251 (MQ=255)
                                                      gTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGG‑CTAAGAGGCGGAAATGCCGGCATGCAAAGGG‑GGT‑TT‑TGAAATGAttagcgcggagttgccgtag                                                                                                                                                                                             >  1:153103/1‑232 (MQ=255)
                                                           cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATc                                                                                                                                                                                                          <  1:204442/233‑1 (MQ=255)
                                                           cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAAAGCCGTCATGAAAAAGG‑GGG‑TT‑GGGAATGATCAGATcctgtctcttttaaacaat                                                                                                                                                                                        >  1:419512/1‑233 (MQ=255)
                                                           cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGG‑GTAAGCGGGGTCAATGCCGTCATGAAAAGGG‑GGA‑TT‑GGAAAAGATCAGAACGAATTTGCCGTCAAA‑ATg                                                                                                                                                                                       >  1:499585/1‑251 (MQ=255)
                                                            acTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTTAAAGCCGTCATGAAAAAGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGTCGTCAat‑attg                                                                                                                                                                                      >  1:447515/1‑249 (MQ=255)
                                                             cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGaaa                                                                                                                                                                                                                    <  1:83623/221‑1 (MQ=255)
                                                             cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATggg                                                                                                                                                                                     <  1:469600/251‑1 (MQ=255)
                                                             cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGGCCTGAAAAAGC‑GGG‑TT‑GGGAAAGATCAGATCGAATTTTCCGTCAat‑attgg                                                                                                                                                                                     >  1:41337/1‑248 (MQ=255)
                                                             cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAAGC‑GGA‑TT‑GGGAATGATTAGATCGAATTTTCCGTCAat‑attgg                                                                                                                                                                                     >  1:205387/1‑248 (MQ=255)
                                                                 aCTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGGGGGGTCAATGCCGTCCTGAAAAAGC‑GGG‑GT‑GGAAAAGATCAGATTGAATTTGCCGTCAAT‑ATTGGaaaa                                                                                                                                                                                 >  1:403236/1‑251 (MQ=255)
                                                                 aCTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGGCGTCATGAAAAGGCGGGA‑TG‑GGAAA‑GATTAGATCGAATTTGCCGTCAATAATGGGaaa                                                                                                                                                                                  >  1:487768/1‑251 (MQ=255)
                                                                    cATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAacac                                                                                                                                                                              <  1:629992/251‑1 (MQ=255)
                                                                    cATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCCATTCCGTCATGAAAAAGC‑GGG‑TT‑GGGAATGATTAGATCGAATTTTCCGTCCAT‑ATGGGGAAaca                                                                                                                                                                               >  1:451378/1‑250 (MQ=255)
                                                                    cATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAAAGCCGTCATGAAAAGGC‑GGG‑TT‑GGAAATGATCAGATAGAATTTTCCGTCCAT‑ATGGGaaaaaaa                                                                                                                                                                              >  1:579333/1‑248 (MQ=255)
                                                                      tCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCCATGCCGTCCTGAAAAGGC‑GGA‑TT‑GGAAATTATCAGATCGAATTTTCCGTCAAA‑ATTGGAAAACACGt                                                                                                                                                                            >  1:562230/1‑251 (MQ=255)
                                                                                gAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAACTCCGTCATGAAAAAGC‑GGG‑TT‑GGAAATGATCAGATCGAATTTTCCGTCAAT‑ATTGGGAAACACGTCCC‑TGGCaa                                                                                                                                                                  >  1:344922/1‑250 (MQ=255)
                                                                                 aaTATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAacc                                                                                                                                                                <  1:359669/251‑1 (MQ=255)
                                                                                  atatACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGCCATGAAAAAGC‑GGA‑TT‑GGGAATGGTCAGATCGAATTTGCCGTCAAT‑ATTGGaaaaaaaaacacttgcaacca                                                                                                                                                                >  1:578189/1‑234 (MQ=255)
                                                                                      aCCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGaaa                                                                                                                                                                                                                                     <  1:335116/182‑1 (MQ=255)
                                                                                       ccTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAGGC‑GGG‑TT‑GGAAATGATCAGATCGAATTTGCCGTCCAT‑ATGGGGAAACACGTCAC‑TGGCAACCACCgcg                                                                                                                                                          >  1:176572/1‑249 (MQ=255)
                                                                                                     tttCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GGAAAGGGGGGCAATGCCGTCCTGAAAAGGC‑GGG‑GT‑GGGAATGAACAGATCGGATTTGCCGTCAAT‑ATGGGGAAAAACGTCAC‑TGGCAAACACCGTGCCCTGCAACCCAtt                                                                                                                                            >  1:550977/1‑251 (MQ=255)
                                                                                                          cGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAGGC‑GGG‑TT‑GGGAATGATCAGAACGAATTTGCCGTCAAT‑ATGGGAAAAAACGTCAC‑TGGCAACCACCCCGCCTTGCAACCGAttttt                                                                                                                                         >  1:149173/1‑248 (MQ=255)
                                                                                                               ggACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAAGC‑GGG‑TT‑GGGAATGATCAGATCGGATTTTCCGGCCAT‑ATGGGGAAAAACGTCAC‑CGGCAACCACCGCGCCCTGCCACCGATTTTTCTGTaaa                                                                                                                                  >  1:341956/1‑251 (MQ=255)
                                                                                                                gACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTaa                                                                                                                                   <  1:567952/249‑1 (MQ=255)
                                                                                                                  cGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGGCAATGACGGCATGAAAAGGC‑GGA‑TT‑GGGGATGATCAGATCGAAAGTGTCGTCAAA‑ATTGGAAAACAAGTCAC‑TGGCAAACACCGGGCCTTGCAACCGATTTTTCTGTAAAGTg                                                                                                                               >  1:604373/1‑251 (MQ=255)
                                                                                                                                 aaTGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGcc                                                                                                                <  1:153488/251‑1 (MQ=255)
                                                                                                                                  aTGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAaccac                                                                                                                                                              <  1:324452/204‑1 (MQ=255)
                                                                                                                                      tGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGc                                                                                                           <  1:643489/251‑1 (MQ=255)
                                                                                                                                        aTTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCa                                                                                                         <  1:254664/251‑1 (MQ=255)
                                                                                                                                             gCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATggg                                                                                                                                                                                     <  1:538800/171‑1 (MQ=255)
                                                                                                                                              ccAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGc                                                                                                     <  1:216724/249‑1 (MQ=255)
                                                                                                                                              ccAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCg                                                                                                    <  1:634544/250‑1 (MQ=255)
                                                                                                                                                  tCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAat                                                                                               <  1:6991/251‑1 (MQ=255)
                                                                                                                                                  tCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAat                                                                                               <  1:493430/251‑1 (MQ=255)
                                                                                                                                                     gCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAAGC‑GGG‑TT‑GGAAATGATCAGATCGAATTTTCCGTCAAT‑ATGGGAAAAAACGTCAC‑TGGCAACCACCCTGCCTTTCAACCGATTTTCCTTTAAAGTGCGCTTTTCGGGAGCCAACGCCAAGGCGTGAatat                                                                                             >  1:52530/1‑250 (MQ=255)
                                                                                                                                                               gAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAACGGGGTCAATGCCGTCCTGAAAAAGC‑GGG‑TT‑TGAAAAGATAAGATCGAATTTTTCGTCCAT‑ATTGGGAAACACGGCAC‑TGGGAAACACCCAGCCTTGCAAACGATTTTTCTGTAAAAGGGGCGTTTCGGGAGCAAAAGCCAAAGCGGGGATATCCGTCATggtt                                                                                  >  1:287077/1‑248 (MQ=255)
                                                                                                                                                                aaaCGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGAt                                                                                   <  1:160373/249‑1 (MQ=255)
                                                                                                                                                                          gCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCCTGAAAAAGC‑GGA‑TT‑GGGAATGATCAGATCGAATTTTCCGTCAAT‑ATGGGAAAAAACGTCAC‑CGGCAACAACCCTGCCTTTCAACCGATTTTTCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAAGCGTGGATATCCGTCATGATTAAAC‑ggggggt                                                                      >  1:43963/1‑247 (MQ=255)
                                                                                                                                                                              cGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCggg                                                                  <  1:198368/251‑1 (MQ=255)
                                                                                                                                                                              cGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCggg                                                                  <  1:637680/251‑1 (MQ=255)
                                                                                                                                                                               gATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTCGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGTCATGAAAAAGG‑GGG‑TT‑GGGAAAGATCAGATCGAAATTTCCGTCAAT‑ATGGGGAAAAACGTCAC‑TGGCAACCACCGTGCCTTTCAACCGATTTTTCTGTAAAGGGCGCTTTTCGGGAGCCAACGCCAAGGCGTGAATATCCGTCATGAtt                                                                                  >  1:294497/1‑235 (MQ=255)
                                                                                                                                                                                aTGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGcc                                                                                                                <  1:242559/204‑1 (MQ=255)
                                                                                                                                                                                aTGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGCCATGAAAAGGC‑GGA‑TT‑GGGAAAGATCAGATCGAATTTTCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCCTGCCTTGCAACCCATTTTTCTGTAAAGTGCGCTTTGCGGGAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCggggt                                                                >  1:76843/1‑249 (MQ=255)
                                                                                                                                                                                    gcAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGCCGGCATGAAAAAGC‑GGG‑TT‑GGAAAAGATCAGATTGAATTTGCCGTCAAA‑ATGGGAAAAAACGTCCC‑TTGCAAACACCGGGGCTTGCAACCGATTTTTCTTTAAAATGCGCTTTTCGGGAGGCAACGCCAAGGGGGGAATATTCGTCCTGATTAATT‑GGGGGGTCGGGGtgtgg                                                            >  1:588087/1‑250 (MQ=255)
                                                                                                                                                                                       gCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGTAA‑CGGTGTCAATGCCGTCATGAAAAGGC‑GGATTT‑GGAAA‑GATTAGATCGAATTTTCCGTCCAT‑AAGGGGAAAAACGTCAC‑TGGCAACCACCCTGCCTTGCCACCGATTTTTCTGTAAAATGCGCTTTTCGGGAGCCAACCCCAAAGCGGGAAAATCCGTCATGATTAATC‑GGGGGGTCGGGTtggg                                                             >  1:478732/1‑245 (MQ=255)
                                                                                                                                                                                         tCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTg                                                       <  1:541599/251‑1 (MQ=255)
                                                                                                                                                                                               cccGGTTTCAAATACCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAatat                                                  <  1:213320/250‑1 (MQ=255)
                                                                                                                                                                                                   gTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAAGC‑GGA‑TT‑GGAAAAGATCAGATCGAATTTTCCGGCAAT‑ATGGGAAAAAACGTCAC‑TGGCAACCCCCGGGCCCTTCAACCGATTTTTCTGTAAAACGCGCTTTGCGGAAGCCAACGCCAAAGCGGGAAAATCCGGCATGATTAATCGGGGGGTTC‑GGTTGGGGTGGTGAATATAAcat                                             >  1:44034/1‑249 (MQ=255)
                                                                                                                                                                                                    tttCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGc                                            <  1:258680/251‑1 (MQ=255)
                                                                                                                                                                                                             ccATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCAcc                                    <  1:644109/250‑1 (MQ=255)
                                                                                                                                                                                                                     ccTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGcgccg                           <  1:201159/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAATGGCGTCATGAAAAAGC‑GGG‑TT‑GGAAAAGATTAAATCGAATTTGTCGGCAAT‑ATGGGGAAAAACGGCAA‑TGGCAAACCCCCTGCCCTGCAACCGATTTTTCTGTAAAATGCGCTTTTCGGGAGGCAAAGCCAAAGCGGGAATATCCGTCATGAGTAATC‑GGGGGTTCGGGGTGTGCTGTTGGATATAAATGCCAATAAACCCTGcgccgcc                         >  1:380096/1‑249 (MQ=255)
                                                                                                                                                                                                                        gTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGGGTCAAAGCCGTCATGAAAAAGC‑GGA‑TT‑GGGAAAGATCAGATCGAATTTGTCGTCAAT‑ATTGGGAAACACGGCCCTTGGC‑ACCCCCGCGCCTTGCCACCGATTTTCCTTTAAAATGCGCTTTTCGGGAGCCCACGCCAAGGCGTGAAAATCCGTCATGATTAATC‑GGGGGTTCCGGGTGTGCTGTTGTATATAACTGCCAATCACCCCCGcgccgccg                        >  1:452989/1‑249 (MQ=255)
                                                                                                                                                                                                                        gTTTTATCAGTTCGTTAACGGCGGT‑GTAAGAGGGGTCCATGCCGTCATGGAAAGGC‑GGG‑TT‑GGGAATGATCAGATCGAATTTGCCGGCAAA‑ATGGGGAAACACGTCCC‑TGGGAAACACCGTGCCCTGCAACCGATTTTTCTTTAAAATGCGCTTTGCGGGAGCCAAAGCCCAGGCGGGGATATACGTCATGATTAATT‑GGGTGGTCGGGGTGTGGTGGTGAATATAACTGCCAATCAccccccgtctcttt                        >  1:278607/1‑241 (MQ=255)
                                                                                                                                                                                                                                      ttAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac          <  1:598440/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac          <  1:422533/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑AGA‑TTGGGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac          <  1:420280/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gt‑gtAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCcg  >  1:488631/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
TTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGT‑GTAAGCGGTGTCAATGCCGTCATGAAAAGGC‑GGA‑TT‑GGAAATGATCAGATCGAATTTGCCGTCAAT‑ATGGGAAAACACGTCAC‑TGGCAACCACCGTGCCTTGCAACCGATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATC‑GGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCG  >  1/1984852‑1985340

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