Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,014	A→G	P284P (CCA→CCG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,014	0	A	G	94.7%	218.0 / 1.5	75	P284P (CCA→CCG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (43/28);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

TATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAA  >  NC_000913/12771‑13250                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAATGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATc                                                                                                                                                                                                                                        >  1:211237/1‑250 (MQ=255)   ttatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                                                                                                                                                                                     >  1:571416/1‑251 (MQ=255)    taacGAAATCGACGAAGGTGACGGCGAAAAAACCGACGAAGTTCTGGCAACCAACTGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                                                                                                                                                                                     <  2:130886/247‑1 (MQ=255)     atCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGc                                                                                                                                                                                                                                   >  2:286867/1‑251 (MQ=255)                  aaGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCCACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGt                                                                                                                                                                                                                      >  1:639562/1‑251 (MQ=255)                      tGACGGCGAAAAAACCTTCGAAGTGCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCg                                                                                                                                                                                                                   <  2:217711/250‑1 (MQ=255)                        ccctacacgacactctcccgactttGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGAAAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGTAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                                                                                                                                                                <  2:369292/227‑1 (MQ=255)                                 aaaCCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAc                                                                                                                                                                                                       >  1:667448/1‑251 (MQ=255)                                  aaCCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACa                                                                                                                                                                                                      <  1:456975/251‑1 (MQ=255)                                   aCCTTAGAAGTTCTGGCAACCAAAGGTAACACCAAACTGGGTGGTGACGACTCCGACCGCAGTCAGATCAACTACCTGGTTGAAGAATTCAAGAAAGCT‑CAGGGCATTGACCTGCGCAACGATACGCTGGCAATGCAGCGCCTGCAAGAAGAAGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGCCCGCAGTTAACCAGCCGTCAATCACTGCAGATGCGACCGCCACGAAACACATGaaat                                                                                                                                                                                                      >  2:659860/1‑248 (MQ=255)                                         gAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGt                                                                                                                                                                                               >  1:294601/1‑251 (MQ=255)                                                tGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                                                                                                                                                                >  1:369292/1‑227 (MQ=255)                                                    aaccaacGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCg                                                                                                                                                                                     >  1:326817/1‑250 (MQ=255)                                                       caacGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAAc                                                                                                                                                                                 <  1:286867/251‑1 (MQ=255)                                                           ggTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGAT‑AAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCCGAAAAAGCGAAAATCGACCTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGaaccggg                                                                                                                                                                              >  2:306015/1‑245 (MQ=255)                                                            cctaTACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGACCAACTATCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAA‑GCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                                                                                                                                                           <  1:342057/248‑1 (MQ=255)                                                             ggATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                                                                                                                                                           <  1:567715/250‑1 (MQ=255)                                                              gATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGACAGAT‑CAGGGCATTGACCTGCGCAACGATACGCTGGAAAAGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGACAATCGCACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACAAATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                                                                                                                                                           >  2:399771/1‑250 (MQ=255)                                                               aTACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGCAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAA‑GCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                                                                                                                                                           >  2:342057/1‑248 (MQ=255)                                                                 aCCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAACATCGAACTGTCTTCCGCTCAGCAGACAGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGCAACACATAAACATCAAAGAGACTCGTGCGCAACTGGACAGCCt                                                                                                                                                                       >  2:284925/1‑251 (MQ=255)                                                                  cccACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTg                                                                                                                                                                      <  1:49258/251‑1 (MQ=255)                                                                    cACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGt                                                                                                                                                                    >  1:427149/1‑251 (MQ=255)                                                                     cccTGGGTGGGGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAACTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGtt                                                                                                                                                                   <  1:481796/250‑1 (MQ=255)                                                                     aCCTGGGTGGTGCAGACTTCGACAGCCGTCTAATCAACTACCTGGTTGAAGCATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGAACCGCTGGCAATGCAGCGCCTGCAAGAAGCGGCAGAACAAGCGAAAATCGAACTGTCTTCCGCTCCGCAGGCCGCCGTTACCCTTCCGTACATCACTGCCGATGCGACCGGTCCGAAACAAATGACAATCAAAGTGACTCGTGCGAAACAGGAAAGCCTGGtt                                                                                                                                                                   >  2:673158/1‑251 (MQ=255)                                                                      ccTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTg                                                                                                                                                                  <  1:429219/251‑1 (MQ=255)                                                                           gtggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGa                                                                                                                                                                 >  1:214669/3‑247 (MQ=255)                                                                           gtggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGAt                                                                                                                                                             <  1:35365/249‑1 (MQ=255)                                                                             ggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTACGAAACACATTAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTAGAAGATCt                                                                                                                                                           >  2:324918/1‑251 (MQ=255)                                                                                   gACTTCGCCAGCCGTCTGATCAACTACCTGGTTGAAGAATTCACGACCGAT‑CAGGGCATTGACCTGCGCAACGATCAGCTGGCAATGCAGCGCCTGAAAGAAGCGGAAGAAAAAGCGAAAATAGAACTGTCTTCCGCTCAGCAGACCGCCGTTAACCTGCCGTACATCACTGCAGCAGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGAGAAACTGGAAAGCCTGGTTGAAGATCTggaaa                                                                                                                                                      >  2:435838/1‑247 (MQ=255)                                                                                       tCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCg                                                                                                                                                  <  1:399771/250‑1 (MQ=255)                                                                                              ccGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGAGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg                                                                                                                                          >  1:487153/1‑251 (MQ=255)                                                                                                 tCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCGGGTTGAAGATCTGGTAAACCGTTCCATTgcgc                                                                                                                                       >  2:225634/1‑248 (MQ=255)                                                                                                   tGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATCGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCa                                                                                                                                                                                                                                       >  1:85603/1‑153 (MQ=255)                                                                                                   tGATCAACTACCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATCGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCa                                                                                                                                                                                                                                       <  2:85603/153‑1 (MQ=255)                                                                                                          ctcccTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGTCGTTAACCTGCCGTACATCACTGCAGATGCGACAGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGACACTGGAAAGACTGGTGGAAGATCTGGTAAAccataccagtgcaccgctgaagg                                                                                                                              >  2:142895/5‑230 (MQ=255)                                                                                                            accTGGTTGAAGAATTCAAGAAAGATACAGG‑CATGGACCTGCGCAACGATCCGCTGGCCATGCAGCGCCAGAAAGAAGCGGCAGAAAAAGCGAAAAGCGAACTGTCTACCGCTCAGCAGACCGCCGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGtt                                                                                                                            <  2:462373/249‑1 (MQ=255)                                                                                                                        aaTTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGAc                                                                                                                <  1:442594/251‑1 (MQ=255)                                                                                                                         aTTCAAGAAAGAT‑AAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAACCACATGAACATCAAAGTGACTCGAGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTCCCATTGAGCCGCTGAAAGTTGCAcagcagtac                                                                                                                >  1:493258/1‑247 (MQ=255)                                                                                                                           tCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCt                                                                                                             >  1:211645/1‑251 (MQ=255)                                                                                                                                        aGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGt                                                                                                  >  1:495653/1‑250 (MQ=255)                                                                                                                                             aTTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTg                                                                                             >  1:112828/1‑250 (MQ=255)                                                                                                                                             aTTGACCTGCGCAACGATCCGATGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGATCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTAGAAAGCCTGGTTTAAGATCTGGTAAACCGTTCCATTGATCCGCTGAAAGTTGCACTGCAGGAAGCTGGCCttcccggatctgt                                                                                            >  2:86269/1‑239 (MQ=255)                                                                                                                                                gACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATAAAAGTGACTCGTGCGAACCTGGAACGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCGGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGatat                                                                                         >  2:242462/1‑251 (MQ=255)                                                                                                                                                  ccTGCGCAACCATCCACTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAACAGCGAAAATCGAACTGTCTTCAGCACAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAAATCAAAGTGACTCGTGCGAACCTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGACAGTTGCACTGCAGGACGCAGGCCTGTACGTATCTGATATcg                                                                                       >  1:272213/1‑251 (MQ=255)                                                                                                                                                     gcgcAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATcga                                                                                      <  1:27867/249‑1 (MQ=255)                                                                                                                                                     gcgcAACGATCCGCGGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATcgac                                                                                     <  2:571416/250‑1 (MQ=255)                                                                                                                                                      cgcAACGACCCGCTGGCAATGCAGCGCCTGAACGAAGCGGCAGAACAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAACCTGGAAAACCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCCCTGCAGGACGCTGGCCTGTCCGTATCTGATATcgacgc                                                                                   >  2:44759/1‑250 (MQ=255)                                                                                                                                                        caactacccGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTTACTCGTGCGAAACTGGAAAGTCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGacgacg                                                                                 <  1:174271/244‑1 (MQ=255)                                                                                                                                                        cAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGACAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATcgacgccg                                                                                 >  2:83960/1‑248 (MQ=255)                                                                                                                                                           cGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTa                                                                              <  1:284925/251‑1 (MQ=255)                                                                                                                                                                   tGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGtt                                                                      >  1:23761/1‑251 (MQ=255)                                                                                                                                                                     caaagtccgggtcTGAAAGAAGCGGCACAAAAAGCGAAATGCTAACTGTCTTCCGCTCAACAGACCGACGTTAACCTGCCGTACATCACTGCATATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCTTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCTCGGAAAGTTGCACGGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTtgg                                                                    <  2:426652/239‑1 (MQ=255)                                                                                                                                                                      cAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTtggg                                                                   >  1:96310/1‑250 (MQ=255)                                                                                                                                                                            aGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTc                                                               >  2:601551/1‑249 (MQ=255)                                                                                                                                                                                    aaaGAAGCGGCATAAAAAGCGAAAAACGAACTGTCATCGGCTCAGAAGTCCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTAt                                                     <  2:23761/251‑1 (MQ=255)                                                                                                                                                                                      agaagCGGCAGAAAAAGCGAAAATCGCACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGAACATCAATGCAGATGCGTCCGGTCCGAAACACAGGACCATCACAGTGACTCGTGCGAAAATGGAAAGCCGGGGTGAAGATATGGTAAACAGTTCCATTAAGCCGCTGAAATTCGCACTGCAGGACGCGAGCCTGTCTGTATCTGATATCTCCAACGTTATCCTCGTtggtggactgagtcgcgtgc                                                   >  2:7427/1‑237 (MQ=255)                                                                                                                                                                                          gCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAAt                                               <  1:306015/251‑1 (MQ=255)                                                                                                                                                                                           cGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTAGTATGCCAATg                                              >  1:106088/1‑251 (MQ=255)                                                                                                                                                                                                       cGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGaaa                                  <  1:601551/251‑1 (MQ=255)                                                                                                                                                                                                            aTCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGc                             <  1:4341/251‑1 (MQ=255)                                                                                                                                                                                                               gAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTg                           >  1:85554/1‑250 (MQ=255)                                                                                                                                                                                                               gAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAAATGGAAAGCCTGGTTGAAGATCTGGTAA‑CCGTTCCATTGAGCCGCTGCAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCACCGACGTTATCCTCGTTGGTGGTCCGCCTCGGATGCCAAGGGTTCGGAAGAAAGTTGCTg                           >  2:573809/1‑249 (MQ=255)                                                                                                                                                                                                                  cTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCa                                                                                                                         <  1:552430/153‑1 (MQ=255)                                                                                                                                                                                                                  cTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCa                                                                                                                         >  2:552430/1‑153 (MQ=255)                                                                                                                                                                                                                     gcTTCCGCACAGCGGAACGACGTTAACCTGCCGTACATCACTGCGGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCt                     <  2:326817/249‑1 (MQ=255)                                                                                                                                                                                                                      cTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCttt                   >  1:120/1‑251 (MQ=255)                                                                                                                                                                                                                         ccGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTgg                 <  1:659627/250‑1 (MQ=255)                                                                                                                                                                                                                            cTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCagaa                                      >  1:260057/1‑226 (MQ=255)                                                                                                                                                                                                                              cagcagACCGACGTGAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAaga           <  1:312103/251‑1 (MQ=255)                                                                                                                                                                                                                                  agcccGACGTTAACCTGCCGTACAGCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGt                                <  2:551258/223‑1 (MQ=255)                                                                                                                                                                                                                                  agACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGt                                >  1:551258/1‑226 (MQ=255)                                                                                                                                                                                                                                      cGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGt    >  2:544028/1‑250 (MQ=255)                                                                                                                                                                                                                                       gACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTg                                                                                                       >  2:38453/1‑150 (MQ=255)                                                                                                                                                                                                                                       gACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTg                                                                                                       <  1:38453/150‑1 (MQ=255)                                                                                                                                                                                                                                          gTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTaa  >  1:558954/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             TATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGAT‑CAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAA  >  NC_000913/12771‑13250

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