Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,714	C→T	G184G (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,714	0	C	T	91.8%	157.0 / 2.8	65	G184G (GGC→GGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (27/29);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

GTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGG  >  NC_000913/12475‑12958                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      gTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCg                                                                                                                                                                                                                                            >  2:124240/1‑251 (MQ=255)    aaaggcaaataaatggacccgccgCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAATGTACTGGCAACCGTa                                                                                                                                                                                                                                          <  1:502715/232‑1 (MQ=255)       ggCCAGAAAATGGCCCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTAt                                                                                                                                                                                                                                      <  1:17325/251‑1 (MQ=255)        gCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATc                                                                                                                                                                                                                                     >  1:108748/1‑251 (MQ=255)        gCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAACCGTACTATc                                                                                                                                                                                                                                     >  1:274820/1‑251 (MQ=255)          aaaaaaaGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATcgc                                                                                                                                                                                                                                   <  2:491558/248‑1 (MQ=255)             aaaaTGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAACCGTACTATCGAGGt                                                                                                                                                                                                                                >  2:505454/1‑251 (MQ=255)                     aCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGAcc                                                                                                                                                                                                                        >  1:324705/1‑251 (MQ=255)                         ccgcAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑ggg                                                                                                                                                                                                                   <  1:313007/251‑1 (MQ=255)                             aGATCTCTGCTGAAGTGCTGAACAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGCAGTAACACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAATGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑Gggtggg                                                                                                                                                                                                               >  2:183988/1‑250 (MQ=255)                               aactCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑Gggtggtgg                                                                                                                                                                                                             <  1:418961/248‑1 (MQ=255)                                tctCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑Gggtggtggt                                                                                                                                                                                                            >  1:251861/3‑251 (MQ=255)                                   ctgctgAAGTGCTGAAAAAAATGAAGAAACCCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCCCCGTACCGGCATACTTTAACGATGCTCGGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑Gggtggtggt                                                                                                                                                                                                            <  2:108748/248‑1 (MQ=255)                                     gctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAATCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACt                                                                                                                                                                                                         <  1:38890/249‑1 (MQ=255)                                     gctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATcgc                                                                                                                                                                                                                                   >  1:180640/1‑224 (MQ=255)                                     gctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATcgc                                                                                                                                                                                                                                   <  2:180640/224‑1 (MQ=255)                                     gctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCT‑Gggtggtgttctg                                                                                                                                                                                                         <  1:124240/249‑6 (MQ=255)                                      ctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTc                                                                                                                                                                                                      >  1:162467/1‑251 (MQ=255)                                        gAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGCCCT‑GGGTGGTGGTACTTTCg                                                                                                                                                                                                     >  2:26321/1‑250 (MQ=255)                                         acgTGCTGAAACCAATGAAGAAAACCGCTGAAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCCAAGACGCAGGCCGAATCGCTGGTCTGGCATTAAACCGTATCATCAACGCACCGCCCGCAGCTGCTCTGGCTTACGGTCTGTACTACGGTACTGGCAACCGTACTATCGCGGTTTATGCCCTTGGGTG‑TTGTACTTTCGat                                                                                                                                                                                                   >  2:56740/3‑251 (MQ=255)                                                       aTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCg                                                                                                                                                                                      >  1:157497/1‑250 (MQ=255)                                                          aagaaAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCAACGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAAt                                                                                                                                                                                  <  2:14311/251‑1 (MQ=255)                                                              aaaCCGCTGACGATTACCTGGGGGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGCACCGACCGCAGCTGCGCTGGCTTCCGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATcgac                                                                                                                                                                              >  2:235995/1‑251 (MQ=255)                                                                  cGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAAAGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATcgacga                                                                                                                                                                            >  2:505120/1‑249 (MQ=255)                                                                   gcagAAGATTGCCTGGGTGAACCGGTAACCGAAGCTGTTATTACCGTACCGTGATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGGCTGGAAGTAAAACGTATCAGCAACGAACCGACCGCAGCGGCGCTGGGTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATTACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGTCGAAGt                                                                                                                                                                         <  2:82162/248‑1 (MQ=255)                                                                               cTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGAGaaa                                                                                                                                                              >  1:250132/1‑250 (MQ=255)                                                                                     gAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCtt                                                                                                                                                       >  1:173405/1‑251 (MQ=255)                                                                                      aaCCGGTAACTGAAGCTGGTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATACTGCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTc                                                                                                                                                      <  2:497469/251‑1 (MQ=255)                                                                                          ggTAACTGAAGCTGTTATCCCCGTCCCGGCCTACTTTAACGATGCTCAGCGGCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAg                                                                                                                                                  <  2:274820/251‑1 (MQ=255)                                                                                              aCTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCt                                                                                                                                              <  2:324705/251‑1 (MQ=255)                                                                                                          aTCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACgg                                                                                                                                  >  2:464188/1‑251 (MQ=255)                                                                                                            cACCGTAGCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCACCCGTACTATCGAGGTTTATGAACT‑GGGTGGTGGTACTTTCGATATTTCTATCATCGAAATCGACGCAGTTGACGGCGAAAAAACCTTCGAAGTTCAGGCAACCAACGGTg                                                                                                                                >  2:99433/1‑251 (MQ=255)                                                                                                                    cGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGCACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGAGGAGGTACTTTCGATATTTCTATTATCGACATCGACGAAGTTGCAGGCGAAACAACCTTCGAAGTTCGGGCAACAAACgg                                                                                                                                  >  2:35416/1‑241 (MQ=255)                                                                                                                    cGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACgg                                                                                                                                  <  1:35416/241‑1 (MQ=255)                                                                                                                      gcagactttaacgaacctagtcgccAGGCACCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCAACAACGAACCGACCGCAGCGGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTTGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCTAAGTTCTGGCAACCAACGGTGATTCCCACCt                                                                                                                      <  2:200040/227‑1 (MQ=255)                                                                                                                       cATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTg                                                                                                                     <  1:153494/251‑1 (MQ=255)                                                                                                                               aaCGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGgtggtg                                                                                                              >  1:351602/1‑250 (MQ=255)                                                                                                                                   aTGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGAc                                                                                                         <  1:151821/251‑1 (MQ=255)                                                                                                                                      cTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTc                                                                                                      >  1:442567/1‑251 (MQ=255)                                                                                                                                       tCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCACCGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGCGGTGAAGACTTCg                                                                                                     >  2:219083/1‑251 (MQ=255)                                                                                                                                                  cAACCAAAGACGCAGGCCGTATCGCTGGACTGGAAGTAAAACGTATCATCAACGAACCGACCGCAACTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTg                                                                                          <  2:187465/251‑1 (MQ=255)                                                                                                                                                                ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGa                                                                                                    <  1:348355/227‑1 (MQ=255)                                                                                                                                                                ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTGCGa                                                                                                    >  2:348355/1‑227 (MQ=255)                                                                                                                                                                ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTGCTGGCAGCCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGt                                                                            <  1:235995/251‑1 (MQ=255)                                                                                                                                                                         gCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAAtt                                                                   <  1:436727/251‑1 (MQ=255)                                                                                                                                                                          cTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGATTTc                                                                  >  2:444713/1‑251 (MQ=255)                                                                                                                                                                            ggTCTGGAATTAAAACGGTTCATCAACGAACCGACCACAGCTGCGCTGGCTTACGTTCTGGACAAAGGTACTGGCAACCGTACTAGCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTATGATCAACGACCTGGTTGAAGAATTCaa                                                                <  2:466947/251‑1 (MQ=255)                                                                                                                                                                             gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTGGAAGAATTCacg                                                               >  2:410941/1‑249 (MQ=255)                                                                                                                                                                               ttGGAAGTAAAACGTATCATCAACACACCGACCGCAGCTGCGCTGGCTTACGTTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACATTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATTCCCACCTGGGTGGTGAAGACTGCGACAGCCGTCTGATCAACTGCCTGTTTGAAGAATTCaaga                                                              <  2:37624/249‑1 (MQ=255)                                                                                                                                                                                    aGTAAAACGTATCATCAACGAACCGACCGCAGATGCGATGGCGTACGGGCGGGAAAAAGGGACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAGACTGCGACAGACGTCTGATCAACTACCTGGTTGAAGAATTCaagaaagagc                                                        <  2:118273/251‑3 (MQ=255)                                                                                                                                                                                     gTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCa                                                       <  1:79492/251‑1 (MQ=255)                                                                                                                                                                                         aacttATCATCAACGAACCGACCGCAGCTGCGCTTGCTTACTGTCTGGACAAATGTTCTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCCGGGc                                                   <  2:359977/247‑1 (MQ=255)                                                                                                                                                                                          aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGc                                                   >  1:452357/1‑250 (MQ=255)                                                                                                                                                                                          aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGc                                                   >  1:156207/1‑250 (MQ=255)                                                                                                                                                                                           cGTATCAACAACGAACCGTCCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGTCCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCAt                                                 <  1:183988/251‑1 (MQ=255)                                                                                                                                                                                                  tcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGAc                                                                                                         >  2:219370/1‑188 (MQ=255)                                                                                                                                                                                                  tcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGAc                                                                                                         <  1:219370/188‑1 (MQ=255)                                                                                                                                                                                                     aCGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTgcgc                                       >  1:294850/1‑251 (MQ=255)                                                                                                                                                                                                       gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGAGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACAAACGGTGATACACACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCa                                      >  1:50896/1‑250 (MQ=255)                                                                                                                                                                                                          ccgaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGa                                  <  1:464188/251‑1 (MQ=255)                                                                                                                                                                                                                               tACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCt             >  1:499822/1‑251 (MQ=255)                                                                                                                                                                                                                               tACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCA‑GAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGcc              <  1:48052/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAaag        <  1:117177/248‑1 (MQ=255)                                                                                                                                                                                                                                          cAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCgg  >  2:97799/1‑251 (MQ=255)                                                                                                                                                                                                                                          cAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATAAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCgg  >  2:322642/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      GTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCT‑GGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGG  >  NC_000913/12475‑12958

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