Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913388,3360AC32.1% 79.2 / 10.1 54N181T (AAC→ACC) tauDtaurine dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base C (1/16);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGC  >  NC_000913/388088‑388581
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
gCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAAcc                                                                                                                                                                                                                                                      <  1:697914/251‑1 (MQ=255)
                     tgacccGCCACCCTGAGGGGCGATTCTGGAAGCAAAAGAGTTCCTTTCGACGGGGGGTGATACGCTCTGGACCAGCGGTA‑TTGGGGCCTATGAGGCGCTCTCTGTTCCC‑TCCCGCCAGCTGC‑TGAGTGGGC‑TGCTGGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAATTCCCGCAAACCCGAGGAGAAACATCACCGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCggtg                                                                                                                                                                                                                                  <  1:852578/245‑1 (MQ=255)
                                 ggggggggCGATGCGGGCAGATAAAGATTTACTTTCGCCGGGCGGTTATCCGCTCTGGACCAGGGGTA‑TTGGGGCTAAGGGGGCG‑TCTTCGTTCCCCTCCCGCCGGCTGC‑TGAGTGGGC‑TGCGGGCGGAGCTGGATTTCCGAAAATCGTTCCCGGAATA‑CAAATCCCTCAAAACCGAGGAGAAACACCAACGCGGGCGGGGGGCGTT‑CGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATc                                                                                                                                                                                                                     <  2:772673/247‑1 (MQ=255)
                                               gggCAGCTAAAGAGTTACCTCCGACCGGCGGTGATACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑GGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTaaa                                                                                                                                                                                                       <  1:319173/250‑1 (MQ=255)
                                                  cagttAAAGAGTTCCCTTCGACCGGCGGTGATCCGCTCGGGACCAGCGGTA‑TTGCGGCCTAGGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGT‑TGAGTGGGT‑TGCGTGCGGAGAATGATTTCCGTAAATCGTTCCCGAAATA‑AAAATACCGCAAAACCGAGGAGGAACATCACCGCGGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAg                                                                                                                                                                                                    <  1:62656/247‑1 (MQ=255)
                                                           agTTACCTTCGACCGGCGGTGATACGCTCTGGACCCGCGGTA‑TTGCGGCCTCTGGGGCGCTCTCTGTTCCC‑TTCCGCCCGCTGC‑TGGGTGGGC‑TGCGTGCGGGGCATGATTTCCGTAAATCGTTCCCCGAATA‑CAAATACCGCAAAACCGAGGAGGGACATCCACGCTGGCGCGGGGCGGG‑CGCGGAAAACCCGCCGGTGCTACATCCGGGGGTGCGAACGCACCCCGTGGGCGGGAAAAAggaccggttt                                                                                                                                                                                           >  1:573997/1‑243 (MQ=255)
                                                            tttacctgcgacaggcgtggaagcGCTCTGGACCAGCG‑TAGTTGCGGCCTATGAGGCGCTCTGTGTTCCC‑TTCCGCC‑TATGCCTAAGTGGGC‑TGCTTGCTGAGCATGTTATCCGTAAATCGTTCCCGGAATA‑CAAATCCCGCAAAACCGAGGAGGAACATAATCTCTGGCGCGGGGCGGT‑CGCGAAAAGCCCGCCGTTGCTTCAGCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGttt                                                                                                                                                                                           <  1:251496/227‑1 (MQ=255)
                                                                     gACCGGCGGTGATACGTTCTGGACCAGCGGTT‑TTGCGGCCTATGATGCGTTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCTTGCGGAGCATAATTTCCGAAAATCTTTCCCGGAATA‑CAAATACCGCAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAg                                                                                                                                                                                 <  1:36340/251‑1 (MQ=255)
                                                                        cggcggTGATACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCGCTGTTCCC‑TCCCGCCAGCTGC‑TGGGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑AAAATACCGCAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCt                                                                                                                                                                              <  1:803366/251‑1 (MQ=255)
                                                                            ggtgatcctctctgaccctgctgTA‑TTGCGTCTTATGGGGCTCTCTCTGTCCCC‑TTCCGCCAGCTGC‑TGAGGGGGC‑GGCTTGCGGAGCATCATTTCCGAAAATCGTCCCCGTAAAA‑CAAATACCGCAAAACTGAGGAGAAACATCACCCCGGTCGCGGGGCGTT‑CGCGAAAAACCCGCCGTTGCTACATCCGGGGGTGCGAACGCATCCGGTGAGCGGTAAACAGTCGCTGTTTGTGAATGAAGGCTTtac                                                                                                                                                                          <  1:266618/229‑1 (MQ=255)
                                                                                atccGCTCGGGACCGGCGGTA‑GGGCGGCTTAGGAGGCGCTCTCTGTCCCC‑TCCCGCCATCGGC‑TGAG‑GGGCGTGCGGGCGGAGCATGATTCCCGAAAATCGTTCCCGAAATA‑AAAATACCGCAAAACCGAGGGGGAACATCAACGCGGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcg                                                                                                                                                                      <  1:45459/248‑1 (MQ=255)
                                                                                 tACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcgc                                                                                                                                                                     >  1:180481/1‑251 (MQ=255)
                                                                                          taCCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAAttgttg                                                                                                                                                            >  1:266449/2‑251 (MQ=255)
                                                                                                             cTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGTAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAagagg                                                                                                                                          >  1:101345/1‑250 (MQ=255)
                                                                                                                  gggCGCTCTCTGTCCCC‑TCCCGCCAGCTGC‑TGAGGGGGC‑TGCGTGCGGAGCATGATTCCCGAAAATCGTCCCCGGAATA‑CAAATCCCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGaa                                                                                                                                     <  1:191226/250‑1 (MQ=255)
                                                                                                                     cgcTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAAACCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTCCGCGAATTGGTGATGTGAGCGGGCAAGGGAGCGAAGcc                                                                                                                                  >  2:253447/1‑251 (MQ=255)
                                                                                                                           ctGTTCCC‑TTCCGCCAGCTGC‑TGGGTGGGC‑GGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTCGCGCGAGGAAGT‑CGCGAGAAACCCGCCGTTGCTACAACAGGGGGTGCGACCGCATCCGGTGAACGGTAAACAGGCGCTATTTGTGAATGAAGGCTTTACTCCGCGCATTGTTGATGTGAGCGGGAAAGAGAGCGCAGCCttttt                                                                                                                             >  2:745197/1‑247 (MQ=255)
                                                                                                                             gtcccc‑TTCCGCCAGTTGC‑TGAGTGGGC‑TCGGTGCGGAGCATTATTTCCGTAAATCCTTCCCGGAATA‑CAAATCCCGAAAAACCGGGGAGGAACATCAGCCCGGGGGCGGGGCGGC‑CGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAg                                                                                                                          <  1:442972/248‑1 (MQ=255)
                                                                                                                               tcccattCCGCCAGCTGT‑TGAGTGGGC‑TGCGGGCGGAGCATGATTTCCGTAAACCGTTCCCGGAAAA‑AAAATCCCGAAAAGCCGAGGAGGAACACCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAg                                                                                                                          <  1:204432/245‑1 (MQ=255)
                                                                                                                                        gCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTtttgttt                                                                                                                >  1:782375/1‑251 (MQ=255)
                                                                                                                                        gCCAGCTGC‑TGAGGGGGC‑TGCGTGCGGAGCATGATTTCCGTAAACCGTTCCCGGAATA‑CAAATACCGAAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTtttgttt                                                                                                                <  1:204745/251‑1 (MQ=255)
                                                                                                                                               gc‑tgAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCat                                                                                                          >  1:295861/1‑250 (MQ=255)
                                                                                                                                                 cggAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCACCGCTGGCGCGGGGCGGT‑CGCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCatat                                                                                                        <  1:121723/249‑1 (MQ=255)
                                                                                                                                                  tgAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAGTA‑CAACTACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGGGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATc                                                                                                       >  2:213583/1‑251 (MQ=255)
                                                                                                                                                  tgAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATc                                                                                                       >  1:572062/1‑251 (MQ=255)
                                                                                                                                                            cgcggggtagacATGATTTCCGTAAAGC‑TTCCCGGAAAAACAAACCCGGCAAAACCGGGGAGGAACATCAAGGCGGGGGCGGGGCGGT‑GGGGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACGGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCg                                                                                             <  2:87287/240‑1 (MQ=255)
                                                                                                                                                                 gcgGGGCATGATTTCGGAAATTCGTTCCCGGAATA‑CAAATCCCGCAAAACCGAGGAGGAACACCAACGCTGGCGCGAGGCGTT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGtt                                                                                        <  1:845246/251‑1 (MQ=255)
                                                                                                                                                                    gAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGCGTTTc                                                                                      >  1:505298/1‑250 (MQ=255)
                                                                                                                                                                       cATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGCGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCCTATCA‑CCAAACCGGAGTTTCAggg                                                                                  >  2:384071/1‑250 (MQ=255)
                                                                                                                                                                       cATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGGTTCAGGt                                                                                  >  2:771409/1‑251 (MQ=255)
                                                                                                                                                                          gATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTgcg                                                                               >  1:130585/1‑251 (MQ=255)
                                                                                                                                                                          gATTTCCGTAAATCGTTCCCGAACTA‑CAAATACCGCAAAACCGAGGAGGAACATAAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATACGGTGGTGCGCACGCCTCCGGTGAGCGGTACACAGGCGCTGTTTGTGAATGAAGGATTTCCTACGCGAATTGTTGATGTGCGCGAGAAAGAGAGCGACGCCTTGTTAAGTTTTTTGTTTGCCCATATAA‑CCCACCCGGAGTTTCAGGTgcg                                                                               >  2:809555/1‑251 (MQ=255)
                                                                                                                                                                             ttCCGTAAATCGTTCCCGAAATA‑CAAATACCGAAAAACCGAGGAGGAACATCATCGCGGGCGCGGGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTgcgct                                                                             <  1:424401/250‑1 (MQ=255)
                                                                                                                                                                              tCCGAAAATCTTTCCCGGAATA‑AAAATCCCGCAAACCCGGGGAGGAACATCAGCGCTGGCGCGAGGCGGC‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTgcgctg                                                                            <  1:653957/250‑1 (MQ=255)
                                                                                                                                                                                cgaaaagcttgcctgaaaTA‑AAAATACCGAAAAACCGAGGATGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTTGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTgcgctggcg                                                                         <  1:521066/235‑1 (MQ=255)
                                                                                                                                                                                    aaTCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCTCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGcgctggcgctgg                                                                     >  1:268206/1‑251 (MQ=255)
                                                                                                                                                                                     aTCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACACGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGCGTTTCAGGTgcgctggcgctagc                                                                    >  2:258697/1‑248 (MQ=255)
                                                                                                                                                                                              ggAATA‑CAAATACCGCAAAACCGAGGAGGCACCTCAACGCTGGCGCGAGGCGGA‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGCACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGACTTGTTGATGTGAGAGAGAAAGAGAGCGAAGCCTTCTTAAGTTTTTTGTTTGCTCAGATCA‑CCAAACCGGAGTTTCAAGTGCGCTGGCGCTGGCCACCACATg                                                           >  2:38508/1‑251 (MQ=255)
                                                                                                                                                                                                aaTA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGat                                                         >  1:412812/1‑251 (MQ=255)
                                                                                                                                                                                                          aCCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAACAACCCGCCGTTGCTACATCCGGTGGTGCGGACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGATAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCAATATTA‑CCAAACCGGAGTTTCAGGAGCGCAGGCGCTGGCCACCAAATGATATTGCGAtt                                                >  1:308147/1‑251 (MQ=255)
                                                                                                                                                                                                           ccGCAAAACCGAGGAGGAACATCAACGCTAGCGCGAGTCGGT‑CGCGAAAAACACGCCGGTGCTACATCCGGTGGTTCGCACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTAATACGCGAATTGTGGATGTGAGGGCGAAAGAAAGCGTAGCCTTGTTAAGTTTTTTGTTGGCCCATA‑CAGCCACAGCGTAGTTTCAGGGGCgctggcgcgggccacgacactgtagtgcgattt                                               >  2:820482/1‑230 (MQ=255)
                                                                                                                                                                                                                aacccgggacgaaaCACCAACGCGGGCGCGGTGGGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGAt                                          <  2:694258/239‑1 (MQ=255)
                                                                                                                                                                                                                aaaCCGGGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTc                                                                                      <  1:482863/207‑1 (MQ=255)
                                                                                                                                                                                                                aaaCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGACGTTGCTACATCCGGTGGTGCGAACGCATCCGGAGAGCGGTAAACAGGCGCTGTTTGTGAATGAATGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGGTACGTTTTTTGTTTGCCCATATCA‑CCAAACCGGGGATTCAGGTGCGCTGGCGCTGGCAACCAAACGAGATTGCGAGTTGGGAt                                          >  2:749940/1‑251 (MQ=255)
                                                                                                                                                                                                                      aggaggAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCaa                                                              >  1:378558/1‑225 (MQ=255)
                                                                                                                                                                                                                            aaCATCCACGCTGGCGCGAAGAGGT‑CGCTCAAAACCCGCAGTTGCTACATCAGGCGGTGCGAACGCATCCGGTGAGCGGACAACCGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGGAAGAGAGCGAAGCCTTGTTACGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGATTCAGGCGCGGTAGCGGTGACACCCTAATGATAATGCggcatgtggacaacgcgagaac                              >  2:555525/1‑230 (MQ=255)
                                                                                                                                                                                                                               accAACGCTGGCGCGAG‑CGGTTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAg                           <  1:96205/249‑1 (MQ=255)
                                                                                                                                                                                                                                tCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑ACAAACCGGAGTTTCAGGTGGGCTGGCGCTGGCCACCAAATGATATTGAGATTTGGGGTAACCGGGTGAcccggc                          >  2:448422/1‑248 (MQ=255)
                                                                                                                                                                                                                                  aactcgggCGCGAGGCGGC‑CCAGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTGTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCAc                        <  1:150787/245‑1 (MQ=255)
                                                                                                                                                                                                                                              gggCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGcc            <  1:384071/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gggCGGT‑CGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGcc            <  1:771409/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGGT‑CGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCg           >  1:709311/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCGGT‑CGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCg           <  2:180481/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gcaaaaacCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCt    <  1:26758/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGc  >  1:222933/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         aaaaacCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGGTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGc  <  2:266449/250‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
GCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTA‑TTGCGGCCTATGAGGCGCTCTCTGTTCCC‑TTCCGCCAGCTGC‑TGAGTGGGC‑TGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGT‑CGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCA‑CCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGC  >  NC_000913/388088‑388581

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