Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,283	T→G	19.4%	V102G (GTG→GGG)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,283	0	T	G	19.4%	176.9 / 3.7	74	V102G (GTG→GGG)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand):  ref base T (19/37);  new base G (16/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.93e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAGCG  >  NC_000913/1406054‑1406520                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGATGCGCAGTGGTAAATAT‑CCGCAACAAGCGAGCCTGAACTGTTTTG‑GCCAGCCCGGGGAAAAGGGCCCCAAAAAGGGGTTcgattt                                                                                                                                                                                                                                        >  2:1041113/1‑234 (MQ=255)      gCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCCGTGGGAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTTGGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTttttttt                                                                                                                                                                                                                                   >  1:95698/1‑238 (MQ=255)        aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCCGGATCAGCTGCGCAGGGGTAAATAT‑CCGCAACAGGCGGGCCTGAA‑TCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCCAAAAGGGGTTCAGTTTTTTTc                                                                                                                                                                                                                                   >  2:1059886/1‑237 (MQ=255)         tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGCAGGGGTGCAACATGGGGTGCTGGGTAAGCTGCGCAGGGGGAAATAC‑ACGCAACAGGCGAGCCCGAC‑ACTTTTGCCCCCGCCGGGGGGCAAATGCGCCAAAAAGGGGTGCAGTTGTAtttaca                                                                                                                                                                                                                                >  2:984002/1‑235 (MQ=255)         tCCAACGCGTAACCAACCTGCGCCGCGGCGTATCGACACGCTGCAGGTTGATAAGTTCCTCCCCGCCGGATTTCCCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaa                                                                                                                                                                                                                                 <  1:503092/238‑1 (MQ=255)           cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGGAGGGTTGCAGCATGGGGGGCTGGGTAAACGCCGCCGTGGGAAATAT‑CCGCAACAGGGGAGCCTGAA‑CCTTTTGAGCCAACCGGTGGAAGAGTGCGCCAAAAAGGGGtttgtgtttattacaaaa                                                                                                                                                                                                                              >  2:918961/1‑223 (MQ=255)              cgggTAACCAGCGGGCGCCGCAGTGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                           <  1:323391/236‑1 (MQ=255)                             cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAAGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGGGGGAAATAT‑CCGCAAAAGGCGGGCCTGGA‑TCTTTTTGGCCCGCCGGGGGAAGAGGGCCGCAAAACGGGGGTCGGTTTTTTTtaaaaaaacccgggggggggggtg                                                                                                                                                                                                            >  2:456553/1‑221 (MQ=255)                               ccgatgtgTACCGACACGCTGCCACTTGATAATTTCCACACCACCGGATTTATCGACATCATCCCAATAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCg                                                                                                                                                                                                          <  1:734359/232‑1 (MQ=255)                                      gtttCGACCCCCTGCAGCTTGATATTTTCCTCACCCCCGGATTTCCCGCCATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACAGGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                                                                                                                                                                                   <  2:261257/236‑1 (MQ=255)                                      gtttCGACACGCTGCCGCTTGATAATTTCCTCCCCACCGGATTTCTCGAAATCCACCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCTGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                                                                                                                                                                                   <  2:412634/236‑1 (MQ=255)                                      gtatcgatcccgtgcacgttGAAAATTTCCACCCCCCCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGGTGCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAAGAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                                                                                                                                                                                   <  1:361929/221‑1 (MQ=255)                                      ctttccacacCCTTCAACTTGACACTGTCCTCACCCCCGGATTTCCCGGCTTCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTTCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                                                                                                                                                                                   <  1:472530/233‑1 (MQ=255)                                       taccGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                           <  1:70691/211‑1 (MQ=255)                                       tATCGACCCCCTGACGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                                                                                                                                                                                   <  1:394015/238‑1 (MQ=255)                                          tgACACGCTTCAGCTTGAATATTTCCTTCCCCCCGGATTTCTCGACATTCTCCCACTAAGGCAGGCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCGGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGa                                                                                                                                                                                               <  2:247418/238‑1 (MQ=255)                                             cacGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGGTAAGCTGCGCAGTGGGAAATAT‑CCGCAACAGGCGGGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAAGGGGCCGCAAAATGGTGTTCAGTTTTTTTTACCAAGCCCTGGGGGAGGGttttgggaaaaggcgggttt                                                                                                                                                                                            >  2:351887/1‑222 (MQ=255)                                                ggtGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTc                                                                                                                                                                                         <  1:645576/237‑1 (MQ=255)                                                gCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGGGTTCAGTTTTTTTCAACAAACCCTTGCGGATGGGTTTGGTAACGTGCTGgttatt                                                                                                                                                                                          >  1:412634/1‑238 (MQ=255)                                                cgcgCAGCTTTATAATTTCCTCCCCCCCGGATTTCCCGACATCACCCCACTAAGTCCGCCGGGGGAGTTTCGGCGGGAAGGGTTGCAACAGGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTc                                                                                                                                                                                         <  2:681122/236‑1 (MQ=255)                                                 cTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGGGCCTGAA‑TCTTTTGCCCCAACCGGTGGAAGAGGGCCCCAAAATGGGTGTCAGTTTTTTTCAACAAACCCCGGCGGAgggggggtcctttttacaacttgacc                                                                                                                                                                                         >  2:70691/1‑215 (MQ=255)                                                 cTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCa                                                                                                                                                                                        <  1:351887/239‑1 (MQ=255)                                                 cTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCCGTTTTTTTCAACAAGCCCTGGCGGGTGGTTTGCGGaaacggctgttttt                                                                                                                                                                                           >  2:285028/1‑224 (MQ=255)                                                      gCTTTATTATTTCCCCACCCCCGGTTTTTCCGACCTCATCCCACTAGTGCAGCCGCGGGAGTGTCGGCGGGAAGGGTTGCAACAGGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTa                                                                                                                                                                                   <  2:969375/239‑1 (MQ=255)                                                        ttGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTa                                                                                                                                                                                   <  1:1057814/237‑1 (MQ=255)                                                           aataTTTCCTCACCCCCGGATGTGTCGACCGCACCCCACTATGTCAGACGCGGGAGTTTCGGCGTGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCTCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAggg                                                                                                                                                                              <  2:652724/236‑1 (MQ=255)                                                                     caccccCGGATTTCGCGACTGTATCCCACTAAGGCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑a                                                                                                                                                                   <  1:530876/239‑1 (MQ=255)                                                                          ccGGATTTCCCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCg                                                                                                                                                              <  1:34844/239‑1 (MQ=255)                                                                            ggATTTCTCGACATCATCCCACTAAGGCAGCCGCGGGAGTTTCGGGGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCa                                                                                                                                                            <  1:576168/239‑1 (MQ=255)                                                                               tttccggACATCATCCCCGTAAGCCAGCCGGTGGAGTTTCGGCGGGAAGGGTTGCAACAGGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATg                                                                                                                                                          <  2:95698/232‑1 (MQ=255)                                                                                     gACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAAGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCCCAAAATGGTGTTCAGTTTTATTCAACAAACCCTGGCGGATGGTTTGCGGAAAGTGCGTGTTATTCATGGGAAAGGGCGGGGGGAT‑AAAACGCAAGGCAacat                                                                                                                                                   >  2:1010415/1‑239 (MQ=255)                                                                                            tcCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCCACAGGCGGGCCTGAA‑TCTTTTGCGCCAGCCGGGGGAAGGGGGCCGCAAAATGGGGTTCAGGTTTATTCAACAAACCCTGGCGGATGGTTTGCGTAACGTGCTGGTTATTTATGGGAAAGGGGGGGGGGGTAAAATCG‑ATGCCAAAAATGTcccc                                                                                                                                            >  2:810049/1‑237 (MQ=255)                                                                                              cccttagtccAGCCGTTGGGGTTTCGGGGGGAAGGGTTGCAACAGGGGGTGCGGGATATGCGCCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAg                                                                                                                                          <  2:269205/230‑1 (MQ=255)                                                                                                 cTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTTCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTTGCGGATGGGTTTCGTAACGTGCTGATTATTCATTGTAAAGGGCGGGCTGAT‑AAATAGCATGCCACATTTTTCCGCAGcca                                                                                                                                       >  1:466036/1‑237 (MQ=255)                                                                                                 cTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGGTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGGGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCcgccgcat                                                                                                                                       >  1:902279/1‑237 (MQ=255)                                                                                                         gCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGgggc                                                                                                                               >  1:788166/1‑239 (MQ=255)                                                                                                              tGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGGGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGGAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGGgctg                                                                                                                           >  1:813676/1‑238 (MQ=255)                                                                                                                    ttCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGGAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGGTTGCGTAACGTGCTGATTATTCATGGTAAAAGGCGGGATGCT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGTGGGCTGGCTGGcc                                                                                                                    >  1:1012964/1‑239 (MQ=255)                                                                                                                          gggAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAAtt                                                                                                               <  1:174297/238‑1 (MQ=255)                                                                                                                              aGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtgatg                                                                                                          <  1:285028/239‑1 (MQ=255)                                                                                                                                     cAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGTTGGCTGACCGAATTTGATGATGTTCa                                                                                                   >  1:571026/1‑239 (MQ=255)                                                                                                                                      aaCATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                                  <  1:795899/239‑1 (MQ=255)                                                                                                                                      aaCAGGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                                  <  1:456553/239‑1 (MQ=255)                                                                                                                                       aCATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                                  <  1:568573/238‑1 (MQ=255)                                                                                                                                          tGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGCTGGTTTGCGTAACGTGCTGATTTTTCATGGTAAAGGGCGGGATGTT‑AAATCGCAGGCCAATTTTTTCCGCAGCTATGTGGCGGGCTGGCTGACCGAATTTGATAAAGTTAAggaaa                                                                                              >  2:319624/1‑236 (MQ=255)                                                                                                                                           ggggTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTTACCGAATTTGATGATGTTCAGGCata                                                                                             >  1:737769/1‑239 (MQ=255)                                                                                                                                           ggggTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGAGGTTCGGGCatt                                                                                             >  1:841147/1‑238 (MQ=255)                                                                                                                                               tgtttgATAAGCGGC‑CAATGGGAAAAATCCCGCAAGAGGGGAGCTTGAA‑TCTTTGGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTg                                                                                          <  2:841147/233‑1 (MQ=255)                                                                                                                                                   ggtttacgttcgcAGTGGTAAATTC‑CCGAAACGGGCGAGCCTGAA‑TTTTTTGCCCCAGGCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACcg                                                                                     <  2:902279/229‑1 (MQ=255)                                                                                                                                                   ggaataggccgcgggTGGTAAATTC‑CCCGAAAAGGGGGGCCTGAA‑TTTTTTGCGCCAGGCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACcg                                                                                     <  2:1012964/225‑1 (MQ=255)                                                                                                                                                   ggATAAGCTGCGGAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACcg                                                                                     <  1:1041113/239‑1 (MQ=255)                                                                                                                                                        aGCTGCGGCGTTGTAAATTT‑CCGGAACGGGAGGGCTTGAA‑TTTTTGTCGCCCGCCGGTGGAAGGGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTa                                                                               <  2:262320/239‑1 (MQ=255)                                                                                                                                                                  gggTAAATTT‑CCGAAACAGGCGAGCCTGAA‑TTTTTTGAGCCAGCCGGTGGAAGGGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTa                                                                               <  1:486183/228‑1 (MQ=255)                                                                                                                                                                       aaTAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGCCCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCTCATCATgccgct                                                                >  1:92174/1‑237 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCg                                                                                                                                                              >  1:751694/1‑141 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCg                                                                                                                                                              <  2:751694/141‑1 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCa                                                                                                                                                            >  2:372660/1‑143 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCa                                                                                                                                                            <  1:372660/143‑1 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTTTCCGCAGCTATGTGGAGCGTTGGCTGCCCGACTTTGATGATGTTCAGGCATATTG‑ACCGCCGCTACCGCTCCATggcggccacgg                                                           >  1:984180/1‑234 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CCCTACCGCATCATGGCGGCacggg                                                           >  1:922354/1‑235 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGAGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCaactg                                                           >  1:409501/1‑235 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCCTGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTTAGGCATATTGCACCG‑CGCTACCGAATCATAGCGGGAgggg                                                           >  2:202568/1‑239 (MQ=255)                                                                                                                                                                             ccGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGAGGATGTTAAGGCATCTTGCACCG‑CGCCACCGCATCATGGCGGCAACgg                                                           >  1:850261/1‑239 (MQ=255)                                                                                                                                                                              cccAACAGGCGGGCCGAAT‑TTTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCggg                                                          <  2:571026/237‑1 (MQ=255)                                                                                                                                                                                   cgggggtgctttca‑TCTTTGGCACCCGCCGGTGGAGAAGTGCGCCAAAATGGTGTTCAGTTTTATTCAACCAGCCCTGGCGGAGGGTTGGCGTAACTGGCTGAGTATTCATGGTAAAGGGCGGGATGAT‑AAAGCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCgtg                                                     <  2:200846/226‑1 (MQ=255)                                                                                                                                                                                     ggccggctttaa‑TTTTTTTCCCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCGTTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTTATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCgcggtgcaat                                                   <  2:466036/229‑7 (MQ=255)                                                                                                                                                                                     ggCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCTTCATGGCGCCAGCGGGGCGGGtt                                                   >  1:19609/1‑239 (MQ=255)                                                                                                                                                                                           ccTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGAGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGACGGGCGGGTTAcagc                                              >  1:127556/1‑235 (MQ=255)                                                                                                                                                                                             tGAA‑GCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCa                                           <  1:918961/239‑1 (MQ=255)                                                                                                                                                                                                                           tGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGGACTACGTAAAACGGCGCAGGCGAAGCAAGa              >  1:554565/1‑239 (MQ=255)                                                                                                                                                                                                                             ccGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGGGTTACGTGCCACACCGCAAAACGGCGCAGGCGAAACaaaaaa            >  2:1077664/1‑239 (MQ=255)                                                                                                                                                                                                                             ccGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCCCATCATGCGCGCAGCCGGGCGGGTTACTGCGCCGTACGGAAAAcacccaagcggaacaaagaaa            >  2:753946/1‑219 (MQ=255)                                                                                                                                                                                                                             ccGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGAGGATGTTCAGGCATATTGCACCG‑CGCTCCCGCCTCATGGCGGCAGCGGGGCGGGTTACGCCGCACTACGTAAAAAGGCGCAGGGGagaagaaaaa            >  2:484726/1‑230 (MQ=255)                                                                                                                                                                                                                                  aaaTGGGTTTTATTTTTATTTAACAAGCCCGGGGGGGTGGTTTGGCTAACGTGCTGATTATTCAAGGTAAAGGGCGGGATGAA‑AAATCGCATGCCAATATTGTGCGCAGTTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTgg       <  2:737769/239‑1 (MQ=255)                                                                                                                                                                                                                                       gtgtTCAGTTTTATTCAAAAAGCCCTGGCGGATGGTTTGCGTAACGTTCTGATTATTCATGGTAAAGGGCGGTATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAgcg  <  1:984002/239‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                CACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATAT‑CCGCAACAGGCGAGCCTGAA‑TCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGAT‑AAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCG‑CGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAGCG  >  NC_000913/1406054‑1406520

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