Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,830,3250GA29.5% 64.9 / 12.9 44E461K (GAA→AAA) xanPxanthine:H(+) symporter XanP
Reads supporting (aligned to +/- strand):  ref base G (20/11);  new base A (0/13);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 4.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGGTTCTGGGCGGCGCAACGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATC  >  NC_000913/3830075‑3830561
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
aCCGGTTCTGGGCGGCGCAACGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAg                                                                                                                                                                                                                                                <  1:257802/251‑1 (MQ=255)
  cGGTTCTGGGCGGCGCAACGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGaa                                                                                                                                                                                                                                              <  1:117939/251‑1 (MQ=255)
        tGGGCGGCGCAACGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑CGAATGGCTGAAAAACCTGCTCTCCTCCGGGAACGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGGATCTGATTTTCCCACCAGACAAA‑c                                                                                                                                                                                                                                         >  2:471137/1‑249 (MQ=255)
                  aaCGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGTCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTg                                                                                                                                                                                                                            <  1:322860/251‑1 (MQ=255)
                   tcgcgtGTAATGTTTGGAGCCATCGCCGCCTCCGGTTTGCTCATCGTTTCTCGAGAGCCGCTGAACCATCGGGCGATTCTGTCTATCGCGCTGTCGCTGTCGGTTGG‑TCTGGGCGTGTCTCGGCAGCCGCTGATTTTGCAGTTTGCCTC‑TGAATGGTTGAAAAACCTTCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACt                                                                                                                                                                                                                             <  2:284434/244‑1 (MQ=255)
                   aCGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTTCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGGATTACTGCCAT‑CGTGCTGAATCTGGTTTTCCCACCCGAAAAA‑CAGTAATT‑CACTgg                                                                                                                                                                                                                           >  2:298787/1‑251 (MQ=255)
                        ngTAATGTTTGGCACCATCGCCGACTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGGGTGTCTCAGCAGGCGCTGATTTTGGAGTTTGCCCC‑TGAATGGCTGAAAAAACTGCTCTCCTCCGGGATCGCCGGGGGCGGTATTAATGCCAT‑AGTGGTGAATTTGGTTTTCCCAACAGAAAAA‑AAGTAATT‑TACTggggggg                                                                                                                                                                                                                      >  2:252074/2‑246 (MQ=255)
                                 tGGCACCATCTCCGCCTCTGGTGTGCGTTTCGTATCTCCGGAGCCGCGGAACCGTCGGTCGATTCTGATTATCGTGTTGTTGCCGGCGGTTTG‑TCTGGGCGTGTCTCTGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTGCTCCGGGAGCGCCGCGGTCGGTATTACTGCCCT‑TGTGCTGAATCTGATTTTCCCCCCAGAAAAA‑CAGTCATT‑CACTGGCGTCGGTATGATa                                                                                                                                                                                                              <  2:155496/250‑1 (MQ=255)
                                 tGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGGATTTCTGCCAT‑CGGGCTGAATCTGATTTTCCCACCCGAAAAA‑CAGGAATT‑CACTGGCGGCGGTAAGATa                                                                                                                                                                                                              >  2:177833/1‑250 (MQ=255)
                                                  ccGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑GCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGATCTCCTCCGGGATCGCCGCGGGCGGGATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGGAAGATAATTAACGCCGCAAtctc                                                                                                                                                                                             >  1:103674/1‑250 (MQ=255)
                                                             ttCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGTTTATCGCGTTGTCGCTGGCGGTGGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTCCAGTTTGCCCC‑TGAATGGCTAAAAACCCTGCTCTCCTCCGGGATCGCCGCGGGCGTTATTACTCCCAT‑CTTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGGGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAAc                                                                                                                                                                                  <  1:257673/249‑1 (MQ=255)
                                                               cGTTTCTCGGGGGCCGCTGAACCTTCGGGCAATTCTTGTTATCGCGCTGGCGTTGGCGTTTGT‑TCTGGGCGTGTCTCAGCAGCCGCTGATTGTCAAGTTTGCCCC‑GGAATGGCTGAAAAACCTGCTCCCCTCGGGGATCGCCGGGGGCGGTATTACTGCCTT‑CGTGCGGACGCTGGTTTTCCCCCCAGAAAAA‑CAGTATTT‑CACTGGCGGGGGTAGGATATTTACCGCCGCAGTCTCACCTCTTTAACCa                                                                                                                                                                                <  2:287687/250‑1 (MQ=255)
                                                                  ttCTCGTGCTCCGCTGAGCCGTCGGGGTATTCTGATCAGCACGGTGTCTCTGTCGTTTGG‑TTTTGGCGTTTCTCAGCATCCGCTGTTTTTGCAGTTTGCCCC‑GGATTGGTTGAAAACCCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CTTGCTGAATCAGATTTTCCCCCCAGAAAAA‑CAGTAATT‑CGCTGGCGGCGGAAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTc                                                                                                                                                                             <  2:80722/250‑1 (MQ=255)
                                                                         gggCCGCTGACCCGTCGGGCGTTTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGTTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAACCCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCCCCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTg                                                                                                                                                                     <  1:316847/249‑1 (MQ=255)
                                                                                 gacccgtggtgctgttCTGTTTCTCGCGGTGTCTCTGG‑GGTTGGGTCGGGGCGTGTCTACGCCGCCTCTGATTTTGGA‑TTTGCCCCGTGATGGGCTGAAAACCCTGCTCTCCCCAGGGATCCCCGCGGGGGTTATTACGCCCAT‑GGTGTTAAATCGGTTTTCCCCCCCAAAAAAA‑CTGTAATT‑CACTGGCGGCGGTAAGATAATGCCCGCCGCTGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTgc                                                                                                                                                             <  2:41132/237‑1 (MQ=255)
                                                                                     cGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCCCACCTCTTTAACCATTTCCGCCTTGAGGATTGCGCGt                                                                                                                                                         >  1:527684/1‑251 (MQ=255)
                                                                                          ggCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTTCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGGATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTTGCGGCGGTAAGATAATTACCCCCGCAGGCTCACCTCTTTTACCATTTCCGCCTTGAGGATTGCGCGTAAATc                                                                                                                                                    >  2:76649/1‑251 (MQ=255)
                                                                                                         cgcgCTGTCGCGGGCGGTTGG‑TCTGGGCGTGTCCCAGCCGCCGTTGTTTTTGCAGTTTGCCCC‑TGAATGGTTGAAAACCCTGCTCTCCTCGGGGATCGCCGCGGGGGTTATTCCTGCCAT‑CGTGCTAAATCTGATTTCCCCCCCAAAAAAA‑CATTAATT‑CACGGGCGGCGGTAAAAAAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCtt                                                                                                                                     <  1:177833/251‑1 (MQ=255)
                                                                                                               gTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGGCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTTTGTGc                                                                                                                               >  2:79075/1‑251 (MQ=255)
                                                                                                                   gtGGCGTTTGG‑TCTGGGCGTGTCTCAGCACCCGCTGATTTTGAGGTTTGCCCG‑TGAATGGCTAAAAACCCTCCTCCCCCCCGGGACCGCCGCGGGCGGTATTCCTGCCAT‑CGTGCTGAATCTGTTTTCCCCCCCAAAAAAA‑CAGTAATT‑CACGGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCAttt                                                                                                                           <  1:79075/250‑1 (MQ=255)
                                                                                                                                gtgggctttcctcAGCGGCCGCTGTTTTTGCAGTTTCCCCC‑TGATTGGCTAAAACCCCTCCTCTCCTCCGGGCCCCCCGCGGGGGTTATTCCTCCCTT‑CGTCCTGAATCTGTTTTCCCCCCCAAAAAAA‑CAGTAATT‑CCCTGGCGGCGTAAAGATATTTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGaa                                                                                                               <  1:157273/241‑1 (MQ=255)
                                                                                                                                           tcAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAAt                                                                                                    >  2:78102/1‑251 (MQ=255)
                                                                                                                                           cccgcacccgtgaaTTTTGAATTTTGCCCC‑TGAATGGTTAAAAACCCTGCTCTCCCCGGGGATCGCCGCGGGGGGTATTACTGCCAT‑CGTGCTGAATCTAATTTCCCCCCCAAAAAAA‑CAGAAATT‑CACGGGCGGCGGTAAGATAATTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAAt                                                                                                    <  1:403492/238‑1 (MQ=255)
                                                                                                                                                   cGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACCGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTggg                                                                                            >  1:286857/1‑251 (MQ=255)
                                                                                                                                                            tGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCtt                                                                                   >  1:468763/1‑251 (MQ=255)
                                                                                                                                                             gCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTtc                                                                                  >  1:206028/1‑251 (MQ=255)
                                                                                                                                                                               ggCTGAAAACCCTGCTCCCCCCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATTTGATTTTCCCACCAAAAAAA‑CAGTAATT‑CCCTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGc                                                                  <  1:148905/250‑1 (MQ=255)
                                                                                                                                                                               ggCTAAAAACCCTGCTCCCCCCGGGGATCGCCGCGGGCGGTTTTACCGCCAT‑CGTGCTAAATCTGTTTTTCCCCCAAAAAAAA‑CAGTATTT‑CTCTGGCGGCGGTAAGAAATTTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGTAAGCTGCTTCTCTACATTCTCATCGct                                                                 <  1:196497/251‑1 (MQ=255)
                                                                                                                                                                                   gAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGt                                                             >  2:373042/1‑251 (MQ=255)
                                                                                                                                                                                    aaaaccctcctctcccccggtgccgccggGGGCGTTATT‑CTGCCATTCGGGCTGAATCGGATTTTCCCACCAAAAAAA‑CAGTATTT‑CCCTGGCGGCGGTAAAATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTGTTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGtt                                                            <  1:298787/228‑1 (MQ=255)
                                                                                                                                                                                        aCCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGACAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTg                                                        >  2:162617/1‑251 (MQ=255)
                                                                                                                                                                                                ccccccGGGATCGCCGGGGGCGGTATTCCTCCCTT‑GGTGCTGAATTTGTTTTCCCCCCCAAAAAAA‑CAGTAATT‑CACTGGGGGCGGAAAGATAATTCCCGCCGCAGTCTCCCCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGc                                                <  1:556770/247‑1 (MQ=255)
                                                                                                                                                                                                  ctccGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTg                                              >  1:539201/1‑251 (MQ=255)
                                                                                                                                                                                                    ccgggtccgccggGGGCGGTTTTCCTGCCAT‑CTTGCTAAATCTGATTTCCCCCCCAAAAAAA‑CATTAATT‑CCCTGGCGGGGGTAAAATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGc                                            <  1:252074/244‑1 (MQ=255)
                                                                                                                                                                                                    ccGGGATCGCCGGGGGCGTTATTACTGCCAT‑CGTGCTGAATCTGTTTTTCCCACCAGAAAAC‑CAGTAATT‑CACTGGCGGCGGTAAAATAATTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGc                                            <  1:172320/251‑1 (MQ=255)
                                                                                                                                                                                                       ggaccgccggGGGCGTTATTACTGCCAT‑CGGGCTGAATCTGATTTCCCCCCAAGAAAAA‑AAGTATTT‑CACTGGCGGGGGTAAGAAAATTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCctc                                         <  1:471137/247‑1 (MQ=255)
                                                                                                                                                                                                       ggATCGCCGCGGGCGGTATTCCTCCCTT‑CATTGGAAATGTGATTTCCCCCCCAGAAAAC‑CA‑TAATTTCTCGGGCGGCGGAAAGATAATTACCGCCGCAGTCTCCCCTCTTTAACCATGCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCctc                                         <  1:198550/251‑1 (MQ=255)
                                                                                                                                                                                                           cgccgcGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTAtt                                     >  1:552069/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggcggTATTACTGCCTT‑CTGGCTGAATCTG‑TTTTCCCCCCAAAAAAAACAGTAATT‑CCCTGGGGGGGGAAAGATAATTCCCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTtcttct                              <  1:201276/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTCTGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCa                           >  2:410269/1‑250 (MQ=255)
                                                                                                                                                                                                                         ttACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAAcc                       >  2:62167/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcccAT‑GGTGCTAAATCTTATTTCCCCACCAAAAAAA‑CAGTATTC‑CACGGGCGGCGTTAAGATATTTCCCGCCGCAGTCTCACCTCTTTACCCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCt                   <  1:172668/249‑1 (MQ=255)
                                                                                                                                                                                                                              gCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCTg                  >  2:452880/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cggtgttaattcttttttccccACCAAAAACA‑CATTATTT‑CACTGGCGCGGGTAAGATAATTCCCGCCGTAGCCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGa              <  1:39330/232‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTCCCCCCAGAAAAA‑CAGTAATT‑CACTGGCCGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTTGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATc  >  2:190551/1‑249 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
ACCGGTTCTGGGCGGCGCAACGCTTGTAATGTTTGGCACCATCGCCGCCTCCGGTGTGCGTATCGTTTCTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGG‑TCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCC‑TGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCAT‑CGTGCTGAATCTGATTTTCCCACCAGAAAAA‑CAGTAATT‑CACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCCGCCTTGAGGATTGCGCGTAAATCGTGCATAACTCCCTTATGTGCATTTCACGGGATGGAAGACCATGAAATTTATTGGGAAGCTGCTTCTCTACATTCTCATCGCTCTGTTAGTGGCGATCGCTGGCCTCTATTTTCTTCTGCAAACCCGCTGGGGAGCAGAACATATC  >  NC_000913/3830075‑3830561

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