Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,383,219	A→C	21.7%	V248G (GTG→GGG)	ycjW ←	LacI family putative transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,383,219	0	A	C	21.7%	144.6 / 2.6	64	V248G (GTG→GGG)	ycjW	LacI family putative transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (35/11);  new base C (3/14);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 6.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTGTGGTGACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGCCCACTTTCAATTAATGCATCGGTC  >  NC_000913/1382996‑1383440                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               tttcttgtggggctgTTTTCCAGAAATTAAGGGAAATAATAACCCGCAGTTTTGTGGACCTAACTGCTGGTTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCCCCCC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTCCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACg                                                                                                                                                                                                                 <  2:36766/225‑1 (MQ=255)      gtggtgCCGTTTTTCCAGAATTTAAGGGGAATAATAACTCGCAGTTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCCCCCC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGAc                                                                                                                                                                                                            <  2:1005727/239‑1 (MQ=255)      gtggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGAc                                                                                                                                                                                                            <  1:257114/239‑1 (MQ=255)      gtggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGAT‑GCCGCA‑CGCTAAATGCCCATTAATTTCA‑GGCTGTCAGTGGCAAA‑A‑CCCGCCCCAGGCCACgggg                                                                                                                                                                                                              >  2:826559/1‑235 (MQ=255)       tggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAACATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCCC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑CCCCCC‑CGATAAATGCTCCTTAATTTCA‑GTCTTTCAGTCGCAAA‑C‑CCCCCCCCAGGCAACCGTg                                                                                                                                                                                                              >  2:699481/1‑236 (MQ=255)       tggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTg                                                                                                                                                                                                              >  1:1005727/1‑236 (MQ=255)       tggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTg                                                                                                                                                                                                              >  1:913666/1‑236 (MQ=255)       tggtgACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTg                                                                                                                                                                                                              >  1:195797/1‑236 (MQ=255)            aCGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAC‑CTGCTGCGGAAT‑AGCATTATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCC‑GGCTGTCAGTGGCAAA‑A‑ACCGCCTCAGGCAACGGTGACTgggta                                                                                                                                                                                                      >  2:562238/1‑236 (MQ=255)            aCGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCtt                                                                                                                                                                                                      >  1:715489/1‑239 (MQ=255)             cGGTTTTCCAAAAATTAAGGGGAATAATAACCTGCAGTTTTGTTGCCCTAACTCCGGGGTCGGAACATCGAGGCCCCCCGGTGCAGGCGTTAAAATAAATGGCA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGTTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGAAAC‑C‑ACCGCCTCAGGCAACGGTGACTGGCtt                                                                                                                                                                                                      <  2:143788/238‑1 (MQ=255)                           ttaatggaaataatAACTCGCAGCTTTATTGGCCTAACTCCTGCGTCGGAACATCGATGCCGACCGGTGCAGGCGTTAAAATAAATTACA‑GCGTTTCATTG‑C‑TTTAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGtgct                                                                                                                                                                                       <  1:5696/231‑1 (MQ=255)                            gcgtcaggttttttaaagggacgttttttGACCTACCCCCGGCGTCGGAACATCGAGGCCCCCCGGTGCAGGCTTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGtgct                                                                                                                                                                                       <  2:372432/215‑1 (MQ=255)                                            ccccAGTTTTTTTGCCCTAACTCCTGCGTCGGAACATCGATGCCCCCCGGTGCAGGCGTTAAAAAAAATGACA‑GCGTTTCATTG‑C‑TAAACCCCACCCC‑CGTTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAC‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTc                                                                                                                                                                      <  2:913666/236‑1 (MQ=255)                                             cGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGCCTGGCTTAATAATTGCCGTCCTGCTTTTAATGCAGTTTCa                                                                                                                                                                     >  1:225631/1‑239 (MQ=255)                                             cGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCa                                                                                                                                                                     >  1:22467/1‑239 (MQ=255)                                              gagtgtttggttcccaaaCCCCGGCGGCGGAACTCGGATGCCCCCCGGGCCAGGCTTTAAAATAAAGGGCA‑TGGTTTCATTG‑C‑TATAACCCCCCCC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCAt                                                                                                                                                                    <  2:366937/223‑1 (MQ=255)                                               cAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTT‑AAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTTTCAGTGGCAAA‑C‑CCCCCCTCAGGCAACGGTGACTTGCTTAATAATTTCCCGGCTGCTTTTAATGCAGTTTCATgg                                                                                                                                                                  >  1:113246/1‑239 (MQ=255)                                                aGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAC‑C‑ACCGCCTCAGGCAACGGTGATTGGCTTAATAATTTCCGTGCTGCTTTTAATGCAGTTTCAAgg                                                                                                                                                                  >  2:360489/1‑239 (MQ=255)                                                aGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAATTTCATgg                                                                                                                                                                  >  1:317731/1‑239 (MQ=255)                                                           ctcatttgtggggTCGGAACTCAGCTGCCCCCCGGGGGGGGGTTTAAAATAAATGAAC‑GTTTTTCTTTG‑T‑TAT‑‑CCCCCCACGCCCTAATCGGCTGGGGAAT‑AAGAATTTTTTT‑CT‑TGCCC‑GCCGAA‑CGAAAAGGGCTAATTATTTCGG‑GGCTTTCAGGGGCAAC‑C‑CCCGCCGCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGcc                                                                                                                                                       <  2:225631/227‑1 (MQ=255)                                                               aCTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATg                                                                                                                                                                             >  1:372432/1‑213 (MQ=255)                                                                      cGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAAt                                                                                                                                              >  1:894483/1‑237 (MQ=255)                                                                           gAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATtc                                                                                                                                        >  1:843462/1‑238 (MQ=255)                                                                                tgtttccccccgggcgCAGGTGTAAAAATAATTGAAC‑GGTTTTCATTT‑C‑TATACCCCCCCGC‑CGAAAC‑TTGGCGCGAAAT‑A‑CAATTTTTTTCCT‑CGGCC‑GCCGCA‑CAAAAATGGCTAATTATTTTCA‑GCCTGT‑ATGGGCAAACC‑CCCCCCTCAGGCAACGGTGACTGACTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATACCCGCCATCAACAATCCATTCATCACg                                                                                                                                  <  2:105855/224‑1 (MQ=255)                                                                                tggATGCCGCCGGGGCCAGG‑GTTAAAATAAATGGCA‑GGTTTTTAATGTC‑TATACCCCACCAC‑CGCAAC‑TTGCTGGGGAAT‑GAAAATATTTTT‑CT‑CGGCT‑CCCCCA‑AGATAAATGATCATTAATTCCA‑GTCTTACGGTGGAAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTATCATGGGTATAACCGCCATCAACAATCCATTCATCACg                                                                                                                                  <  1:54058/237‑1 (MQ=255)                                                                                tCGATGCCGCCGGGTGCAGGCGTTAAAATAAATGTAA‑GCTTTTCATTG‑C‑TAAACCCCACCAC‑CGTTAA‑CTGCTGGGGAAT‑AGCAATATTTTC‑CT‑CTCCT‑GCCGAA‑CGATAAATGTTCATTAATTTCA‑GTCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACg                                                                                                                                  <  1:895011/239‑1 (MQ=255)                                                                                cggtggccgccGTGTGCTGGAGTTAAAATAAAGTACA‑GTTTTTCATTG‑C‑TATCCCCCACACC‑CCCTAA‑CTTGTGGAGAAT‑AACAATTTTTTC‑CT‑CTGCT‑GCACCA‑CGATACATGCCTATTATTTTCA‑GGCGGTCGGTGGCAAC‑C‑CCCGCCTCAGGCAACGGTGACTGGTTTAATAATTGCCGTGCTGCTTTTAAGGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACg                                                                                                                                  <  2:633911/234‑1 (MQ=255)                                                                                  gttcccccccgtgccgggCTTTAAAATAAATG‑CA‑GGGTTTTTTTT‑T‑TATAACCCCCCCC‑CGCAAA‑ACGTGTGGGAAA‑AGAAATTTTTTT‑CTCCTGCC‑GCC‑CA‑CGAAAAAAGCTCAATTATTTCAGGGTTGACGGGGGCAAC‑C‑CCCGCCTCAGGAAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCAGGGGTATAACCGCCATCAACAATCCATTCATCACGCa                                                                                                                                <  2:1071405/223‑1 (MQ=255)                                                                                  gATGCCGCCGGGTGCAGGCGTTAAAATAAATGCCA‑GCTTTTCATTG‑C‑TATAACCCACCCC‑CGTTAA‑CGGCGGGGAAAA‑AGAAATATTTTT‑CT‑CTGCT‑GCCGCC‑CGAAAAATGCTAATTAATTTCG‑GGCTGCCAGTGGAAAC‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCa                                                                                                                                <  1:173697/239‑1 (MQ=255)                                                                                  gATCCCGCCCGGT‑CACGGGTTAAAAAAAATTGCACACTTTTCTTTC‑T‑TATACCCCCCCCC‑CGTTAA‑CTGCGGGGAAAT‑AGAAATTTTTTG‑CC‑CTGCT‑GCCGCA‑GAATAAATGCTCTTTATTTTCG‑GGTTGTAGGTGGAAAC‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCa                                                                                                                                <  1:375020/239‑1 (MQ=255)                                                                                   aTGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATc                                                                                                                                             >  1:503164/1‑225 (MQ=255)                                                                                                 cggcttaaaaaaaaaTGACA‑GCTTTTCTTTG‑T‑AATCCCCCCCCCC‑CTCTAC‑CTTCTG‑GGAATAAGAAATTTTTTT‑CG‑CTGCC‑CCGACA‑CGATAATTGCTTATTAATTCCA‑GGCTTTCGGTGGCAAC‑A‑CCCCCCCCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGTTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATgcg                                                                                                                <  2:195797/232‑1 (MQ=255)                                                                                                   gCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGGATAACCGCCATCAACAATCCATTAATCACGCCC‑TGCAAAATTATgcgc                                                                                                               >  2:891731/1‑238 (MQ=255)                                                                                                        aaaaTAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATACCCGCCATAAACAATCCATTCATCAC‑CACATGCAATATTATCCCCAGCCaa                                                                                                         >  2:608835/1‑238 (MQ=255)                                                                                                            ataaaggcactcttttCATTG‑C‑TAAACCCCCCCCC‑CGTTAA‑CGGCTGCGGAAT‑GGCAATATTTTT‑CC‑CTGCT‑GCCGCA‑CAATAAATGCTCATTAATTGCA‑GGCTGTCAGTGGAAAC‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGct                                                                                                      <  1:670695/226‑1 (MQ=255)                                                                                                             aaaaaggacacGTTTTATTT‑CTTATAACCCCCCCC‑CGTAAA‑CGGC‑GCGGAAA‑AAGAAAATTTTT‑CC‑CCCCC‑GCCCGACCGAAAATTGCTATTTAATTGCA‑GGCTTTCAGGGGAAAA‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGct                                                                                                      <  1:360489/229‑1 (MQ=255)                                                                                                                    a‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTACTAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCACCAATCCATTAATCACGCCA‑TGCAATATTATGCCCAGCCAGGCCCTGCTtata                                                                                             >  2:122410/1‑239 (MQ=255)                                                                                                                       cGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAATTGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTCATGCAGTTTCATGGGTATAACCTCCATCAACACTCCATTCATCCCGCTC‑TCCATTATTATGCCCCGCCAGGccccttttaaaac                                                                                           >  2:762929/1‑227 (MQ=255)                                                                                                                         tttCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAAtatt                                                                                                                     >  1:633673/1‑211 (MQ=255)                                                                                                                                  c‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTAATCCCGCAC‑TGCAATATTATGCGAAGCCAGCCTCTGCTTATAACCATTAAcccg                                                                                 >  2:125838/1‑239 (MQ=255)                                                                                                                                             ccac‑cGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATCAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGACATCACCAATCCATTCATCACGCAC‑TGCAATATTATACGCAGCCAAGCTCAGATTATAACCATTAACGCGACCAACCGaa                                                                       >  2:995571/1‑239 (MQ=255)                                                                                                                                             ccac‑cGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCCC‑TGCAAAATTATGCCCAGCCAGGCCCTGCTTATAACCATTAACGCGACCAACCaaa                                                                       >  2:718478/1‑239 (MQ=255)                                                                                                                                             ccac‑cGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGaa                                                                       >  1:28139/1‑239 (MQ=255)                                                                                                                                              cac‑cGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATgcg                                                                                                                >  1:795950/1‑197 (MQ=255)                                                                                                                                                    cTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCCCTGCTTTTAACCATTAACCCGATCAACCGAAACATga                                                                 >  2:147844/1‑238 (MQ=255)                                                                                                                                                    cTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGAAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGACAAACCGAAACATgg                                                                 >  2:859080/1‑239 (MQ=255)                                                                                                                                                     cctaactttgtccaaaatagtatatGTTT‑CT‑CGTCTGGCCGCA‑AGTAAAAGCCTCATAAATTTCG‑GGCTGTATGT‑GCAAA‑CTCCCCCCTAGGGCACCGGTGATGGGCTTAATACTTGCCGGTCTGCTTTTAATGCATTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACAt                                                                   <  2:979231/218‑1 (MQ=255)                                                                                                                                                       a‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACa                                                              >  1:704354/1‑239 (MQ=255)                                                                                                                                                             tgcGGAAA‑AGAAATTTTTTT‑CT‑CTGTG‑GCCCAA‑CGATAATTGCTTATTAATTACG‑GGTTGCCAGGGGAAAA‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGTTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAg                                                         <  1:122410/239‑1 (MQ=255)                                                                                                                                                                            tttttt‑CC‑CTCCC‑GCCCCA‑AAATAAATCCTCTTAAATTTCG‑GGTTGTCGGGGGAAAA‑C‑CCCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCATTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATgcaggc                                           <  2:715489/238‑1 (MQ=255)                                                                                                                                                                                   ct‑ctGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATaa                                                                                            >  1:970802/1‑184 (MQ=255)                                                                                                                                                                                         ct‑gcCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTACCGCGATCAACCGAAACATGGACATCAACCGGAGCATGCAGGCAGGCAATAttt                                >  1:836277/1‑239 (MQ=255)                                                                                                                                                                                         ct‑gcCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACACGATCAACCGAAACATGGCCATCAAGCGGGGCATGCAGGCAGGCATTAttt                                >  1:744329/1‑239 (MQ=255)                                                                                                                                                                                             ccGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCtt                                                                                                                                                                                   >  2:533831/1‑90 (MQ=255)                                                                                                                                                                                             ccGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCtt                                                                                                                                                                                   <  1:533831/90‑1 (MQ=255)                                                                                                                                                                                              cGCC‑CGATAATTGCTTATTAATTTAG‑GGCTGGCGGTGGAAAA‑C‑CCCGCCGCAGGCAACGGTGCCTGGCTTAATAATTGCCGTGCTGCTTTTAATGCATTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGa                             <  1:30618/238‑1 (MQ=255)                                                                                                                                                                                                  acaaTAAATGCATATAAATTTGG‑GGTTGCCAGGGGAAAA‑C‑CCCGCCTCAGGCAACGGTGCCTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGcc                         <  1:718478/236‑1 (MQ=255)                                                                                                                                                                                                               aTTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACACTCCATTCATAACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAAAATGGCCCTCACACGGAGCACGCAGGCAGGCAGTATTTTGTTGCCCACTTTCAATTaa            >  1:77751/1‑239 (MQ=255)                                                                                                                                                                                                               aTTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGCCCACTTTCAATTaa            >  1:528065/1‑239 (MQ=255)                                                                                                                                                                                                                ttAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATgcg                                                                                                                >  1:840158/1‑138 (MQ=255)                                                                                                                                                                                                                   ttttGA‑GGTTGTAGGTGGAAAC‑C‑TCCGCCTCGGGAAAGGGGGATTGGCTTAATAATTGCCGTGCTGTTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGCCCACTTTCAATTAATGCa        <  2:843462/238‑1 (MQ=255)                                                                                                                                                                                                                          ggCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACTCGATCAACCGAAACATGGACATCAAACGGGGCAGACAGGCAGGCAATATTTTGATGCCCACTTTCAATTACTGCATCGGTc  >  1:1032845/1‑239 (MQ=255)                                                                                                                                                                                                                                       aaa‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGCCCACTTTCAATTAAt           >  1:314118/1‑217 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               ATTCTGTGGTGACGGTTTTCCAGAAATTAAGCGGAATAATAACTCGCAGCTTTGTTGACCTAACTCCTGCGTCGGAACATCGATGCCGCCCGGTGCAGGCGTTAAAATAAATGACA‑GCGTTTCATTG‑C‑TATAACCCACCAC‑CGCTAA‑CTGCTGCGGAAT‑AGCAATATTTTT‑CT‑CTGCT‑GCCGCA‑CGATAAATGCTCATTAATTTCA‑GGCTGTCAGTGGCAAA‑C‑ACCGCCTCAGGCAACGGTGACTGGCTTAATAATTGCCGTGCTGCTTTTAATGCAGTTTCATGGGTATAACCGCCATCAACAATCCATTCATCACGCAC‑TGCAATATTATGCGCAGCCAGGCTCTGCTTATAACCATTAACGCGATCAACCGAAACATGGACATCAAGCGGAGCATGCAGGCAGGCAATATTTTGATGCCCACTTTCAATTAATGCATCGGTC  >  NC_000913/1382996‑1383440

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