Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,228,227	G→C	17.8%	pseudogene (381/1287 nt)	mdtQ ←	pseudogene, putative channel/filament proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,228,227	0	G	C	17.8%	111.8 / 4.1	55	pseudogene (381/1287 nt)	mdtQ	pseudogene, putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (17/26);  new base C (8/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCA  >  NC_000913/2227994‑2228452                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               ttcattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGTTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCCGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGA‑GGG‑GTTTCTcgccgtc                                                                                                                                                                                                                                           >  2:1366354/1‑235 (MQ=255)   cattttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTccg                                                                                                                                                                                                                                         <  1:1268961/239‑1 (MQ=255)    attttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTccgc                                                                                                                                                                                                                                        <  2:490828/239‑1 (MQ=255)        tcCCCTCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGAGATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                    <  1:740733/239‑1 (MQ=255)        tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                    <  1:1050857/239‑1 (MQ=255)        tcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGCGCTACCCTGCCAGCCCGCAATTGGC‑GGG‑TTTTGTCTCGTTTccccgc                                                                                                                                                                                                                                     >  2:524944/1‑238 (MQ=255)        gccccGCGACATCGGTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCAAACCTTAAAATGAGGAAGTAGTCCCGTTCTTATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCGATCTGCTGCAAGACGGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCgccgcc                                                                                                                                                                                                                                    <  2:731897/238‑1 (MQ=255)               gACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGT‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑tt                                                                                                                                                                                                                            <  2:651489/239‑1 (MQ=255)                  aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTaa                                                                                                                                                                                                                         <  2:1040693/239‑1 (MQ=255)                  aTCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTaa                                                                                                                                                                                                                         <  1:718385/239‑1 (MQ=255)                       tGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTa                                                                                                                                                                                                                    <  2:1272767/239‑1 (MQ=255)                           ctgctgCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTAccc                                                                                                                                                                                                                 <  2:699243/238‑1 (MQ=255)                             gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGCCCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAgg                                                                                                                                                                                                              <  2:275186/239‑1 (MQ=255)                             gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAgg                                                                                                                                                                                                              <  1:744827/239‑1 (MQ=255)                              ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTGTTCGCTTTCCCCCCCCCCGCC‑TTTACCCGTcccccgc                                                                                                                                                                                                              >  2:849964/1‑236 (MQ=255)                              ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCCTAccccagc                                                                                                                                                                                                             >  1:658867/1‑235 (MQ=255)                               tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGc                                                                                                                                                                                                             <  1:1394819/238‑1 (MQ=255)                               tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTCTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGc                                                                                                                                                                                                             <  1:754416/238‑1 (MQ=255)                                gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTTCGCTGCCAGCCAGCAATTGGC‑GGG‑TTTTCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCgg                                                                                                                                                                                                           <  1:459330/239‑1 (MQ=255)                                  cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCggg                                                                                                                                                                                                          <  1:812763/238‑1 (MQ=255)                                                          gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCCCCCGTGC‑TTTAACCGTACCCAAGCGGGCAGTAA‑ACTCCGCCCGGTT‑Ctttt                                                                                                                                                                                >  1:33122/1‑237 (MQ=255)                                                            ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCGATT‑CTTTccc                                                                                                                                                                              <  1:293422/238‑1 (MQ=255)                                                             tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTCGCGGTCCGCCGCCCGTGC‑TTTTACCGTACCCAGGCGGGCAGTAA‑CCCCCCCCCGGTT‑CTTTTCCCa                                                                                                                                                                            >  1:454244/1‑239 (MQ=255)                                                                     ttCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCCCCCGGGC‑TTTAACCGTACCCCGGCGGGGCGTAA‑CCTCCCCCCGATT‑CTTTcccccagttttca                                                                                                                                                                    >  1:1065543/1‑230 (MQ=255)                                                                         aCTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTTTCGTTGCCCCCCCCCCGC‑TTTAAACCGGCCCAAGTCGGCAGACA‑CCCCCCCCCCATTACTTTcccagcagacgaaat                                                                                                                                                                  >  2:416005/1‑226 (MQ=255)                                                                              agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGGAA‑CCTCCGCCCGATT‑CTTTCCCCCG‑ATATCGAG‑ATG‑CCAg                                                                                                                                                         >  1:1159355/1‑238 (MQ=255)                                                                              agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGGTCCCCCCCCCGTGC‑TTTTACCCGACCCCGGCGGGCCGTAA‑CCCCCCCCCCGTT‑CTTTTCCCCC‑GTATCCCG‑ATG‑cccg                                                                                                                                                         >  2:300563/1‑236 (MQ=255)                                                                              agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTGTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCCGTAA‑CCTCCGCCCCATT‑CTTTCCCCCG‑ATATCgcgcagccaa                                                                                                                                                           >  1:1060747/1‑229 (MQ=255)                                                                                           ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTGCTCGCCTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGGAA‑ACTCCGCCCGATT‑CTTTTCCCAG‑ATATCGCG‑ATG‑CCAGCCCCCC‑GCTAAG‑cc                                                                                                                                         >  1:892628/1‑239 (MQ=255)                                                                                                atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTCGCGGTTCCCCGCCCGTGC‑CTTTACCCTACCCCAGCGGGCCATAA‑CCTCCCCCCCGTT‑CTTTCCCCCA‑AGATCCCG‑AAG‑CCCACCCCCC‑GTTTAG‑ccccaaa                                                                                                                                    >  2:1247715/1‑234 (MQ=255)                                                                                                       gggggaccggcgcggATATTGTTGTCCTGCTCTTTTCCTATTTGCGGCAGCACCGTTTTGACCGCCGCC‑GGGTTTGCCCCCCCCCGGTACGGGCGGCCTCCGCGGCCAGCCAGCAATGGGC‑GGG‑TTTGCGCGCGTCCCCCCCCCCGTGC‑TTTAACCGTCCCCGGGCGGGCAGAAC‑CCCCCCCCCGATT‑CTTCCCCCAG‑ATAGCGAG‑ATC‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGt                                                                                                                              <  2:431115/224‑1 (MQ=255)                                                                                                          cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTTTTCGCCGTCCGCCACCCCGGC‑TTTAACCCGACCCCGGCGGGGAGTGA‑ACTCCGCCCGGTT‑CTTTTCTCAG‑ATAATGAG‑GCG‑CGACCCCGCC‑CCTTAG‑CCCCAAGAACCGCTggg                                                                                                                          >  2:76914/1‑238 (MQ=255)                                                                                                          cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTGCTCGCGTTCCGCCGCCCCGGC‑TTTTACCCTAACCAGGCGGGCCGTAA‑CCTCCGCCCGGTT‑CTTTTCCCCG‑ATATCGAG‑ATG‑CCCGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTTgt                                                                                                                          >  1:509571/1‑239 (MQ=255)                                                                                                                ggcgcaatatgggTTTTCTGCTTTTTTTCTTTTTGCTGAAGGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACGGGGGGGCTACGCTGCCACCCAGCAATGGGC‑GGT‑TTTGCTCGCGTCCCGCCGCCCGTGC‑TTTACCCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTAcc                                                                                                                     <  2:1159355/227‑1 (MQ=255)                                                                                                                               ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTTCTCCCCGTTTGCCGCCCCCGC‑CTTTACCCTACCCCCGCGGGCCGTAA‑ACCCCCCCCCCCT‑CTTTTCCCAG‑CTATCCAC‑ATGCCCAGCCCCCC‑GTAACC‑CCAAAAAACCCTTTGGT‑TGCCCCGGGCCGGGC‑Gtca                                                                                                     >  2:10902/1‑236 (MQ=255)                                                                                                                                         tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTCGAGGTCCGCCCCCCCGGCTTTTAACCTTACCCAGCGGGGGATTAA‑CC‑CCGCCCCCTT‑CTTTCCCCCA‑ATATCGAG‑ATG‑CCACCCCCCC‑GTTAAG‑CCCAAAATACCCTTGGTGTACCACGGGCCCGGC‑GGACCTGA‑GGcccgg                                                                                         >  2:704923/1‑236 (MQ=255)                                                                                                                                         tttCGATTTGCGGCAAGACCGTGTTTACCGCCGCTGGGGTTTG‑CCCCTCCCAGTACGGGCGGGCTTCGCGGCCAGCCGGAAATGGGC‑GGG‑TTTGCTCGCTTCCCGCCGCCCGTGC‑TTTAACCGTACCCGGGCGGGCAGTAC‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGa                                                                                        <  2:1065543/239‑1 (MQ=255)                                                                                                                                                           ccGTTTTTACCCCCCCCTGGGTTTC‑CACCCCCCTGTGCGGGCGGCCTCCCCCCCCAGCCAGCAATTGGC‑GGG‑TTTGCCCCCTTTCCGCCCCCCGTCC‑TTAACCCTTCCCCGGGCGGGCAGTAC‑CCCCCGCCCGATT‑CTTTCCCCAG‑ATACCGAG‑ATG‑CCGGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGc                                                                     <  2:33122/239‑1 (MQ=255)                                                                                                                                                              tgtTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGG‑TTTCCCCCCCTTCCCCGCCCCCGGT‑TTTAACCG‑ACCCCGGCGGGGCGTAACCCCCCGCCCGTTT‑TTTTCCCCAG‑‑TAACCAGAATG‑CCACCCCCCC‑GTTAAC‑CCAAAAGACCCGTGGGTGTACC‑CGGGCCGGTC‑GGCACCGC‑GCCC‑CGAACGTTCCGAGC‑ACACGGccc                                                                   >  2:679777/1‑238 (MQ=255)                                                                                                                                                              tgtTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCCACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGGTCCCCCGCCCGTGG‑TTTAAACGTACCCAAGCGGGCAGGTA‑ACTCCCCCCGATT‑CTTTTCCCCG‑ATATCGAG‑GAG‑CCAGCCCGCC‑CGTAAG‑CCCAAAATGCCGCTGGGGTACCAACGGCCGGTC‑GTACCTGC‑GGGC‑GGGTCGTTCCGGGCGAAACGGcc                                                                    >  2:416611/1‑238 (MQ=255)                                                                                                                                                                     cgccgcCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑CTTAACCCGACCCAGGCGGGCCGTAA‑ACCCCCCCCGGTT‑CTTTTCCCCG‑ATATCCCG‑ATG‑CCCACCCCGC‑GTTAAG‑CCCAAAAGACCGTTGGTGTACCCCGGGCCCGGC‑GGCCCCGC‑GGCC‑GGgttgtttcgcagaaaagggcccacttaa                                                             >  2:515561/1‑216 (MQ=255)                                                                                                                                                                         gccTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCCCCCGGTCCCCCCCCCCTGC‑CTTAACCCCACCCAGGCCGGCAGCAA‑CCCCCGCCCCCGT‑TTTTCCCCAC‑ATATCGCG‑TTG‑CCCACCCCCC‑GCTAAG‑CCCCAAACACCCCTGGTGTTCCACCGGCCCGTC‑CGCCCTGC‑GGCC‑CGACAGTTCTGAGC‑aaaaagcgcccctaagaact                                                        >  2:779470/1‑221 (MQ=255)                                                                                                                                                                          ccTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTGCTCCCCGTCCTCCCCCCCGGC‑TTTAACCCCACCCCGGCGGGCCGTAA‑CATCCGCCCCCTT‑CTTTCCCCAC‑GTATCCAG‑ATG‑CCAGCCCGCC‑CGTAAA‑CCCAAAATACCCCTTGTGGACCACGGGGCCGGC‑GGACCCGC‑GGCCGGGATCGTCCA‑AGC‑AACCGCCCCCCTTACGCCTTc                                                       >  2:410459/1‑238 (MQ=255)                                                                                                                                                                                         tCCCATTACAGGCGGGCTCCGCTGCCAGCAAGCAATGGGC‑GGG‑TTTGCCGGCGTTCCGCCGCCCGTGC‑TTTACCCGTCCCCAGGCGGGCAGTAA‑CCCCCGCCCTATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTg                                       <  1:1015008/239‑1 (MQ=255)                                                                                                                                                                                          cccGGTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGC‑GGG‑TTGGCTCGCGTCCCGCCGCCCGTGC‑TTTAACCGTCCCCAGGGGGGCAGTAC‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGc                                      <  2:509571/239‑1 (MQ=255)                                                                                                                                                                                             aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCg                                                                                                                                                                                         <  2:212451/100‑1 (MQ=255)                                                                                                                                                                                                aCAGGCGGGCTACGCTGCCAGCCACCAATTGGC‑GGG‑TTTCCTCGCGTCCCGCGGCCGGTGC‑TTTAACCGTCCCCGGGCGGGCAGTAC‑CCCCCGCCCGATT‑CTTTCCCCAG‑ATACCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGc                                   <  1:1247715/236‑1 (MQ=255)                                                                                                                                                                                                    gCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCCCCCGTTCCCCCCCCCGCGC‑TTTTAACCTACCCCGGCGGGCAGGAA‑CCTCCCCCCCATT‑TTTTTCCCCG‑ATATCGGG‑AAG‑GCCCCCCGCC‑GTTTAC‑CCCAAAAAACCCGTGGTGCACCACCGGGCGGTCGGTACCCTC‑GTCC‑GGGAGCTACGGCAC‑AAACGGCCCCCTTAAGCCCTCTGTCG‑CAAttttttacacgcccttt                             >  2:1000316/1‑228 (MQ=255)                                                                                                                                                                                                    gCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTGTCCCGGTCCGCCGCCCGTGC‑TTTAACCGCAACCAGGCGGGGCGGAA‑CCTCCGCCCCGTT‑CTTTCCCCCG‑ATATCGAG‑ATG‑CCACCCCCCC‑GTTTAG‑GCCAAAATACCCTTGGGGTACCACGGGCCGGGC‑GTGCCCGG‑GGCC‑GCGGCGTTCAGAGC‑AAAACGCCCCACTAAGCCCTCCTCCG‑GCCTTCTTTTCCggccccac                             >  2:831934/1‑231 (MQ=255)                                                                                                                                                                                                    gCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTCGCGGTCCCCCGCCCCTGC‑CTTAACCGTACCCACGCGGGCAGTAA‑ACCCCGCCCCATT‑CTTTCCCCCG‑CTATCGAG‑AGG‑CCAGCCCCCC‑GGTAAGACCAAAAG‑AACGTTGGTGTACCACGGGGCCGCC‑GGAACCGC‑GGCC‑CGATCGTTCCTAGC‑AAACCGCCCCATTAACCCTCCTGCCGCCCATTTTTTTccctcctac                              >  2:22901/1‑231 (MQ=255)                                                                                                                                                                                                       ggCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTTGTCCCGTTTCGCCCCCCGCGC‑CTTAACCGTACCCCGGCGGGGAGTAA‑CCTCCCCCCCATT‑CTTTccc                                                                                                                                                                              >  1:212451/1‑100 (MQ=255)                                                                                                                                                                                                       ggCTACGCTGCCAGCCAGCAATTGGC‑GGG‑GTTTCTCCCGTTCCCCCGCCCCGGC‑TTTTAACCTACCCCGGCGGGGAGGAA‑CCTCCCCCCCACT‑TTTTTCCCCC‑ATAACGAG‑AAG‑CCCGCCCCCCGGTTAAC‑CCCAAACAACCGTT‑GTGTACCACCGGCCCGTC‑GGACCTGCGGGCC‑GGACGGTCCAGGCC‑AAACGCCCCCATTACCCTTTCTCCGG‑CCATTTTTTGCCGcccaaactc                           >  2:334400/1‑230 (MQ=255)                                                                                                                                                                                                                      cagcAATTGGC‑GGG‑TTTGCGCGCGTCCCGCCCCCCGTGT‑TTTAACCGTACCCGGGCGGGCATTAC‑CCCCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCt          <  1:853612/239‑1 (MQ=255)                                                                                                                                                                                                                           atgggC‑GGG‑TTTGCTCGCGTCCCGCCGCCCGTGC‑TTTACCCGTCCCCAGGCGGGCAGTAC‑CCCCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa     <  1:300563/236‑1 (MQ=255)                                                                                                                                                                                                                           aggggC‑GTG‑TTTGCCCCCTTCCCGCCCCCCGTAC‑TTTACCCGTCCCCAGGCGGGCACTAA‑CCTCCCCCCTATT‑C‑TTCCCCAGTATACCGAG‑CGG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCGTATCCGCCGAAAAGTCGATCTGCGGa     <  2:645117/236‑1 (MQ=255)                                                                                                                                                                                                                               gC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGCTCTGCGGACCa  >  1:1086643/1‑238 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               TTCATTTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTG‑CCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGG‑TTTGCTCGCGTTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAGTAA‑CCTCCGCCCGATT‑CTTTCCCCAG‑ATATCGAG‑ATG‑CCAGCCCGCC‑GTTAAG‑CCAAAAGTACCGTTGGTGTACCACGGGCCGGTC‑GTACCTGC‑GGCC‑GGATCGTTCAGAGC‑AAACGGCCCCATTAAGCCTTCTGCCG‑ACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCA  >  NC_000913/2227994‑2228452

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