Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,909	C→T	G249G (GGC→GGT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,909	0	C	T	88.3%	136.8 / 8.0	63	G249G (GGC→GGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base T (25/28);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

CATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTC  >  NC_000913/12666‑13151                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   catcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGAc                                                                                                                                                                                                                                              <  2:390946/250‑1 (MQ=255)    caACGAACCGACCGCAGCTGCGCTGGGTTACGGTCTTGACAAAGGTACTGGCAACCGTACTATCGCGGTTGATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTGCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTg                                                                                                                                                                                                                                           <  1:84621/250‑1 (MQ=255)     aaCGACCCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCCACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATAGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAACCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTgcg                                                                                                                                                                                                                                         >  2:206573/1‑251 (MQ=255)        gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCACCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGCGAAG‑ACTTCGACAGCCGTCCGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCa                                                                                                                                                                                                                                       >  2:443930/1‑250 (MQ=255)                 gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCt                                                                                                                                                                                                                             <  1:612287/251‑1 (MQ=255)                 gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGGATTGACCTGCGCAACGATCCGCt                                                                                                                                                                                                                             <  2:102943/251‑1 (MQ=255)                              cTTACGGTCTCGACAAAGGTACTGGCAACCGTACTATCGCGGTTTGTGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGGTACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGCTGAAGAATTCAAGAAAGCTCAGGGGATTGACCTGCGCAACGCTCCGCTGGCAATGCAgcgc                                                                                                                                                                                                                >  2:274912/1‑251 (MQ=255)                                     tCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGCCGAAGTTGACGGCGAACAAACCTTCGAAGTTCTGGCCACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGTATTGACCTGAGCAACGATCCGCTGGCAATGCAGCGCCTgca                                                                                                                                                                                                           >  2:543454/1‑247 (MQ=255)                                         gACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGc                                                                                                                                                                                                     <  2:567615/251‑1 (MQ=255)                                             aaGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCa                                                                                                                                                                                                 <  2:557460/251‑1 (MQ=255)                                                  aCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGaaaa                                                                                                                                                                                            <  2:625801/251‑1 (MQ=255)                                                     ggCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGaaaaa                                                                                                                                                                                           <  1:456812/249‑1 (MQ=255)                                                             tACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCg                                                                                                                                                                                 >  1:62058/1‑251 (MQ=255)                                                                tATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGGATTGACCAGCGCAACGATCCGCTGGCAATGCAGCGCCTGCAAGAAGCGGCAGAAAAAGCGAAAATCGAAc                                                                                                                                                                              >  2:137802/1‑251 (MQ=255)                                                                     cgGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCt                                                                                                                                                                         >  2:6800/1‑251 (MQ=255)                                                                           aTGACCTGGGTGGTGGTACTTTCGATAGTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCt                                                                                                                                                                   <  2:249344/251‑1 (MQ=255)                                                                               ccTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGAAAAAACCTTCGAAGATCTGGCAACCAACGGTGATACCCACCTGGGAGGTGAAG‑ACTTCCACAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGACGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                                                                               >  1:275980/1‑251 (MQ=255)                                                                               ccTGGGTGGTGGTACTTTCGATATTTCTATTAGCGAAATCGACGAAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                                                                               <  2:281324/251‑1 (MQ=255)                                                                                 ggggTGGGGGTTCTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACTGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagcag                                                                                                                                                             <  2:336696/250‑1 (MQ=255)                                                                                   ggtggtggtACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGAc                                                                                                                                                           >  1:257232/1‑251 (MQ=255)                                                                                           tACTTTCGATATTAATATTATCGAAATCGACGAAGTTGCCGGAGACACAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGAGGTGCAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGtt                                                                                                                                                    >  2:419170/1‑250 (MQ=255)                                                                                            aCTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTa                                                                                                                                                   <  2:490508/250‑1 (MQ=255)                                                                                             atttCGATATTTCTATTATCGAAATCGACGAAGTAGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAAc                                                                                                                                                 <  1:294944/250‑1 (MQ=255)                                                                                                          tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAAt                                                                                                                                                                                                                       <  1:500646/168‑1 (MQ=255)                                                                                                          tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAAt                                                                                                                                                                                                                       >  2:500646/1‑168 (MQ=255)                                                                                                            ttatCGACATCCGCGAAGTTGACGACGCAAATACCTTCGAAGTTCTGAAAACCAAAGATGATACCCAGATGGGAGGTAACG‑CCTTCGACAGCCGTCTGCTCAACTCCATGGTTGAAGAATTCAAGAAAGCTCAGGGCATTGACCTGCGCAACGATCCGCTGGCACTGCAGCGCCTGAAAGTAGCGGCAGCAAAAACGAAAATCGACCTGTCTTCCGCTCAGCAGCCCGACGTTCACCTTCCGAACATCACt                                                                                                                                  >  2:309762/1‑251 (MQ=255)                                                                                                            ttatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                                                                                  >  1:484449/1‑251 (MQ=255)                                                                                                                cGAAATCGCCGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGACGTACATCACTGc                                                                                                                                >  2:394471/1‑249 (MQ=255)                                                                                                                   aaTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGatg                                                                                                                           >  2:286925/1‑249 (MQ=255)                                                                                                                    aTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTGTATTTACCTGCGCAACGATCCGCTGTCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGatgc                                                                                                                          <  1:591604/251‑4 (MQ=255)                                                                                                                       gacgaAGTTGACGGCGAAAAAACCTTCGACGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGAc                                                                                                                       >  2:185390/1‑251 (MQ=255)                                                                                                                         cgaATTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTTCCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTACCCTGCCGTACATCGCTGCAGATGCGACCg                                                                                                                     <  1:482719/251‑1 (MQ=255)                                                                                                                              ttGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGATGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGACTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCg                                                                                                                >  1:375631/1‑251 (MQ=255)                                                                                                                                   ggCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGACGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCGCTGCAGATGCGACCGGTCCGAAaca                                                                                                           <  2:349468/251‑1 (MQ=255)                                                                                                                                    gCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTACCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAacac                                                                                                          <  2:162937/251‑1 (MQ=255)                                                                                                                                                 cGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAACGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGCTGCGACCGGTCCGAAACACATGAACATCaa                                                                                               >  1:359086/1‑249 (MQ=255)                                                                                                                                                     gTTCTGTCAACCAACGGTGAGACCCACTTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATGGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACGGTCTTCCGCGCAGCAGACCGACGTTAACCTGCCGTACATCACTGGAGATGCGACCGGTCCGAAACACAGGAACATCAAAGttac                                                                                         <  2:62050/251‑4 (MQ=255)                                                                                                                                                     gTTCTGGCAACCAACGGTGATACCCCCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGGTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGAc                                                                                         <  2:53247/251‑1 (MQ=255)                                                                                                                                                          ggCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCg                                                                                      <  1:6800/249‑1 (MQ=255)                                                                                                                                                                 aacGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACt                                                                             >  1:79777/1‑251 (MQ=255)                                                                                                                                                                       gatccccACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCGGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTGACTCGTGCGAAACTGGAAAg                                                                       <  2:359086/247‑1 (MQ=255)                                                                                                                                                                       aaaaCCCACCTGGGTGGTGAAG‑ACTTCGAATGCCGTCTTAGCAACTTCCTGGGTGAAAAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAATTGACTCGTGCGAAACTGGAAAg                                                                       <  2:80329/248‑1 (MQ=255)                                                                                                                                                                           cccACCTGGGGGGAGAAGCACT‑CGCCACCAGTCTGATCACCTACCTGGTTGAAGAATTCAGGAACGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGACAAAGCGAAAATCGAACTGTCTTCCGCTAAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACCTGAACATCAAAGTGCCTCGTGCGAAACTGGGAAGCCTg                                                                   >  2:621188/1‑251 (MQ=255)                                                                                                                                                                                cTGGGTGGTGAAG‑ACTTCGCCAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATAAGGGTCTTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTg                                                               >  2:62082/1‑250 (MQ=255)                                                                                                                                                                                   ggtggTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGa                                                           <  1:41494/248‑1 (MQ=255)                                                                                                                                                                                         gAAG‑ACTTCGACAGCCGTCTGATTAACTACCTGGTTGAAGAATTCAATAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTTCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTgg                                                      <  1:393122/250‑1 (MQ=255)                                                                                                                                                                                           aG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTaa                                                   <  1:394471/251‑1 (MQ=255)                                                                                                                                                                                                 tCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAAcc                                                <  1:352166/249‑1 (MQ=255)                                                                                                                                                                                                    aCAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGCAGATCTGGTAAACCGTTcc                                           >  1:605826/1‑251 (MQ=255)                                                                                                                                                                                                      aGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCAt                                         >  2:5916/1‑251 (MQ=255)                                                                                                                                                                                                          gTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                       >  1:65167/1‑249 (MQ=255)                                                                                                                                                                                                            cTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTG‑CAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGcc                                   >  2:418949/1‑250 (MQ=255)                                                                                                                                                                                                              gATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTG‑CAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGcc                                   <  1:418949/248‑1 (MQ=255)                                                                                                                                                                                                                 cAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTg                               >  1:171869/1‑250 (MQ=255)                                                                                                                                                                                                                        cTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGAAACGATCCGCTGGAAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGCCCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGc                       >  2:412961/1‑251 (MQ=255)                                                                                                                                                                                                                          tgTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGAGCCGCTGGCAATGCAGGGCCTGAAAGAAGCGTCCGAAAAATCGAAAATCGAACTGTGTTCCGCTCAGCAGTCCGATGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACAGCAAAGTGACGCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCa                      <  1:581646/249‑1 (MQ=255)                                                                                                                                                                                                                                 gaaTTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAgg               >  1:35673/1‑250 (MQ=255)                                                                                                                                                                                                                                  aatgcAAGAAAGATCAGGGTATTGCCCTGCGCAACGATCCGCTGGCAATGCAGCTCCTGAAAGAAGCGGCAGAAAAAGCGAAACTCTAACTGTCTTCCGCTCATCAGACCGACGTTAACCTGCCGTACATCACTGCAGAGGCGACCGGTCCGAAACACATGAACAGCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGTGCCGCTGAAAGTTTCCCTGCAggcc             <  2:222136/247‑3 (MQ=255)                                                                                                                                                                                                                                    ttCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGc           >  1:366800/1‑251 (MQ=255)                                                                                                                                                                                                                                        agaaagaTCCGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGc       >  1:517080/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaagaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCt     <  1:116228/251‑1 (MQ=255)                                                                                                                                                                                                                                           aagaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGGAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTg    >  2:216500/1‑251 (MQ=255)                                                                                                                                                                                                                                             gaTCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGACTGTc  >  2:182391/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTC  >  NC_000913/12666‑13151

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