Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,241,337	A→C	23.6%	G245G (GGT→GGG)	cvrA ←	putative cation/proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,241,337	0	A	C	23.6%	159.3 / 8.4	75	G245G (GGT→GGG)	cvrA	putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (40/14);  new base C (2/18);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 7.14e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAATCCAGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGCTCTCATGATGCTGGATCATCGCAATTAGGGTAATCGTCAG  >  NC_000913/1241104‑1241570                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            aGCAATCCAGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGAGGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCCGAATACCGc                                                                                                                                                                                                                                         >  1:295628/1‑238 (MQ=255)    aaTCCAGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTa                                                                                                                                                                                                                                       <  1:1002879/237‑1 (MQ=255)       cgggCAAAGGCCGAAAGGGGAGGGGGGAAAAATGCCAGCCATGCGGACAATCAGAGGCGCGGAATGGCTATGGGCAGCAGATCGGTTGGGTTAACCAACAGCCGCAGGACCAAGAAAGTGGCGATTTGCGCCA‑GCCAGGTGAGGCCGTCGAAATGTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCACAAAGATACA‑CCGGCAGAATGCCGCTGCGTTc                                                                                                                                                                                                                                  <  2:467417/236‑1 (MQ=255)       ccAGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTc                                                                                                                                                                                                                                  <  1:1281299/239‑1 (MQ=255)       ccAGAAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTc                                                                                                                                                                                                                                  <  1:1129670/239‑1 (MQ=255)         aGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGCGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAACAGGGCGA‑TT‑ACCCAGCAGACACCCGCCCAGATACA‑CCGCCAGAATACCGCTACCTTcc                                                                                                                                                                                                                                 >  1:1138384/1‑238 (MQ=255)            aaaGACCGAAAGGGGACGGGCAAAGAATATCAACCATGCGGAAAAAATGAGCCCCGGAAGTGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACAGGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAgcg                                                                                                                                                                                                                             <  2:272297/239‑1 (MQ=255)              gtcagattgggttaagagaaaggaaTATCACCCAAGGGGAAAAAAGGAGCGCCGGAAGGGCAAAAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACAGGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAgcgc                                                                                                                                                                                                                            <  2:575735/216‑1 (MQ=255)               tACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCag                                                                                                                                                                                                                          <  2:761772/239‑1 (MQ=255)                    aaaGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGACAAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCCCAGG‑AGtt                                                                                                                                                                                                                    >  2:9126/1‑239 (MQ=255)                     aaGGGGACGGGCAAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑a                                                                                                                                                                                                                  <  1:490646/239‑1 (MQ=255)                          gACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAACGCAGC‑AGTT‑AATgag                                                                                                                                                                                                             >  1:219224/1‑239 (MQ=255)                           cagggCGAAGAATATCAGCCAGTCGGAAAAAAGGGCGGCCGGAATGGCAATAGGAAGCAGATCGGTTGGGGTAACCGACAGCCCCAGCAACAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATgag                                                                                                                                                                                                             <  2:365528/236‑1 (MQ=255)                           aggggCGAAGAATATCACCCATGCGGACAAAATGAGCGCCGGAATGGCAAAAGGCAGCGGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATgag                                                                                                                                                                                                             <  2:689714/236‑1 (MQ=255)                           aCGGGGGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATgag                                                                                                                                                                                                             <  2:1057687/238‑1 (MQ=255)                           aCGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCCCCGT‑AGTT‑AATgaga                                                                                                                                                                                                            >  1:912191/1‑239 (MQ=255)                                        atcatcCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCCGAATACCGCTACCTTCCAGCGCAGT‑CGTT‑AATGAGAAAATCAGAAtct                                                                                                                                                                                                >  1:689714/1‑236 (MQ=255)                                         tcatcCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCCCAGCTACA‑CCGCCAGAATACCGCTACCTTCCACCGCCGA‑GGTT‑AATCGGAAAATCAGAATAcc                                                                                                                                                                                               >  1:413956/1‑238 (MQ=255)                                         tcatcCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTT‑GGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAAAACCGc                                                                                                                                                                                             >  1:422549/1‑239 (MQ=255)                                               aTGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AAt                                                                                                                                                                                                                >  1:575735/1‑215 (MQ=255)                                               aTGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCCCTGCCTTCCAGCGCGGT‑AGTT‑AATAACAAAATCGGAATACGGCCCCtt                                                                                                                                                                                        >  1:25578/1‑239 (MQ=255)                                                tGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCCGAATACCGCTACCTTCCAGCGCAGG‑AGTT‑AATGAAAAAATAAAAAAACCGCCCCTTa                                                                                                                                                                                       >  2:468449/1‑239 (MQ=255)                                                         aaTGGGCCCCGAATTGGCAATGGGCAGCGGATCGCTTGGGTAAACCAACGGCCCAAGCACCAGGAACAGGGCGATTTGCGCCA‑GCCAGGCGGGGCCGTCGAAATTGTGCAGGAGGCCGT‑GGCGGTGGCGAATAGGGCGA‑TT‑CCCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTA‑AATGAAAAAATCAGAATACCGCCACTTAATGCCAGCa                                                                                                                                                                              <  2:295628/239‑1 (MQ=255)                                                           tggccgcCGGATCGGCAAGAGGAGTCAGACCCTTGGGGTAAACCACAACCCCCAGCACCGGGAAAAGGGCGTTTTGCACAA‑GACGGGCGAGCCCTCCGAAATTTTGCAGGATCCCGT‑AGCGGTTGCAAATGGGGCGA‑TT‑ACCAAGCAGAAAACCGCAAAGATCCC‑CCGCCAGAATCCCGCTACCTCCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAAt                                                                                                                                                                            <  2:852446/235‑1 (MQ=255)                                                           aaatgtgtaaaaaggaaagaGGCAGAGGATCGCTGGGGTTAACCAACAGCCCCAGCACCAGGAACATGGGGATTGGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑CCCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCaa                                                                                                                                                                             <  1:1022290/219‑1 (MQ=255)                                                               cgcCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGCATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCAAGCAATGGat                                                                                                                                                                        >  1:205940/1‑239 (MQ=255)                                                               cgcCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGAACAGATACC‑CCGCCAGAATACCGCTACCTTCCAGCGCCGA‑AGTT‑AATGAGAAAAACAGAAAACCGCAACTTAATACAAGCAATGGaa                                                                                                                                                                        >  2:374634/1‑238 (MQ=255)                                                               cgcCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAACAGAAACCCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAACGCAGT‑AGTT‑AATGAAAAAATCAGAAAACCGCCCCTTAATGCCaaaaagtgaa                                                                                                                                                                        >  2:1050299/1‑230 (MQ=255)                                                                 ctgaaaTGGCAATAGGCAGCAGATCGCTTGGGTTACCCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGAAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAACACCGCCACTTAATGCCAGCAATGGata                                                                                                                                                                       <  1:1050299/234‑1 (MQ=255)                                                                         gaaaatggcaggaagtggtttGTGTTACCCAACAGCCCCACCCACGGAAACGGGGCGTTGTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTGTGAAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑CCCCAGCAGAAAACCGCACAGATAAC‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTA‑AATGAGAAAATCAGAATCCCGCCACTTAATGCCAGCAATGGATATAATCcg                                                                                                                                                                <  2:205940/218‑1 (MQ=255)                                                                            atgggcacggatcggCTTGGTTTAACCACCACCCCAACCCCCGGGAAGATGGCGATTTGCGCCA‑GCCGGGGGGGGCCTCAGAAATTTAGAGGGAGGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAACCCGCAAAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCggg                                                                                                                                                           <  2:862060/225‑1 (MQ=255)                                                                            aTAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCAACTTAATGCCAGCAATGGATATAATcctgtt                                                                                                                                                             >  1:543381/1‑233 (MQ=255)                                                                               ggCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCACCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGGATTAATCCGGCGGGCAg                                                                                                                                                        >  1:1081576/1‑239 (MQ=255)                                                                                   gcagATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATAAAATCCGGCGGGCAGGGCg                                                                                                                                                    >  2:291434/1‑239 (MQ=255)                                                                                   gcagATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGGCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGAGATAATCCGCCGGGCAggggg                                                                                                                                                    >  1:1155057/1‑237 (MQ=255)                                                                                      gctggtgtGGTTTACCCAACCCCCCCAGCCCCAGACAGAGGGCGATTTGCGCCA‑GCCGGGGGGGGCCGGCAAAATTTTGCAGGATGCCGG‑GGCGGTTGCGAATCGGGCAA‑TA‑CCCCAGCAGAAAACCGCATAGATACC‑CCGCCAGAATCCCGCTACCTACCAGCGCAGT‑AGTT‑AATGAGAAAACCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATg                                                                                                                                                 <  2:631774/232‑1 (MQ=255)                                                                                         cGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGAAACA‑CCGCCAGAAAACCGCTACCTTCCAACGCAGC‑AGTT‑AATGAGGAAATAAGAAAACCGCACCTTAATGCCAGCAAAGGCaataat                                                                                                                                                                   >  2:1022290/1‑218 (MQ=255)                                                                                                    cccaccaccccccgaccaaGGAACATGGCGTTTTGCGCCA‑GCCGGGCGAGGCTGTCAAAATTTTGAAGAGTGCCGT‑GGCGGTTGCAAATGGGGCGA‑TC‑ACCCAGCAGAAACCCGCACAGATACC‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTg                                                                                                                                   <  2:698304/221‑1 (MQ=255)                                                                                                    cccAACCCCCCCAGCCACAGGAACGTGGCTTTTGGCCCCA‑GCGGGGGGGGGCGGTAGAAATTTTGCAGGTTGCTGT‑AGTGGTTGCGAAGGGGGCGA‑TT‑CCCCAGCAAAAACCCGAACAAATACA‑CCGCCAGAATCCCGCTACCTCCCAGCGCAGT‑AGTA‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTg                                                                                                                                   <  2:1155057/238‑1 (MQ=255)                                                                                                    caccccagccccagcccCGGGAACAGGGTGTTTTG‑GCCAGGGGGGGCGGGGGCGTCGAAATTGGGAGGGAGGCCGG‑AGGGGTTGCGAAT‑GGGCGAGTT‑CCCCAACAGAAACCCCCAAAGAAACA‑CCGCCAGAATCCCGCTCCCTTCCAGCGCAGG‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTg                                                                                                                                   <  2:1276599/234‑1 (MQ=255)                                                                                                       aaacgCCCCAGCACCAGGAACGTGGCGATTTGCCCCA‑GCCAGGGGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCag                                                                                                                                <  1:374634/235‑1 (MQ=255)                                                                                                           gCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAACTTTTGCAGGATGCCGT‑AGCGGCTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CGGCCAGCATACCGCTACCTTCCAGGCCAGG‑AGTT‑AATGAAAAAATCACAATACACCCACTTACTGCCCGCAATGGATTAAATAAGGCGGCAGGGGCGATACggatcaccctttgaataaaa                                                                                                                            >  2:1119339/1‑221 (MQ=255)                                                                                                                    cgaggaaggggggggttggcgcCC‑GCGAGGGGGGGCCGCAGAAGTTGGGCAGGAT‑CCGTGAGGGGTTGCAATTGGGGCGT‑TC‑CCCAAGCAAAAACCCCAAAAGATACC‑CCCCCAGAATCCCGCTCCCTACCACCGAAGT‑ATTA‑AAGAATAAAGTCAGATGCCCGCCCCTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAAccg                                                                                                                   <  2:260843/221‑1 (MQ=255)                                                                                                                      aGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGAACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCCGAATACCGCCACTTAATGCCAGCAATGGATCTAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGCAGTAAATACccgc                                                                                                                  >  1:146764/1‑238 (MQ=255)                                                                                                                       ggAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTCCTGCCAGCAATGGATATAAACCGGCGGGCAGGGCGAGTCGGTTGATCATTTGCAGTTGTAAATGAcccgc                                                                                                                 >  1:539148/1‑235 (MQ=255)                                                                                                                            tgggCGATTTGCCCCA‑GCCGGGGGTGCCTCTGAAATTTTTGCGGGCGTCCGT‑CGGGGTTGCGAGGGCGTCGA‑TC‑CCCCAGCAAAAATCAGCACAGTGACC‑CGTCAAGACTACCGCTCCTTTCCGGC‑CAGG‑TGTTAAATAAGAAAAACAGAATCCCGCCCCTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAg                                                                                                           <  2:159241/237‑1 (MQ=255)                                                                                                                            aTGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCaa                                                                                                            >  1:859735/1‑238 (MQ=255)                                                                                                                             tGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAg                                                                                                           >  1:889056/1‑238 (MQ=255)                                                                                                                                                           tCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCCATGCGGTTGATCATTTGCAGTAGTAAATAACCCCCGCCAAGACCTATTACAATGACCAGACCAAAtttgtg                                                                             >  2:1052428/1‑235 (MQ=255)                                                                                                                                                            aaattttttGGAGG‑TGCCGT‑GGGGGTTGGCAAAGGGGCGG‑TTCCCCCAGAAAAAACCCGAACAAAACCC‑CCGCCAAAACCCCCCCACCTCCAAGCGCAGT‑AGTT‑AATGAAAAAAACAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGc                                                                            <  2:1322422/234‑1 (MQ=255)                                                                                                                                                              aaaatttttgagggtgcgcgtgCGGTTGGGAAAGGGGGGA‑TT‑CACCACGCGAAAACCCCAAAG‑TACACCCGCCAGAATCCCGCTACCTTCAAACGCAGTGTGTA‑AAGGAAAAATACAGAAACCCGCCACTTAATGCCAGCAATGGATATAACCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGc                                                                            <  2:1081576/218‑1 (MQ=255)                                                                                                                                                                                  tggggtcggaattAGGGCAA‑TA‑CCCCAGAAAAAACCCGCAAAGACACC‑CCGCCAGAATCCCGCCCCCTGCCAGCGCAGT‑AGTT‑AATGAAAAAACCAAAATCCCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACAt                                                        <  1:39663/226‑1 (MQ=255)                                                                                                                                                                                      gTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGc                                                                                                                                                                                   >  1:131717/1‑111 (MQ=255)                                                                                                                                                                                      gTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGCCAATGAAAATCAAGc                                                   >  1:801918/1‑239 (MQ=255)                                                                                                                                                                                      gTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGCCAATAAACAcccccc                                                   >  2:938304/1‑236 (MQ=255)                                                                                                                                                                                      gTTGCAAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATCCCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGc                                                                                                                                                                                   <  2:131717/111‑1 (MQ=255)                                                                                                                                                                                         gCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCag                                                                                                                                                                                                                          >  2:834482/1‑71 (MQ=255)                                                                                                                                                                                         gCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCag                                                                                                                                                                                                                          <  1:834482/71‑1 (MQ=255)                                                                                                                                                                                         gCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGAt                                                                                                                                         >  1:968864/1‑150 (MQ=255)                                                                                                                                                                                                 ggCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTa                                                                                                   >  1:1254618/1‑180 (MQ=255)                                                                                                                                                                                                 ggCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGCAAATCCAACTAATAtttcc                                         >  2:1098891/1‑235 (MQ=255)                                                                                                                                                                                                 ggCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCAAGCTAATATTGct                                         >  1:139983/1‑238 (MQ=255)                                                                                                                                                                                                  gcta‑TT‑CCCCAACAGAAAACCGCAAAAATCCC‑CCGCAAGAATCCCGCTACCTTCAAGCGCAGT‑AGTT‑AATAAAAAAATCAAAATCCCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGct                                         <  1:291434/234‑1 (MQ=255)                                                                                                                                                                                                      tttaccccAGAAAAAACCCCCACAAACACC‑CCGACAAAATCCCGTGCCCTCCCAGCGAATT‑AGTA‑AATAAAAAAACCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGCTCTCa                                     <  1:1254934/234‑1 (MQ=255)                                                                                                                                                                                                             cagcagAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTaaa                                                                                                                         >  1:822823/1‑148 (MQ=255)                                                                                                                                                                                                                   aaaaCCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGCTCTCATGATGCTGGatca                        >  1:414870/1‑239 (MQ=255)                                                                                                                                                                                                                          cacaGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCACGACCACTTACAATGCCCCGTCCAAATTGATGCAGAATATCGACAAAGAACATCCAGCTAATATTGCTCTCATGATGCTGGGAAATCGCaaa                 >  2:894562/1‑238 (MQ=255)                                                                                                                                                                                                                          cacaGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAAAATTGCCCTCATAATGCTGGATCAACGCAAt                 >  2:215469/1‑239 (MQ=255)                                                                                                                                                                                                                                     ccgccgAAATCCCGTTACCTGCAAGCGCAGA‑ATGT‑AAGGAAAAAATCAAAATGCCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGCTCTCATGATGCTGGATCATCGCAATTAGGGTAATc       <  1:370265/239‑1 (MQ=255)                                                                                                                                                                                                                                      cGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCACATTATTGCAGAATATCGACACTGAAAATCCAGCTAATATTGCTCTCATGATGCTGGATAATCGCAATTAGGGTAATc       >  1:66631/1‑238 (MQ=255)                                                                                                                                                                                                                                          aGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGCCAATAAACATCCAGCAAAAATTTCTCTCTGGATGCGGGGTCAACGCATTTAGGGGaacagc     >  2:495324/1‑232 (MQ=255)                                                                                                                                                                                                                                          aGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGAc                                                                 >  2:477289/1‑176 (MQ=255)                                                                                                                                                                                                                                          aGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGAAGAAAATCGCAAATGCACATCAAGCTAAAATTCCTTCCATGATGCTGGATCATCGCAATTAGGGTAATCTTCaa  >  2:444345/1‑238 (MQ=255)                                                                                                                                                                                                                                          aGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGGAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTTGTGCAGAAAATCGACAATGAACATCCAGCTAATATTGCTTTCATGTTACTGGAAAAACAAAATTAGGGAAATCGGCaa  >  2:710071/1‑238 (MQ=255)                                                                                                                                                                                                                                          aGAATACCGCTACCTTCCAGCGCAGT‑AGTA‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGAc                                                                 <  1:477289/176‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            AGCAATCCAGCAAATACCGAAAGGGGACGGGCGAAGAATATCATCCATGCGGACAAAATGAGCGCCGGAATGGCAATAGGCAGCAGATCGCTTGGGTTAACCAACAGCCCCAGCACCAGGAACATGGCGATTTGCGCCA‑GCCAGGCGAGGCCGTCGAAATTTTGCAGGATGCCGT‑AGCGGTTGCGAATAGGGCGA‑TT‑ACCCAGCAGAAAACCGCACAGATACA‑CCGCCAGAATACCGCTACCTTCCAGCGCAGT‑AGTT‑AATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGGCGGGCAGGGCGATGCGGTTGATCATTTGCAGTAGTAAATAACCGCCGCCAAGACCTATTACAATGCCCAGACCAAATTGTTGCAGAATATCGACAATGAACATCCAGCTAATATTGCTCTCATGATGCTGGATCATCGCAATTAGGGTAATCGTCAG  >  NC_000913/1241104‑1241570

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