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.

GCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGC  >  NC_000913/388088‑388581
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
gCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAAcc                                                                                                                                                                                                                                                      <  1:697914/251‑1 (MQ=255)
                     tgacccGCCACCCTGAGGGGCGATTCTGGAAGCAAAAGAGTTCCTTTCGACGGGGGGTGATACGCTCTGGACCAGCGGTATTGGGGCCTATGAGGCGCTCTCTGTTCCCTCCCGCCAGCTGCTGAGTGGGCTGCTGGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAATTCCCGCAAACCCGAGGAGAAACATCACCGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCggtg                                                                                                                                                                                                                                  <  1:852578/245‑1 (MQ=255)
                                 ggggggggCGATGCGGGCAGATAAAGATTTACTTTCGCCGGGCGGTTATCCGCTCTGGACCAGGGGTATTGGGGCTAAGGGGGCGTCTTCGTTCCCCTCCCGCCGGCTGCTGAGTGGGCTGCGGGCGGAGCTGGATTTCCGAAAATCGTTCCCGGAATACAAATCCCTCAAAACCGAGGAGAAACACCAACGCGGGCGGGGGGCGTTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATc                                                                                                                                                                                                                     <  2:772673/247‑1 (MQ=255)
                                               gggCAGCTAAAGAGTTACCTCCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCGGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTaaa                                                                                                                                                                                                       <  1:319173/250‑1 (MQ=255)
                                                  cagttAAAGAGTTCCCTTCGACCGGCGGTGATCCGCTCGGGACCAGCGGTATTGCGGCCTAGGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGTTGAGTGGGTTGCGTGCGGAGAATGATTTCCGTAAATCGTTCCCGAAATAAAAATACCGCAAAACCGAGGAGGAACATCACCGCGGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAg                                                                                                                                                                                                    <  1:62656/247‑1 (MQ=255)
                                                           agTTACCTTCGACCGGCGGTGATACGCTCTGGACCCGCGGTATTGCGGCCTCTGGGGCGCTCTCTGTTCCCTTCCGCCCGCTGCTGGGTGGGCTGCGTGCGGGGCATGATTTCCGTAAATCGTTCCCCGAATACAAATACCGCAAAACCGAGGAGGGACATCCACGCTGGCGCGGGGCGGGCGCGGAAAACCCGCCGGTGCTACATCCGGGGGTGCGAACGCACCCCGTGGGCGGGAAAAAggaccggttt                                                                                                                                                                                           >  1:573997/1‑243 (MQ=255)
                                                            tttacctgcgacaggcgtggaagcGCTCTGGACCAGCGTAGTTGCGGCCTATGAGGCGCTCTGTGTTCCCTTCCGCCTATGCCTAAGTGGGCTGCTTGCTGAGCATGTTATCCGTAAATCGTTCCCGGAATACAAATCCCGCAAAACCGAGGAGGAACATAATCTCTGGCGCGGGGCGGTCGCGAAAAGCCCGCCGTTGCTTCAGCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGttt                                                                                                                                                                                           <  1:251496/227‑1 (MQ=255)
                                                                     gACCGGCGGTGATACGTTCTGGACCAGCGGTTTTGCGGCCTATGATGCGTTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCTTGCGGAGCATAATTTCCGAAAATCTTTCCCGGAATACAAATACCGCAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAg                                                                                                                                                                                 <  1:36340/251‑1 (MQ=255)
                                                                        cggcggTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCGCTGTTCCCTCCCGCCAGCTGCTGGGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATAAAAATACCGCAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCt                                                                                                                                                                              <  1:803366/251‑1 (MQ=255)
                                                                            ggtgatcctctctgaccctgctgTATTGCGTCTTATGGGGCTCTCTCTGTCCCCTTCCGCCAGCTGCTGAGGGGGCGGCTTGCGGAGCATCATTTCCGAAAATCGTCCCCGTAAAACAAATACCGCAAAACTGAGGAGAAACATCACCCCGGTCGCGGGGCGTTCGCGAAAAACCCGCCGTTGCTACATCCGGGGGTGCGAACGCATCCGGTGAGCGGTAAACAGTCGCTGTTTGTGAATGAAGGCTTtac                                                                                                                                                                          <  1:266618/229‑1 (MQ=255)
                                                                                atccGCTCGGGACCGGCGGTAGGGCGGCTTAGGAGGCGCTCTCTGTCCCCTCCCGCCATCGGCTGAGGGGCGTGCGGGCGGAGCATGATTCCCGAAAATCGTTCCCGAAATAAAAATACCGCAAAACCGAGGGGGAACATCAACGCGGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcg                                                                                                                                                                      <  1:45459/248‑1 (MQ=255)
                                                                                 tACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcgc                                                                                                                                                                     >  1:180481/1‑251 (MQ=255)
                                                                                          taCCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAAttgttg                                                                                                                                                            >  1:266449/2‑251 (MQ=255)
                                                                                                             cTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGTAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAagagg                                                                                                                                          >  1:101345/1‑250 (MQ=255)
                                                                                                                  gggCGCTCTCTGTCCCCTCCCGCCAGCTGCTGAGGGGGCTGCGTGCGGAGCATGATTCCCGAAAATCGTCCCCGGAATACAAATCCCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGaa                                                                                                                                     <  1:191226/250‑1 (MQ=255)
                                                                                                                     cgcTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAAACCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTCCGCGAATTGGTGATGTGAGCGGGCAAGGGAGCGAAGcc                                                                                                                                  >  2:253447/1‑251 (MQ=255)
                                                                                                                           ctGTTCCCTTCCGCCAGCTGCTGGGTGGGCGGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTCGCGCGAGGAAGTCGCGAGAAACCCGCCGTTGCTACAACAGGGGGTGCGACCGCATCCGGTGAACGGTAAACAGGCGCTATTTGTGAATGAAGGCTTTACTCCGCGCATTGTTGATGTGAGCGGGAAAGAGAGCGCAGCCttttt                                                                                                                             >  2:745197/1‑247 (MQ=255)
                                                                                                                             gtccccTTCCGCCAGTTGCTGAGTGGGCTCGGTGCGGAGCATTATTTCCGTAAATCCTTCCCGGAATACAAATCCCGAAAAACCGGGGAGGAACATCAGCCCGGGGGCGGGGCGGCCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAg                                                                                                                          <  1:442972/248‑1 (MQ=255)
                                                                                                                               tcccattCCGCCAGCTGTTGAGTGGGCTGCGGGCGGAGCATGATTTCCGTAAACCGTTCCCGGAAAAAAAATCCCGAAAAGCCGAGGAGGAACACCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAg                                                                                                                          <  1:204432/245‑1 (MQ=255)
                                                                                                                                        gCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTtttgttt                                                                                                                >  1:782375/1‑251 (MQ=255)
                                                                                                                                        gCCAGCTGCTGAGGGGGCTGCGTGCGGAGCATGATTTCCGTAAACCGTTCCCGGAATACAAATACCGAAAACCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTtttgttt                                                                                                                <  1:204745/251‑1 (MQ=255)
                                                                                                                                               gctgAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCat                                                                                                          >  1:295861/1‑250 (MQ=255)
                                                                                                                                                 cggAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCACCGCTGGCGCGGGGCGGTCGCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCatat                                                                                                        <  1:121723/249‑1 (MQ=255)
                                                                                                                                                  tgAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAGTACAACTACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGGGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATc                                                                                                       >  2:213583/1‑251 (MQ=255)
                                                                                                                                                  tgAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATc                                                                                                       >  1:572062/1‑251 (MQ=255)
                                                                                                                                                            cgcggggtagacATGATTTCCGTAAAGCTTCCCGGAAAAACAAACCCGGCAAAACCGGGGAGGAACATCAAGGCGGGGGCGGGGCGGTGGGGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACGGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCg                                                                                             <  2:87287/240‑1 (MQ=255)
                                                                                                                                                                 gcgGGGCATGATTTCGGAAATTCGTTCCCGGAATACAAATCCCGCAAAACCGAGGAGGAACACCAACGCTGGCGCGAGGCGTTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGtt                                                                                        <  1:845246/251‑1 (MQ=255)
                                                                                                                                                                    gAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGCGTTTc                                                                                      >  1:505298/1‑250 (MQ=255)
                                                                                                                                                                       cATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGCGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCCTATCACCAAACCGGAGTTTCAggg                                                                                  >  2:384071/1‑250 (MQ=255)
                                                                                                                                                                       cATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGGTTCAGGt                                                                                  >  2:771409/1‑251 (MQ=255)
                                                                                                                                                                          gATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcg                                                                               >  1:130585/1‑251 (MQ=255)
                                                                                                                                                                          gATTTCCGTAAATCGTTCCCGAACTACAAATACCGCAAAACCGAGGAGGAACATAAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATACGGTGGTGCGCACGCCTCCGGTGAGCGGTACACAGGCGCTGTTTGTGAATGAAGGATTTCCTACGCGAATTGTTGATGTGCGCGAGAAAGAGAGCGACGCCTTGTTAAGTTTTTTGTTTGCCCATATAACCCACCCGGAGTTTCAGGTgcg                                                                               >  2:809555/1‑251 (MQ=255)
                                                                                                                                                                             ttCCGTAAATCGTTCCCGAAATACAAATACCGAAAAACCGAGGAGGAACATCATCGCGGGCGCGGGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgct                                                                             <  1:424401/250‑1 (MQ=255)
                                                                                                                                                                              tCCGAAAATCTTTCCCGGAATAAAAATCCCGCAAACCCGGGGAGGAACATCAGCGCTGGCGCGAGGCGGCCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctg                                                                            <  1:653957/250‑1 (MQ=255)
                                                                                                                                                                                cgaaaagcttgcctgaaaTAAAAATACCGAAAAACCGAGGATGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTTGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctggcg                                                                         <  1:521066/235‑1 (MQ=255)
                                                                                                                                                                                    aaTCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCTCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGcgctggcgctgg                                                                     >  1:268206/1‑251 (MQ=255)
                                                                                                                                                                                     aTCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACACGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGCGTTTCAGGTgcgctggcgctagc                                                                    >  2:258697/1‑248 (MQ=255)
                                                                                                                                                                                              ggAATACAAATACCGCAAAACCGAGGAGGCACCTCAACGCTGGCGCGAGGCGGACGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGCACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGACTTGTTGATGTGAGAGAGAAAGAGAGCGAAGCCTTCTTAAGTTTTTTGTTTGCTCAGATCACCAAACCGGAGTTTCAAGTGCGCTGGCGCTGGCCACCACATg                                                           >  2:38508/1‑251 (MQ=255)
                                                                                                                                                                                                aaTACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGat                                                         >  1:412812/1‑251 (MQ=255)
                                                                                                                                                                                                          aCCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAACAACCCGCCGTTGCTACATCCGGTGGTGCGGACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGATAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCAATATTACCAAACCGGAGTTTCAGGAGCGCAGGCGCTGGCCACCAAATGATATTGCGAtt                                                >  1:308147/1‑251 (MQ=255)
                                                                                                                                                                                                           ccGCAAAACCGAGGAGGAACATCAACGCTAGCGCGAGTCGGTCGCGAAAAACACGCCGGTGCTACATCCGGTGGTTCGCACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTAATACGCGAATTGTGGATGTGAGGGCGAAAGAAAGCGTAGCCTTGTTAAGTTTTTTGTTGGCCCATACAGCCACAGCGTAGTTTCAGGGGCgctggcgcgggccacgacactgtagtgcgattt                                               >  2:820482/1‑230 (MQ=255)
                                                                                                                                                                                                                aacccgggacgaaaCACCAACGCGGGCGCGGTGGGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGAt                                          <  2:694258/239‑1 (MQ=255)
                                                                                                                                                                                                                aaaCCGGGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTc                                                                                      <  1:482863/207‑1 (MQ=255)
                                                                                                                                                                                                                aaaCCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGACGTTGCTACATCCGGTGGTGCGAACGCATCCGGAGAGCGGTAAACAGGCGCTGTTTGTGAATGAATGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGGTACGTTTTTTGTTTGCCCATATCACCAAACCGGGGATTCAGGTGCGCTGGCGCTGGCAACCAAACGAGATTGCGAGTTGGGAt                                          >  2:749940/1‑251 (MQ=255)
                                                                                                                                                                                                                      aggaggAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCaa                                                              >  1:378558/1‑225 (MQ=255)
                                                                                                                                                                                                                            aaCATCCACGCTGGCGCGAAGAGGTCGCTCAAAACCCGCAGTTGCTACATCAGGCGGTGCGAACGCATCCGGTGAGCGGACAACCGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGGAAGAGAGCGAAGCCTTGTTACGTTTTTTGTTTGCCCATATCACCAAACCGGAGATTCAGGCGCGGTAGCGGTGACACCCTAATGATAATGCggcatgtggacaacgcgagaac                              >  2:555525/1‑230 (MQ=255)
                                                                                                                                                                                                                               accAACGCTGGCGCGAGCGGTTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAg                           <  1:96205/249‑1 (MQ=255)
                                                                                                                                                                                                                                tCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCAACAAACCGGAGTTTCAGGTGGGCTGGCGCTGGCCACCAAATGATATTGAGATTTGGGGTAACCGGGTGAcccggc                          >  2:448422/1‑248 (MQ=255)
                                                                                                                                                                                                                                  aactcgggCGCGAGGCGGCCCAGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTGTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCAc                        <  1:150787/245‑1 (MQ=255)
                                                                                                                                                                                                                                              gggCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGcc            <  1:384071/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gggCGGTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGcc            <  1:771409/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCg           >  1:709311/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCGGTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCg           <  2:180481/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gcaaaaacCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCt    <  1:26758/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGc  >  1:222933/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         aaaaacCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGGTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGc  <  2:266449/250‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
GCATACCGATGTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGC  >  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: