Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,1920GC31.6% 77.3 / 21.2 57A293G (GCG→GGG) araGarabinose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (13/26);  new base C (17/1);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.15e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTAAAGACAGTAATGGC  >  NC_000913/1983959‑1984436
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            
tgATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATCTAAACCTGACCGGCGGTGA‑GTTTCGTCCC‑GCC‑CA‑ACATGGCCTTTC‑TTAATTTGCTACGCCCCGCCCCTACCCGAC‑CAAAc                                                                                                                                                                                                                                            >  1:506381/1‑251 (MQ=255)
tgATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAAACTGACCGGCGGGGA‑ATTTCGGCCC‑CCC‑AA‑ACATA‑CCTTTTCTTAATTCGCCACGCCCCCCCCCCACCCGAc‑ccaa                                                                                                                                                                                                                                             >  1:10582/1‑247 (MQ=255)
  aTTACACAACCGCCGAGCACCTGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTTCCGCCC‑CCC‑AA‑ACATG‑CCTTTTATTAAATCGCTACACCCCCCCCCCAcccaaccacaaa                                                                                                                                                                                                                                            >  1:623609/1‑240 (MQ=255)
       acaACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGCGA‑TTTTCGTCCC‑GCCCAA‑AAAT‑‑CCCTTCATTTAATCGCTACACCCCGCCCCCACCAACC‑CAAACAC‑CCca                                                                                                                                                                                                                                     >  1:588299/1‑250 (MQ=255)
         aaCCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GATCATCATGCCTGCGGCAATGGCGGGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGTTGA‑TTTGCGCCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCATAC‑CAAACAG‑CCcaaca                                                                                                                                                                                                                                  <  2:89032/251‑1 (MQ=255)
               cgAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGG‑TTTGCGTCCC‑GCC‑CA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCCGAC‑CAAAAAA‑CCCAACAA‑ttttc                                                                                                                                                                                                                           >  1:201062/1‑249 (MQ=255)
               cgAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATT‑CCCTTCATTAATTTGCTACCCCCCCCCCCCACCAGAC‑CCAACAG‑CCCCACAA‑ATTcc                                                                                                                                                                                                                           >  1:184567/1‑250 (MQ=255)
                gAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGGGG‑TTTTCGGCCC‑GCC‑AA‑ACAAG‑CCTTTCATTAATTCTCTACGCCCCCCCCCCACCCGAC‑CCAACAA‑CCCCACAA‑ATTcac                                                                                                                                                                                                                          >  1:400109/1‑251 (MQ=255)
                  gCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGAATTTG‑GTCCC‑CCC‑AA‑ACATG‑CCTTTCATTAAATCTCTACGCCCCCCCCCTAACAGAC‑CAAACAA‑CCCCACAA‑TTTccccc                                                                                                                                                                                                                        >  1:703020/1‑247 (MQ=255)
                      aTGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGG‑TTTGCGTACC‑GCC‑AA‑AAATA‑CCCTTCCTTAATATGCCACGCCCCCCCCCAACCCACC‑CAAACAGACCCAA‑AA‑TTTcacccatcc                                                                                                                                                                                                                    >  1:377663/1‑246 (MQ=255)
                        gtttgctgctgtcACTGATGGTGATATTGTCGCGAACGGATTGCACGGTAATAATGCCTTCCGCTTGGCGTTCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATTTAGAACGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑CCATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTTACCACttcg                                                                                                                                                                                                                   <  2:666320/238‑4 (MQ=255)
                               tCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGG‑TTTTCGTCCC‑GCC‑AA‑ACATT‑CCTTTCTTTAATTCGCTACGCCCCCCCCCTACCAGAC‑CAAACAG‑CCCCACAA‑TTTTACCACCGCG‑CACCGccc                                                                                                                                                                                                          >  1:161044/1‑250 (MQ=255)
                                             tgttatggTCGCTAACGGAGGGCCCGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCGGCTGCAATGGCGTGGCTCGTGTTACGAATATCGATCGGCTTTTGGCCGAAATAAACCTGACCGTCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTCCGCCCCGCCCCTACCAGTC‑CAAACAG‑CCCAACAA‑TTTCCCCACTGCG‑AACCGCCAGACTTATTGGCGt                                                                                                                                                                                             <  2:872362/243‑1 (MQ=255)
                                                   tGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑AAATG‑CCTTTTATTAAATCGCTACGCCCCCCCCCCCCCCGAA‑CAAACCG‑CCCCACAA‑ATTTACCCCCGCG‑AACCCCCCGACTTATTTGCGTacgcccc                                                                                                                                                                                      >  1:490387/1‑248 (MQ=255)
                                                       gcgAACGGAGTGAACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGGTGGTCGATATAAACCTGACCGGCGGTGG‑TTTTCGTCCC‑CCC‑AA‑AAAAG‑ACTTTCATTAATTCGCTCCGCCCCCCCCCTTCCCGGCACAAACA‑‑CCCAACAA‑ATTTACCACTACG‑AACCGcccaaattatttgccgtcgcacgactg                                                                                                                                                                                  >  1:446081/1‑226 (MQ=255)
                                                        cgAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACCTT‑CCCTTCCTTAAATCGCTACGCCCCCCCCCCTCCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCCCTTCG‑CACCGCCCGACttttttgcggtcgcaacccatg                                                                                                                                                                                 >  1:173799/1‑231 (MQ=255)
                                                             ggAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGATTACGAATATCGCTCGGCTGTTGGTCGATATAAACATGACCGGCGGTGA‑GTTGCGTCCC‑GCC‑CA‑ACATG‑GCTTTTTTTTATTCGCTACGCCCCCCCCCCACCACAA‑CAAACAA‑CCCAACAA‑TTTCACCACCGCG‑AACCGTCACACTTATTGGCGTACGA‑ACGCC‑Tgggggtt                                                                                                                                                                          >  1:243000/1‑246 (MQ=255)
                                                                       ggAATAAGTCCTTCCTCTTTGCCATCTTCCGG‑CATAAGCATCTGGCCTGCGGCAATGGCGTGCCTCGTTTAACGATTATCGATCGTCGGTTGGTGCATTTAAACCTTACCGCCGGTCA‑TTTCCTTCCC‑GCC‑AA‑ACATC‑CCTTTTATTAATTCTCTCCGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATc                                                                                                                                                                <  2:180352/251‑1 (MQ=255)
                                                                            aatGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTTCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTTCGCCCCCCCCCTACCAGAA‑CAAACAA‑CCCAACAA‑TTTTACCACAGCG‑AACCCCCCCACTTCTTGGCGTACGCCAC‑CC‑TGGTGCTTTCACAGCATccaaa                                                                                                                                                            >  1:446169/1‑246 (MQ=255)
                                                                             atGCCCTCCGTTTTGCGATCTCCCGGGCA‑GGGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTGGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTCCGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACg                                                                                                                                                          <  2:490387/251‑1 (MQ=255)
                                                                                    ccGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGGCCC‑GCC‑CA‑ACATG‑CCTTTCATTAATTTGCTACGCCCCCCCCCTACCACAC‑CAAAAAA‑CCCAACAAATTTC‑CCACCTCG‑AACCGCCCCACTTATTTGCGTACGC‑ACGCC‑TGGTGcctttccagactccaaacggaggagc                                                                                                                                                   >  1:343336/1‑226 (MQ=255)
                                                                                       cTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTTCGTCCC‑GCC‑AA‑ACAAG‑CCCTTTATTAATTTGCTACGCCCCCCCCCTACCCGAC‑CCAACAG‑CCCCACAA‑TTTTACCACTGCG‑AAACGCCCGACCTATTGGCGTACGC‑ACGCC‑CGGGGCTTTtacaaccacaaaaactaagccgtcc                                                                                                                                                >  1:296859/1‑230 (MQ=255)
                                                                                       cTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTTCGGCCC‑GCC‑AA‑AAAAG‑CCTTTTATTTATTTGCTAAGCCCCGCCCCCAACAGAA‑CAAAAAA‑CCCCACAA‑ATTCAACACATCG‑AAACGCCCGAATTATTTGCGTAAAC‑AAGCC‑CGGGGGTTTCACAACATCAAggagtagacggtcc                                                                                                                                                >  1:445430/1‑238 (MQ=255)
                                                                                          tGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGctcctc                                                                                                                                              <  1:327442/250‑1 (MQ=255)
                                                                                                  tCCGGGCA‑GAGCATCATGCCTGCGGCAAAGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATa                                                                                                                                        <  1:375093/248‑1 (MQ=255)
                                                                                                   ccGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCCCCCCCACCCGAC‑CAAACAG‑CCCAACAA‑TTTTACCACTGCG‑AACCGCCCGACCTATTTGGGTACGC‑ACGCC‑CGGTGCTTTTACAGCATCAAAACGTAAGCGCTCCTCGCCCTAACt                                                                                                                                     >  1:75561/1‑250 (MQ=255)
                                                                                                      gggcagagatacatgctggCGGAAATGGCGGGGCACGTATTACGAATAGCGATCGGCTGTTGGGCAATAAAAACCTGACCGGGGGTGA‑TTTGCGTCCC‑GCC‑AA‑AAACG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTAcg                                                                                                                                  <  2:298258/233‑1 (MQ=255)
                                                                                                                  atGCCGGCGGCAAGGGCTTGGCTCGGTTTACTAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTGCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGca                                                                                                                     <  1:312993/251‑2 (MQ=255)
                                                                                                                           gCAATGGCGTGGCTCGGTTTACGAATATCGTTCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGCCCC‑GCC‑AA‑ACTTG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑c                                                                                                           <  1:604147/251‑1 (MQ=255)
                                                                                                                              atgtcGTGGCCCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGTTGA‑TTTGCGCCCC‑CCC‑AA‑ACATG‑CCTTTCATTAATTCGCTCCGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGcc                                                                                                        <  1:415668/247‑1 (MQ=255)
                                                                                                                                         cGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCCCCCCTACCAGAC‑CAAACAA‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCCGACTTATTTGCGGTCGC‑ACGGC‑TGGTGCTTTTACAGCAACAAGACGGAGGCGCTTCTCTCCCTTACTAAGCGG‑TTGGCAGCCGTAGATATAGGCCAA‑GTCGCGcc                                                                                             >  1:170550/1‑251 (MQ=255)
                                                                                                                                               aCGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGG‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCCCCCCTACCAGAA‑CAAAAAG‑CCCAACAA‑TTTCACCCCTGCG‑GACCGCCAGGCTTATTGGCGTACGC‑ACGCC‑TGGGGGTTTCCCAGCATCAAGACGTAAGCGCTCCTCGCCATAACTACGCGGGTTGC‑AGCCGTTGATAT‑CGCCAAAGTcaccgccgaccc                                                                                       >  1:252538/1‑240 (MQ=255)
                                                                                                                                                    atccgatcgTCTTTTGGTCGATATAACCCTGCCCGGCGGTGA‑TTTGCGCCCC‑GCC‑AA‑ACATG‑CCTTTCATAAATCCGCTCCCCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTCCACCCCTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGcc                                                                                  <  2:75561/248‑1 (MQ=255)
                                                                                                                                                     accgttcggtgttgggtcaaaataACCCTGACCGGCGGGGA‑TTTGCGCCCC‑GCC‑AA‑AAATG‑CCTTACTTAAATTCGCTACGCCCCGCCCCTCCCAGAC‑CAAACAG‑CCCAACAA‑TTCCCCCCCTGCG‑AACCGCCAGACTTATTGGCGTACGC‑CCGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCt                                                                                 <  2:257435/230‑1 (MQ=255)
                                                                                                                                                           cgGCTGTTGGTCGATATAACCCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCAcc                                                                            <  1:539748/250‑1 (MQ=255)
                                                                                                                                                             gctgtggttcgatgaaaccctcacgggcggTGA‑TTTGC‑TCCCATCC‑AA‑ACATT‑CCTTTCATTAATTCGCTCCCCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACGTCG‑AACCGCCAGACTTATTGGCGTACTC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAgc                                                                         <  2:795684/226‑1 (MQ=255)
                                                                                                                                                             gCTGTTGGTCGATATAAACCTGACCGGCGGTTA‑TTTGCGTCCC‑GCC‑AA‑ACAGG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGACCGACCATCGCCTGCACCAgc                                                                         <  1:799319/251‑1 (MQ=255)
                                                                                                                                                               tGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑CA‑ACATG‑CCTTTCATTAATTCGCTACACCCCGCCCCTACCAAAC‑CAAAAAA‑CACCACAA‑ATTCACCACCGCG‑AAACCCCCGACTTTTTAGCGTACCA‑ACCCCCTGG‑GCTTTCACAGCATCCAGAAGTAAGCGCCCCTCGCCACAACAACGCGG‑GTGCCAGACGTAGATAT‑CCCCAATGTGCCGCCCGACCACCGCCTGccaccgcac                                                                       >  1:836481/1‑243 (MQ=255)
                                                                                                                                                               tGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGACCC‑CCC‑AA‑AAATG‑CCTTTCATTAATTAGCTACGCCCCGCCCCTAACCGAC‑CAAAAAGGCCCACAAA‑TTTCAC‑ACTGCGGAACC‑CCAGACTTAATGGCGTACGC‑ACGCC‑CGGGGCTTTCACAACAACCAGACCTAGGCGCTCCACCGCAACACCACGCGG‑TTGCCAGCCGTAGCAAT‑AGCCCATGTTGCGCCCggacaacacccgcaacagcac                                                                       >  1:836458/1‑231 (MQ=255)
                                                                                                                                                                             acccggcccgggggggA‑TTTGCGTCCC‑GCC‑AA‑ACAGC‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGAG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACct                                                         <  1:671645/236‑1 (MQ=255)
                                                                                                                                                                             aaCCTGACCGGCGGTGA‑TTTCCGCCCC‑GCC‑AA‑CCATG‑CCTTCCATTATTTCCCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCCCTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACct                                                         <  1:640170/251‑1 (MQ=255)
                                                                                                                                                                              aCCTGACCGGCGGTGA‑TTTGAGTCCC‑GCA‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACGGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑GGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACctg                                                        <  1:619794/251‑1 (MQ=255)
                                                                                                                                                                               ccTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑TCATG‑CCTTTCATTAATTCGCTACCCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACctgc                                                       <  1:403510/251‑1 (MQ=255)
                                                                                                                                                                                        gggTGA‑TTTGCTCCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTCCGCCCCGCCCCTACCAGAC‑AAAACAG‑CCCAACAA‑TTTCCCCACTGCT‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCg                                              <  1:794931/250‑1 (MQ=255)
                                                                                                                                                                                        gggTGA‑TTTGCGTCCC‑‑CC‑AACACATG‑CCTTCCATTAATGCGCTCCGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑GTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCg                                              <  2:645381/250‑1 (MQ=255)
                                                                                                                                                                                         ggTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCgg                                             <  1:510655/251‑1 (MQ=255)
                                                                                                                                                                                                 tGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTAACAGCCTCAAGACGTAAGCGCTCCTCGCCCTAACTACGCGG‑TTGCCAGCCGTAGATAT‑AGCCAATGTCGCGCCCGACCATCGCCCGCCCCAGCCGGTCGGGGTCAACCTGCTCCATACCGGTAAGGGt                                      >  2:275199/1‑251 (MQ=255)
                                                                                                                                                                                                   cGTCCC‑ACC‑AA‑ACATC‑CCTTACATAAATTCGCAACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑GGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGttt                                    <  2:201062/251‑1 (MQ=255)
                                                                                                                                                                                                     ggcgcccctccaaacactcttttCTTTAATTCTCTAGCCCCCCCCCCTCCCAGCC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑ACCCGCCAGACTTATGGGCGTACGC‑ACGCC‑TGGTGCTTCCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTTCATATCGGTAAAgg                                       <  2:252538/230‑1 (MQ=255)
                                                                                                                                                                                                      ccc‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑GCGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCAACAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTg                                  >  1:415196/1‑250 (MQ=255)
                                                                                                                                                                                                                  cATG‑CCTTTCATTAATTCACTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTATGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCACCCTGCTGCATCTCGGGAAGGGGTTTGACATAACGTc                        >  2:274602/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTTAATCCCCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTaaa              <  1:819715/251‑1 (MQ=255)
                                                                                                                                                                                                                              ttAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTAGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTTCATATCGGTAAAGGGTTTGACATAACGTCCATCTTTaa               >  2:585361/1‑249 (MQ=255)
                                                                                                                                                                                                                               taatcctcTACGCCCCGCCCCTCCCGGAC‑CAACCAG‑CCCAACAG‑TTCCACCCCTGCG‑AACCGCCAGACTTATTGGCGGACGC‑ACGCC‑TGGTGCTTTCCCAGCATCAAGACGTAGTCCCTCCTCTCCATAACTACTCGG‑GTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTACTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTAAAg             <  2:69535/243‑1 (MQ=255)
                                                                                                                                                                                                                               aatttcccccccccccgccctTACCAAGC‑CAACCAC‑CCAAACAA‑TTACACCACTGCA‑GACAGCCAGCCTGATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGGCGTATGCGCTCCTCGCCATAACTACGCGG‑GTGCGAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTAAAGa            <  2:588299/239‑1 (MQ=255)
                                                                                                                                                                                                                                aaTTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAAGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATTTCGGGAAAGGTTTTGACATAACTCCCATCTTTAAAGAc           >  2:858418/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAAGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCAGACAATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGGAAAGGTTTTGACATAACGTCCATCTTTAAAGACAGt        >  1:587147/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tACGCCCCGCACCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCCCTGCG‑ACACGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGGGCTTTCACAGCCTCAAGCCGTAAGCGCTCCTCGCCATAAATACGCGG‑TTGCCAGCCGTAGATAA‑CGCCAATGTCGCCCCCGCCCATCGCATGCGCCAGCGAGACGTGGTCAAACTGCTGCATATCAGGCACGGTTTTCACATAACGTCCAACTTGAAAGACCGTAAtt    >  2:213949/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gccccgccccTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGAAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCCCCAGCGCCTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTAAAGACAGAAATGGc  >  2:216939/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                            
TGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGTCGCGAACGGAGTGCACGGGAATAATGCCTTCCGCTTTGCGATCTTCCGGGCA‑GAGCATCATGCCTGCGGCAATGGCGTGGCTCGGTTTACGAATATCGATCGGCTGTTGGTCGATATAAACCTGACCGGCGGTGA‑TTTGCGTCCC‑GCC‑AA‑ACATG‑CCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGAC‑CAAACAG‑CCCAACAA‑TTTCACCACTGCG‑AACCGCCAGACTTATTGGCGTACGC‑ACGCC‑TGGTGCTTTCACAGCATCAAGACGTAGGCGCTCCTCGCCATAACTACGCGG‑TTGCCAGCCGTAGATAT‑CGCCAATGTCGCGCCCGACCATCGCCTGCACCAGCGCGTCGTGGTCAACCTGCTGCATATCGGTAAAGGTTTTGACATAACGTCCATCTTTAAAGACAGTAATGGC  >  NC_000913/1983959‑1984436

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