Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913115,7970CA26.1% 68.9 / 12.7 46E435* (GAA→TAA) hofBT2SSE family protein
Reads supporting (aligned to +/- strand):  ref base C (28/6);  new base A (0/12);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 4.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGC  >  NC_000913/115550‑116045
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
tATCTCCACGCCTCTGAGCGCAATTGATGGCGCTTCGCTTCAGGCTTAGTGCGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCAATAGCATCCCCTCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCCGAGTTGCTTACTCGCCATGCGGCATCACCAGTACGCGGATTCACTCTTCAAAGGTGGTTCAGCCTTGCTCTACGGCAAGGCAGCCGTTTTcca                                                                                                                                                                                                                                                          >  2:420761/1‑248 (MQ=255)
           cACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATCGCCTCCCATCTTGCGCATTGCAGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGcc                                                                                                                                                                                                                                              >  2:257449/1‑250 (MQ=255)
             cTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCAt                                                                                                                                                                                                                                             >  1:281386/1‑250 (MQ=255)
              tgggcAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCAGCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGGTGCTTACTCGCCATGCGGCATCCCCAGTTCGCGGATTAACTCTTCAAAGGTTGTTTAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATa                                                                                                                                                                                                                                            <  1:312944/247‑1 (MQ=255)
                   aGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                                                                                                                                                                                                     >  2:483513/1‑251 (MQ=255)
                     aaTTGATGGCGATTCGCTTCAGGCTTAGTGGGGTGACCATCTGTTGATGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCACATAGCATACCATCTTGCGCATTGCCGTTGCAGGCTATGCCATGCCAGCGCAAGCGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTACTTATTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCCGCAGTTTTCCAAAAGCGTACGCATACCCGCCTg                                                                                                                                                                                                                                   >  2:434777/1‑251 (MQ=255)
                                  tcgcttacggcttagacggcaaaCCATTTGTTGCAGTAGTGCCATAAGCAATAAAGTACGGCTCTCCACCACAAGAATCCCATCTTGCGCATTGCCGTCGCCGCTTATGCCATGCCAGCGCCAGCGTTGCTTACTAGCCATGCGGCCTACCCAGTACGCGCAATAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTAGGCCTACCCGCCTGTCGTGCGTGCGtt                                                                                                                                                                                                                      >  2:312224/22‑251 (MQ=255)
                                        cAGGCTTAGTGGGGTACCCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCACCCCATCTGGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCAATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCTTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                                                                                                                                                                                                <  1:112369/251‑1 (MQ=255)
                                            cTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAtt                                                                                                                                                                                                            >  1:354622/1‑251 (MQ=255)
                                                gTGGGTTAACCATCTTGTGCTGTGGTGCCAAAAGCAGAAAAGTACGGCTCTCGCCCCAAAGCACCCCATCTTGCGCTTTGCCTCCGCCGGTTATCCCATCCCGGCGCCAGGGTTGCTTACTCCCCATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAAGGTGTTTAAGCCTTGCTCAACGGCCGGGCTGCCGTTTTCAAAAGGCGTTCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                                                                                                                                                                                                        <  1:388160/251‑1 (MQ=255)
                                                    ggTAACCCATCTGTTGCTGTAGTGCCATAGGAAGTAAATGACGGCTCTCTGCCCAAAGCACCCCATCTTGCGCATTGCCGTCCCGGGTTATGCCATGCCAGCGCCAGAGTTGTTTACCCGCCATTCGGCACCCCCATTACGCGGATTAACTTTTAAAAGGTGGTTAAGCTTTGCTCAACGGCCTGGCAGCCGTTTAAAAAAAGCGTTCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTc                                                                                                                                                                                                    <  1:57561/251‑1 (MQ=255)
                                                        aCCATCTGTTGCTGTAGGGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTAACAGTCGGTAt                                                                                                                                                                                               >  2:101715/1‑251 (MQ=255)
                                                         ccATCTGTTGCTGTAGTGCCATAAGCAGTACCGTACGGCTCTCTGCCAATAGCATCACATCTTGAGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCACCAGTACGCGGATCAACTCTTCAACGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCCAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                                                                                                                                                                                                        >  2:300967/1‑242 (MQ=255)
                                                                                  gCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGAc                                                                                                                                                                      >  1:405426/1‑251 (MQ=255)
                                                                                     gTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGc                                                                                                                                                                                                                                      >  1:329818/1‑185 (MQ=255)
                                                                                     gTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCAGCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGc                                                                                                                                                                                                                                      <  2:329818/185‑1 (MQ=255)
                                                                                              ggCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATggg                                                                                                                                                          >  2:183079/1‑251 (MQ=255)
                                                                                                             aGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGACAGCGCCAGAGCTGCTTACTCGCAATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTACAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt                                                                                                                                             >  1:437253/1‑250 (MQ=255)
                                                                                                                aTCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAgg                                                                                                                                         >  2:418716/1‑251 (MQ=255)
                                                                                                                 ccccATCTTGCGCATGGCCGTCGCCGGTGATGCCATGCCAGCGCAAGTGTTGCTTACCCCCCAGGCGGCATCCCCAGTCCGCGGATAACCTCTTCAAGGGTGTTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGc                                                                                                                                        <  1:53555/250‑1 (MQ=255)
                                                                                                                            gcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCAATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAGTAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCAt                                                                                                                             >  1:437962/1‑251 (MQ=255)
                                                                                                                                    gTCGCCGGTTATGCCATGCCAGCGACAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACTCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGTAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCg                                                                                                                     >  2:296300/1‑251 (MQ=255)
                                                                                                                                          gtttATGCCATGCCAGCCCCAGAGTGGTTTACCCGCCATGCGGCACCCCCAGTACGCGGATAACCTCTTAAAGGGTGGTAATGCTTTTCTCAACGCCCAGTCTGCCTTTGTCAAAAACCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAg                                                                                                               <  1:326723/249‑1 (MQ=255)
                                                                                                                                          ggatATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGCTTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGCCCAGGCAGCCGTTTTCAAACAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTCGCGGGAATAGGGTGACGAATGAACGGCGTTCAGGGCAGAACTTTAAATAAAGCCGTACTAccctacacaccgaggc                                                                                                                 >  2:9189/4‑235 (MQ=255)
                                                                                                                                             tatcctttcccaGCCCCAGTGTTCCTTACTCCCCAGGCGTCATCCCCAGTCCGCGATTTACCTTTTCAAGGTTGTTTAGGCTTTTCTCTACGCCCAGGCAGCCTTTTACAAAAAGCGTACACTTCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCa                                                                                                             <  1:63954/241‑1 (MQ=255)
                                                                                                                                             tATGCCATGCCAGCGCCAGAGTTTCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCCGCCGTTTTCACCAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCCACGTCGGTATTAGCGGAAATCAGCTTACGAATGACCGGCGTTATGGGCAGAACGTCAAATAAGGCCGTACGGCCATAAAAACCGGGGGAGCa                                                                                                             >  2:252350/1‑250 (MQ=255)
                                                                                                                                              atgccatgccAGCGCCAGATTTGCTTACTCGCCAGGCGGCACCCCCAGTCCGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTCCGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCAGCTCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTAGGGGCAGATCTTGAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCaa                                                                                                            <  1:199968/250‑1 (MQ=255)
                                                                                                                                               tgccatgccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATg                                                                                                          >  1:10307/1‑251 (MQ=255)
                                                                                                                                                        ggcgcaagggTTGCTTACTCCCCAGGCGGCACCCCCAGTCCGCGGATTAACCCTTAAAGGGTGTTTAAGCTTTGCTCTCCGCCCAGGCGCCCTTTTTCAAAAACCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAAcc                                                                                                 <  1:483513/242‑1 (MQ=255)
                                                                                                                                                        aGCGCCAGAGTTGCTTCCGCGCCATGCGGCACCCCCAGTACGCGGATTACCTCTTAAAAGTTGTTTAAGCCTTGCTCTAGGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCTGTACGACCATAAAAACCGTGGTAGCAATGTACACAAcc                                                                                                 <  1:368991/251‑1 (MQ=255)
                                                                                                                                                             cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCggggg                                                                                            >  2:159806/1‑249 (MQ=255)
                                                                                                                                                                            cccaggcgccacccccAGTACGCGGATAAACTCTTAAAAGGTGTTTAAGCCTTGCTCTACGCCCAGGCACCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCa                                                                             <  1:452285/239‑1 (MQ=255)
                                                                                                                                                                                  gggcccccccAGTACGCGGATTAACTCTAAAAAGGTGTTAAAGCCTTGCTCTGCGGCCGGGCGGCGGTTTTAAAAAGGCGAACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAgcggcg                                                                       <  1:345733/245‑1 (MQ=255)
                                                                                                                                                                                               ccgcgAATTACCTCTTAAGAGGGGGTTAAGCCTTGCTCTCCGGCCGGGCGGCCGTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATACGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACAtt                                                          <  1:183079/250‑1 (MQ=255)
                                                                                                                                                                                                  cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGCACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGACTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAATCCGAGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCCTACATTGTc                                                       >  1:25556/1‑251 (MQ=255)
                                                                                                                                                                                                     aTTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGCGGGGCAGCGGCGATGGCCATACATTGTCTgg                                                    >  1:9863/1‑251 (MQ=255)
                                                                                                                                                                                                              caaaaGGTGGTAAACCCTTGCTCTCCGGCCGGGCGGCCGTTTTAAAAAGGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGAt                                           <  1:123620/248‑1 (MQ=255)
                                                                                                                                                                                                                aagggtggtTAAGCCTTGCTCTACGGCCAGGCAGCCGATTTCAAAAAGCGTACGCATACCAGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGGCGGGATTAGCGGAAATACGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAGAAGGCCGGACGACCAGAACAACCGTTGTAGCAATGTACACAACCGGGTCCCTGCCAGGTGGGAAGCGGCGATGGTCATTCATTGTCTGGAATGTGGTTgg                                         >  2:258632/4‑251 (MQ=255)
                                                                                                                                                                                                                  agtgggtTAATCCTTGCTCCCCGCCCAGGCGGCGTTTTTAAAAAAGTGTACGAATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTGATGGGCAGAACTTCAAAGAAGGCCGTACGACCATAAAAACCGGGGTAGCAATGTACACAACCGGGGGCCTGCCAGTGGGGCATCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGc                                       <  1:147204/247‑1 (MQ=255)
                                                                                                                                                                                                                     tggtTACGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGGGGATGGGCTcc                                    >  2:180975/1‑251 (MQ=255)
                                                                                                                                                                                                                       gttagccctgcctctGCGGCCAGGCAGCTGTATAAAAAAAGCGCACGCCTACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTcccc                                  <  2:40561/240‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                              >  2:244481/1‑251 (MQ=255)
                                                                                                                                                                                                                           aGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTtgc                              >  2:158920/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAt                 >  1:54390/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cccgtgctaGCCGTTTTAAAAAAGCGTCCGAATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaa                  <  1:44070/242‑1 (MQ=255)
                                                                                                                                                                                                                                          ccAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAAtgt               >  2:55283/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACaaagctcgga          >  2:44070/1‑242 (MQ=255)
                                                                                                                                                                                                                                                        acAAAAAGCGTCCCATACCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGc  <  1:215468/249‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
TATCTCCTCGCCACTGCGCAGAATTGATGGCGATTCGCTTCAGGCTTAGTGGGGTAACCATCTGTTGCTGTAGTGCCATAAGCAGTAAAGTACGGCTCTCTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGC  >  NC_000913/115550‑116045

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