Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,228,256	A→C	22.2%	pseudogene (354/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,256	0	A	C	22.2%	79.1 / 4.3	61	pseudogene (354/1287 nt)	mdtQ	pseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base A (32/12);  new base C (4/12);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 2.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_000913/2228024‑2228491                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAgg                                                                                                                                                                                                                                                   <  1:187512/239‑1 (MQ=255) gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAgg                                                                                                                                                                                                                                                   <  1:1304874/239‑1 (MQ=255)   tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGc                                                                                                                                                                                                                                                  <  2:43019/238‑1 (MQ=255)      cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATTTAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCGCTTTTTCTATTTGCTGCAAGACCGGGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTC‑CCCAGGC‑ggcg                                                                                                                                                                                                                                             <  2:242440/239‑3 (MQ=255)                  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                                                                                                                                                                                                                 <  1:1269034/238‑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)                                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)                                                                                                                                                                                                                                                  tA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:741399/1‑237 (MQ=255)                                                                                                                                                                                                                                                  tA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  1:715508/1‑237 (MQ=255)                                                                                                                                                                                                                                                   a‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:77098/1‑236 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTA‑ACGCCGCCTGGGTTTGCCACTCCCAGTACA‑GGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTT‑G‑CTC‑GC‑GTTCCGCCG‑CCC‑GTGC‑TTTAACCGTA‑CCCAGGC‑GGGCAGT‑AACCTCCGCCCG‑ATTC‑TTTCCCCAGATATCGAG‑AT‑GCCAGCCCGCCGTT‑AAGCCAAAAGTACCGTTGG‑TGTA‑CCA‑CGGGCCGGTC‑GTACCTGCGGCCGGATCGTTCAGAG‑CAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_000913/2228024‑2228491

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