Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,228,272	A→C	27.3%	pseudogene (338/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,272	0	A	C	27.3%	102.0 / 7.5	63	pseudogene (338/1287 nt)	mdtQ	pseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base A (37/6);  new base C (3/17);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 7.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGT  >  NC_000913/2228041‑2228505                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑T‑CTC‑GC‑GTTCCCCCC‑CCCGGGGC‑TTAACCCGTC‑CCCAGCC‑GGG‑AGT‑AACccccc                                                                                                                                                                                                                                               >  1:248805/1‑237 (MQ=255)      aaTATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑T‑CTC‑CC‑GGTCCCCCC‑CCC‑CGGC‑CTTAACCCTA‑CCCCGGG‑GGGGAGT‑AAACTCCTCCCC‑At                                                                                                                                                                                                                                        >  1:544568/1‑239 (MQ=255)           cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGCTT‑T‑CTC‑GC‑GTTCCCCCG‑CCC‑GCGC‑TTTTACCCTT‑CCCCGGC‑GGGCCGC‑AAACTCCCCCCG‑GTTC‑ttt                                                                                                                                                                                                                                  >  1:105349/1‑239 (MQ=255)           cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTAAGCTGCCAGCCAGCAATTGGC‑GGGCTT‑T‑CTC‑GC‑GGTCCCCCG‑CCC‑CTGC‑TTTTAACGCA‑CACCGGG‑TGGCAGT‑AAACCCCCCCCC‑ACTC‑ttt                                                                                                                                                                                                                                  >  2:779984/1‑239 (MQ=255)             gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑TTC‑GC‑GCTCCGCCG‑CCC‑CTGC‑TTTAAACGTA‑CCCAAGC‑GGGCCGT‑AACCTCCGCCCC‑AGTC‑tttt                                                                                                                                                                                                                                 >  1:393417/1‑237 (MQ=255)             gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCt                                                                                                                                                                                                                                                  <  1:1097723/223‑1 (MQ=255)               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTcccc                                                                                                                                                                                                                              <  1:1067846/239‑1 (MQ=255)               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTccc                                                                                                                                                                                                                               <  2:316400/238‑1 (MQ=255)               ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTccc                                                                                                                                                                                                                               <  1:1269034/238‑1 (MQ=255)               tgccACACCTTCAACTGAAGAAGTAATCCCGGTCTGATAGAGCTGGCGATCGGTCGCGATAATGGTGGTCTGCTCTTTTGCTATTTGCTTCAAGACCGTGTTTA‑CCGACGCCTGGGTTTGCCACACCCAGTTCA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTGcccc                                                                                                                                                                                                                              <  2:218677/237‑1 (MQ=255)                  cacaCCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGTC‑GGGGTT‑G‑CTC‑GC‑GGTCCCCCG‑CCC‑CGGC‑TTTAACCTCA‑CCCACGG‑GGGGAGC‑TACATCCGCCCG‑GTTC‑TTTCCCCa                                                                                                                                                                                                                             >  2:161552/1‑237 (MQ=255)                    cacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCCCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGG‑AAACTCCGCCCG‑ATTC‑TTTCCCCAGa                                                                                                                                                                                                                           >  1:1117300/1‑237 (MQ=255)                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACTCTGCCAGCCAGCCATTGGC‑GGGGTT‑T‑GTC‑CC‑GGTCCGGCGCCCC‑GTCC‑TTTAACCGGC‑CCCGGCG‑GGGGACA‑AACCCCCGCCGC‑TTTC‑TTTccccctcactgc                                                                                                                                                                                                                      >  2:751771/1‑229 (MQ=255)                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCCGGC‑GGGCAGT‑AACCCCCGCCCC‑GTTT‑TTTCCCCCGATATCgag‑a                                                                                                                                                                                                                  >  1:1261303/1‑239 (MQ=255)                          cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑T‑C‑C‑GCGGTTCCGCCC‑CCC‑CTGC‑TTTTACCGCACCCCCGGG‑GGGCACG‑ACCCCCCCCCCG‑ATCT‑TTCCCCCAGaattgtga                                                                                                                                                                                                                    >  2:588789/1‑232 (MQ=255)                           aaCTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAGTGGAGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑ATCGGC‑GTTCC‑CCG‑CCC‑CCTC‑TTTTACCATA‑CCCCGGCGGGGCGTT‑ACCCCCCGCCCC‑ATCT‑TTTccccgcaatcgcgat                                                                                                                                                                                                                   >  2:510482/1‑228 (MQ=255)                                 agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GGTCCCCCC‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCCGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCCGATATCGAG‑AT‑GCCAg                                                                                                                                                                                                           >  1:556770/1‑238 (MQ=255)                                     gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑CTC‑GC‑GTTCCGCCC‑CCC‑GTGC‑TTTAACCGCA‑CCCATGC‑GGGCCGT‑AACCTCCGCCCC‑CTTC‑TTTTCCCCGAGATCGTG‑AT‑GCCAGCccg                                                                                                                                                                                                       >  1:1159706/1‑238 (MQ=255)                                     gTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑CGTTCGCCG‑CCC‑CGGC‑TTTAACCCCA‑CCCCGGC‑GGGGAGT‑AACCTCCCCCCC‑CGTC‑TTTTCCCACATATCCAG‑AT‑GCCAGcccc                                                                                                                                                                                                       >  2:1233532/1‑237 (MQ=255)                                             gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCCCCG‑CCC‑GGGC‑TTTTACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCCCCCG‑ATTT‑TTTCCCCAGATATCGAG‑AT‑GCCCGCCCGCCGTT‑TAg                                                                                                                                                                                              >  1:615416/1‑238 (MQ=255)                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑GTC‑GC‑GCTTCGCCC‑CCC‑GTGCTTTTACCCGTC‑CCCAGCG‑GGGCAGG‑CA‑CTCCGCCCG‑GCTCTTTTCCCCATAT‑TCGAG‑AT‑GGTCGCCCGCCGTT‑atac                                                                                                                                                                                             >  2:628952/1‑235 (MQ=255)                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑CTC‑GC‑GGTCCGCCG‑CCC‑GGGC‑TTTTACCCTA‑CCCCGGG‑GGGGCAT‑AAACTCCCCCCC‑GTTC‑TTTTCCCCCATATCGGG‑GT‑TGCAGCCCCCCGTT‑ACGc                                                                                                                                                                                             >  2:1309124/1‑238 (MQ=255)                                              ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑CC‑GTTCCGCCC‑CCC‑GGGC‑TTTAACCCTA‑CCCACGG‑GGGCCGG‑AACCCTCCCCCG‑GTTC‑TTTTCCCCCATATTGTGAAT‑GCAAGCCCCCCGTT‑AAGc                                                                                                                                                                                             >  1:4625/1‑239 (MQ=255)                                                       aGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCCCCG‑CCC‑GTGC‑TTTTACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCCCCCG‑ATTC‑TTTCCCCCGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCCAAAATAc                                                                                                                                                                                    >  1:86242/1‑238 (MQ=255)                                                             cGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCCTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTTCCCCGATATCGAG‑AT‑GCCAGCCCCCCGTT‑AAGCCAAAAGTACCGTTGg‑t                                                                                                                                                                            >  1:1286388/1‑239 (MQ=255)                                                              gATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑ACCCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGg‑tg                                                                                                                                                                           <  1:241165/239‑1 (MQ=255)                                                                        gcaagtGGTGTTCGCCGCTTTTTCTATTGCCCGAACGCCCGTGTTAG‑ACCCCGCCTGGGTTGCCCCCCCCCAGTACG‑GGCGGGCTGCCCGGCCAGCCAGCATTTGGC‑GGGTTT‑G‑CGC‑GC‑GTTCCGCCC‑CCC‑GTGC‑TTTACCCGTA‑CCCAGGC‑GGGCGGT‑ACCCTCCGCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑Cggg                                                                                                                                                                <  2:1002013/233‑1 (MQ=255)                                                                        gATAATGTTGTCCTCCTTTTTTTCTTTTTGCTGCAAGACCGTTTTTA‑ACGCCGCCGGGGTTGGCCCCCCCCAGTACG‑GGCGGGCTACGCTGCCAGCCAGCAATGGGG‑GGGTTT‑G‑CTC‑GC‑TTTCCGCCG‑CCG‑GTGC‑TTTAACCGTC‑CCCAGGC‑GGGCAGT‑ACCCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑Cggg                                                                                                                                                                <  2:1261303/238‑1 (MQ=255)                                                                              ggTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACCCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGTTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTC‑CCCAGGC‑GGGCAGT‑ACCCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑g                                                                                                                                                        <  1:1130764/239‑1 (MQ=255)                                                                                  ttCTGCTTTTTTTCTATTGTCGGCAAGCCCGTGTTTA‑CCGCCGCCGGGGTTCGCCACCCCCAGTACA‑GGCGGGCTACGCTCCCGCCCAGCAATGGGG‑GGGTTT‑G‑CTC‑GC‑GTCCCGCCG‑CCC‑GTTC‑TTTACCCGTC‑CCCAGGG‑GGGCAGT‑ACCCTCCCCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTAcc                                                                                                                                                    <  2:615416/239‑1 (MQ=255)                                                                                    cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCCGT‑AACCTCCGCCCG‑ATTC‑TTTTCCCCGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCt                                                                                                                                                   >  1:860516/1‑238 (MQ=255)                                                                                            tttgcgtttGCTGAAAG‑CCGTTTTTACCCGCCGCCGGGGTTGGCCACGCCCA‑CACACGGGGCGTTACGCTGCCAG‑CAGGCAGTGGT‑GGGTTTGG‑CTC‑GT‑GTCCCGCGG‑CCG‑ATGC‑TTGACCAGGA‑CCATGGC‑GGG‑AGTGACCCTCCGCCCA‑CTCT‑TTTCCCCGAATACCGTG‑AG‑TCCTGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGa                                                                                                                                          <  2:763338/233‑1 (MQ=255)                                                                                            tttCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GCGC‑TTTTACCGTA‑CCCAGGC‑GGGGCGT‑AACCCCCGCCCG‑CTTC‑TTTCCCCCGATATCCAG‑AT‑GCCCACCCGCCGGT‑AAGCCCAAAGTAACCCTGG‑TGTA‑CCA‑CCGGGCGGCC‑GTACCCTCGGccacg                                                                                                                                          >  2:335867/1‑236 (MQ=255)                                                                                                 aTTGGCTGAAAGACCGTTTTAC‑ACGCCGCCTGGGTGCGCCACCCCCAGTACG‑GGCGGGCTACCCTGCCCGCCAGCAATGGGC‑GGGTTT‑G‑CTC‑GC‑GTCCCGCCC‑CCC‑GTGC‑TTTAACCGTC‑CCCGGGC‑GGGCAGT‑ACCCCCCGCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGtt                                                                                                                                     <  2:849228/239‑1 (MQ=255)                                                                                                    tgctgcAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTTCCAGCCCGCAATTGGC‑GGGGTT‑T‑CTC‑GC‑CTTCCGCCC‑CCC‑CGGC‑TTTAACCCTA‑CCCCGGC‑GGGCCGG‑AACCTCCCCCCC‑CTTC‑TTTTCCCCGATATCTCG‑AT‑GCCCGCCCCCCCCT‑AAACCCAAAGTAACGCTTG‑GGTG‑CCC‑CCGGCCCGCC‑CTACCCGCGCCCCGCTCGTTCac                                                                                                                                  >  2:1185461/1‑238 (MQ=255)                                                                                                            gACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑GTC‑CC‑GTTCCCCCG‑CCC‑CTGC‑TTTAACCCTA‑ACCCGGC‑GGGGCGT‑AACCTCCCCCCG‑CTTC‑TTTTCCCAGATATCCAG‑AT‑ACCCGCCCCCCGTT‑AAGCCCAAAAATCCCTTGG‑TGGA‑CCC‑CCGGCCGGTC‑GGAACTGCGGCCGGAGCGGTCAGAG‑CAAACg                                                                                                                         >  2:372653/1‑239 (MQ=255)                                                                                                              ccGTGTTTA‑ACGCCGCCGGGTTTTGCCACCCCCAGTACG‑GGCGGGCTCCGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTACCCGTC‑CCCAGGC‑GGGCAGT‑ACCCTCCGCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGc                                                                                                                       <  1:682587/239‑1 (MQ=255)                                                                                                                 tgtTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑G‑GCC‑CC‑GGTCCCCCG‑CCC‑CGGC‑CTTATACGCA‑CCCCGGC‑GGGCCGG‑AACCCCCCCCCC‑GATC‑TTTCTCCAGATATCCCG‑AT‑GGCCCCCCCCCTCT‑AAACCAAAAGTACCCTTGG‑TGTG‑CCA‑CAGGCCGGGC‑GTCACCGTCGCCGGGTAGTTCCAAGCCAAACGGcc                                                                                                                      >  2:180379/1‑238 (MQ=255)                                                                                                                         cgccgcCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑CTC‑GC‑GGTCCCCCC‑CCC‑CTGG‑TTTAACCCTA‑CCCCGGC‑CGGCAGG‑AACCTCCCCCCC‑CTTC‑TTTTCCCCCATATCCTG‑ATGGCCACCCCCCCGTT‑AACCCAAAA‑TACCGTTTG‑GTTA‑CCC‑ACGGCCGGTCGGTCCCCGC‑GCCGGGTCGTTCACGG‑CACACAGcccccctaaa                                                                                                              >  2:671070/1‑232 (MQ=255)                                                                                                                              ccTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTT‑G‑TCC‑CC‑GTCCGCCCC‑CCC‑GGGC‑TTAACCCGAA‑CCCCGGC‑GGG‑AGT‑GACCTCCCCCCGCTTTC‑TTTCCCCCAATA‑CCAG‑GT‑GCCCGCCCCCCCTTAAAGCCCAAA‑TACCCTTGGTTGTA‑‑CA‑CCGGCCCGTC‑GTACCCTCGGCCTGGCTGCTCAGAA‑CAAACACCCCCAATAAACCTTc                                                                                                         >  2:1303864/1‑238 (MQ=255)                                                                                                                              ccTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑T‑CCC‑GC‑GGTCCGCCG‑CCC‑GGGC‑TTTTACCCTA‑CCCAGGC‑GGGGAGT‑AAACCCCGCCCG‑GTTC‑TTTCCCCCGATATCGAG‑AT‑GCCAGCCCCCCCGT‑AAGCCAAAAATAACCTTTG‑GGTA‑CCA‑CGGGGCGGGC‑GTACCCCCGGCCGGAGCGTTCAGAG‑GCAACGGCCCCCCTTAACCTTCt                                                                                                        >  2:505640/1‑239 (MQ=255)                                                                                                                              ccTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCCCCG‑CCC‑CCGC‑TTTTACCCGA‑CCCACGC‑GGGGCGT‑AAACTCCCCCCC‑GTTC‑TTTTCCCCGATATCGAG‑AT‑GCCCACCCCCCGCT‑AAGCCCAAAATAACGTTGG‑TGTA‑ACACCGGGCCGG‑C‑GTACCCGCGGCCCGGACGGTCAGAG‑CAAACGGGCCCCTTTAACCTTCt                                                                                                        >  2:496678/1‑239 (MQ=255)                                                                                                                                   ttttGCCCCCCCCATAACA‑GGCGGTCTACGCTGCCAGCCAGCAATTGGG‑GGGTTT‑C‑CTG‑GC‑GTCCCGCCC‑CCC‑GTGC‑TTTAACCGTC‑CCCGGGG‑GGGAAGT‑ACCCTCCCCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AG‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCg                                                                                                    <  2:544568/237‑1 (MQ=255)                                                                                                                                        ccACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑G‑CTC‑CC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCCTA‑CCCAGGC‑GGGCAGT‑AAACTCCCCCCG‑ATTC‑TTTCCCCCGATATCGAG‑AT‑GCCCACCCGCCCTT‑AAACCCAAAGTACCCTTTG‑GGTACCCC‑CGGGCCGG‑C‑GGACCTGCGGCCGGGTCGTTCCGAG‑CAACAGGCCCCACTAAGCCTTCTGCCGCCAttt                                                                                              >  2:230716/1‑239 (MQ=255)                                                                                                                                             cccAGTACA‑GGCGGTCTCCGCTGCAGGCCAGAAATTGGC‑GGGTTT‑C‑CTC‑GT‑GTTCCCCCG‑CCC‑GGTT‑TTAAACCGTC‑CCCAGGG‑GGGCAGT‑ACCCTCCCCCCG‑ATTT‑TTCCCCCAGATAGCGAG‑AG‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                                                                         <  2:1267132/239‑1 (MQ=255)                                                                                                                                                        gCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑T‑GTC‑GC‑GGTTCGCCG‑CCC‑GTGC‑GTTTACCGTC‑CCCAGGC‑CGGGAGT‑AACCTCCCCCCG‑CTTT‑TTTTCCCACGGATCGAG‑AT‑GCCCGCCCCCCGTT‑AAACCCAAAATAACGTTGG‑GGGT‑CCA‑CGGGCCGGGC‑GTGCGCGCGGCCCGGGCGGTCACGG‑CCAACGGGCCCCCTAACAACTCTCTCGACCTTTTTTGCTGCCCCatt                                                                                >  2:168018/1‑237 (MQ=255)                                                                                                                                                        gCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTT‑T‑CTC‑GC‑GGTTCGCCG‑CCC‑CTGC‑TTTTACCGTA‑CCCCCGC‑GGGGAGT‑AACCTCCGCCCC‑ATTC‑TTTTCCCCGGTATCCCG‑AA‑GCCAACCCGCCGTT‑AAGCCAAAAGGACCGCTGG‑TGTA‑CCA‑CGGGCCGGGC‑GTGCCTGCGGCCCGGTCCGGCTGAG‑CCAACTGCCCCCCTAAGCCTCCTCCCGACCTTTTTTTCCGCTccctc                                                                                >  2:1023423/1‑235 (MQ=255)                                                                                                                                                                        gccagcATTGGGG‑GGGTTG‑C‑CTC‑GC‑TTTCCGCCG‑CCC‑TTGC‑TTTACCCGTC‑CCCAGGC‑GGGCATT‑ACCCCCCCCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGAt                                                               <  1:1233532/239‑1 (MQ=255)                                                                                                                                                                                   tgggggTTT‑GCCCC‑GC‑GTTCCGCCC‑CCC‑GTGC‑TTAAACCGTC‑CCCAGGG‑GGGAAGA‑ACCCCCCGCCCA‑ATTT‑TTCCCCCAAAAATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa                                                        <  1:163008/234‑1 (MQ=255)                                                                                                                                                                                                      gTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTACGCCTTCTGCCGACATTTTTTGCCGCCCCATCTCCGCCGCAAAGTCCACCTGCGGACCATCCTGACTGGcg                                         >  2:1094475/1‑235 (MQ=255)                                                                                                                                                                                                         ccgccc‑CCC‑GTGC‑TTAAACCGTA‑CCCAGGC‑GGGAAGT‑AACCCCCCCCCG‑ATTC‑TTCCCCCAGATACCGAG‑AT‑GCCGGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCAGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGcc                                   <  1:335867/239‑1 (MQ=255)                                                                                                                                                                                                                       c‑TTTAACCGTC‑CCCAGGC‑GGGCAGT‑ACCCCCCGCCCG‑ATTC‑TTCCCCAAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCt                       <  1:505640/239‑1 (MQ=255)                                                                                                                                                                                                                        ttttAACCGTC‑CCCGGGC‑GGGCAGT‑ACCCCCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AG‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCt                       <  1:230716/238‑1 (MQ=255)                                                                                                                                                                                                                                gTA‑CCCAGGC‑GGGCAGT‑ACCCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCACAACTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTACAGAG‑CAACCGGCCCCATTAAGCCTTCCGCCGCCCTTTTTTGCCGCCCCATATCCGCCGAAAAGTCGATCTGCGGACCACCCTGAGCGGCACCTGCCTTCGCCGCGCCTTCaactc                >  2:166032/1‑237 (MQ=255)                                                                                                                                                                                                                                gTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                >  2:934381/1‑238 (MQ=255)                                                                                                                                                                                                                                catcccaagcc‑GGCCAGA‑ACCCCCCGCCCG‑TTTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                <  1:1224294/229‑1 (MQ=255)                                                                                                                                                                                                                                 tA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                >  1:715508/1‑237 (MQ=255)                                                                                                                                                                                                                                 tA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                >  2:741399/1‑237 (MQ=255)                                                                                                                                                                                                                                  a‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                >  2:77098/1‑236 (MQ=255)                                                                                                                                                                                                                                     ccgggC‑GGGCAGT‑ACCCGCCCCCCG‑ATTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGa           <  1:557141/236‑1 (MQ=255)                                                                                                                                                                                                                                         gC‑GGGCAGT‑ACCCTCCGCCCG‑TTTC‑TTCCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATg        <  1:372653/238‑1 (MQ=255)                                                                                                                                                                                                                                             ggCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc       >  1:278877/1‑236 (MQ=255)                                                                                                                                                                                                                                               cAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGt  >  2:842661/1‑239 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    GCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGT  >  NC_000913/2228041‑2228505

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: ‑