Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,237	C→T	R25R (CGC→CGT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,237	0	C	T	96.4%	169.3 / ‑2.4	56	R25R (CGC→CGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (36/18);  total (37/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTA  >  NC_000913/12006‑12478                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   tctcCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAACAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCGCTCCTCGTGTGCTGGAGAACGCCGa                                                                                                                                                                                                                                   >  2:464360/1‑249 (MQ=255)    cccccTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCg                                                                                                                                                                                                                              >  1:54722/1‑251 (MQ=255)    cccccTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCg                                                                                                                                                                                                                              >  1:277892/1‑251 (MQ=255)        cTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCg                                                                                                                                                                                                                          >  1:77293/1‑251 (MQ=255)                    ggttACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTc                                                                                                                                                                                                              >  1:251490/3‑251 (MQ=255)                          gACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCAt                                                                                                                                                                                                        <  1:501075/251‑1 (MQ=255)                                      tagtAAACCGCAGTGAGTGAGTCTGAAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCAGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGAACAACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCCCGCCTTCTATCATTGCCTATAccc                                                                                                                                                                                             >  1:173801/1‑250 (MQ=255)                                                 agtgagtgAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGaa                                                                                                                                                                                  >  1:1286/1‑250 (MQ=255)                                                  gtgagtgaGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAAc                                                                                                                                                                                <  1:139079/251‑1 (MQ=255)                                                        gagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAActctggt                                                                                                                                                                          <  1:464360/251‑4 (MQ=255)                                                         agtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTCCCTATACCCAGGATGGTGAAACTCTGGtt                                                                                                                                                                         <  1:533352/251‑1 (MQ=255)                                                              gCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAAGTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCGCGTGTGCTGTAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCa                                                                                                                                                                   <  1:83223/251‑1 (MQ=255)                                                               cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTTGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATTCCCAGGATGGGGAAACTCTAGTTGG‑ttag                                                                                                                                                                  <  1:192684/251‑4 (MQ=255)                                                                 aaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGcc                                                                                                                                                                >  1:652953/1‑251 (MQ=255)                                                                      aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCAAATGATAACCGAATATCTAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACATGGGTACTACCAACTCTTGTGTAGCGAATATGGACGGCACCACTCCTCGTGTGCTGGAGAACGCCGCAGTCGATCGCTCCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCT‑GTTGGGTCAACCGGCt                                                                                                                                                            >  2:485011/1‑250 (MQ=255)                                                                       tGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGACACTCTGGTTGG‑TCAGCCGGCTa                                                                                                                                                           >  1:129902/1‑250 (MQ=255)                                                                                   gTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGGGAAGACGTTTAAATGAATAAAATCATTGGTAAAGAACTAGGGACTACCGACTCTTGTGTAGCGATTATGGATGGCACAACTCAGCGTGTGCTGGAGCACACCGAAGGCGATAGCACCACACCTTCTATCATTGCCTATAACCAAGATAGTGAAACTCTGGTTGG‑TCAGCAGGATAAACGtaagacagg                                                                                                                                              >  2:150810/1‑243 (MQ=255)                                                                                    ttGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTg                                                                                                                                             >  2:321035/1‑251 (MQ=255)                                                                                    ttGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTATGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTg                                                                                                                                             >  1:451867/1‑251 (MQ=255)                                                                                           cAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAAcc                                                                                                                                      >  1:389693/1‑251 (MQ=255)                                                                                                gTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCa                                                                                                                                  >  1:430425/1‑250 (MQ=255)                                                                                                gTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATACAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCa                                                                                                                                  >  1:537419/1‑250 (MQ=255)                                                                                                   tCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTTGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCaaaa                                                                                                                               <  2:251490/250‑1 (MQ=255)                                                                                                    cGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAac                                                                                                                              >  1:334932/1‑250 (MQ=255)                                                                                                       cccTATTAAAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCATCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAacaccc                                                                                                                          >  1:307285/1‑249 (MQ=255)                                                                                                                    tCACAACCAAATGATGACCGAATATATAGTGGAGACGTTTCGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGCACGCCGCAGGCGATCGCACCACGACTTCTATCATTGCCTCTACCCAGGAGGGTGAAACTCTGGTTAG‑TCAGCCGGCTAAACGTAAGCCCGTGCCGAACCCGAAAACCACCATGTTTGCGTTTa                                                                                                              >  2:326692/1‑250 (MQ=255)                                                                                                                      acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAAc                                                                                                           >  1:102822/1‑251 (MQ=255)                                                                                                                                gatgaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCaa                                                                                                                                 >  1:232268/1‑219 (MQ=255)                                                                                                                                gatgaCCGAATATATAGTGGAGACGGTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑GCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCaa                                                                                                                                 <  2:232268/219‑1 (MQ=255)                                                                                                                                    aCCGAATATATAGTGGAGACGTTTAGAGGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGAAGGCACCACTCC‑CGCGTGCTGGAGAACGCCGAAGGCGATAGCACCACGCCTTCTCTCATTGCCTATACCCAGGATGGTGAAAATCTAGTTGG‑TCAGCCGGCTAAACGTCCGGCAGTGACGCACCCGCACAAAACTCTGTTTGCGATTAAACGCCTGATTGGTcgcc                                                                                            >  2:25964/1‑251 (MQ=255)                                                                                                                                               aGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTCTACCCAGGATGGTGAAACTCTGGTTGG‑TCAACCGGCTAAACGTCAGGCAGTGACGAACACGCAAACAACCCAGTTTGCGTTTAAAAGCCTGATTGGCCGCCGCTTCCAggg                                                                                  >  2:359242/1‑250 (MQ=255)                                                                                                                                                     gaCGTTTAGATGGGCAAAATAATTGGTATCGACCTGGGTACTACCCACTATTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGCTGGTAAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGGGACGAACACGCAAAACACCCTGTTTGCGATTAAAAGCCTGAGTggccgacaattacaggcaga                                                                               >  2:665322/1‑230 (MQ=255)                                                                                                                                                        gTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGt                                                                         >  1:211786/1‑251 (MQ=255)                                                                                                                                                                  gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                               <  1:154278/251‑1 (MQ=255)                                                                                                                                                                      aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGttt                                                           >  1:635786/1‑251 (MQ=255)                                                                                                                                                                            tGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatc                                                      <  1:160228/250‑1 (MQ=255)                                                                                                                                                                             ggTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatca                                                     <  1:503846/250‑1 (MQ=255)                                                                                                                                                                              gTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATg                                                   >  1:271009/1‑251 (MQ=255)                                                                                                                                                                                     ccTGGGTACTAACAACTCTTGTGTAGCGATTATGGATGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTAAGGCAGTGACGAACACGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCAAGGACGAAGAAGAACAGCGTGAAGGTTCCATCATGCCGTTc                                             >  1:450239/1‑250 (MQ=255)                                                                                                                                                                                      cTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGCCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCa                                            >  1:314579/1‑250 (MQ=255)                                                                                                                                                                                        gggTACTACAAACTCTTGTGTAGCGTGTATTGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCTATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCaaaa                                         <  2:113325/251‑1 (MQ=255)                                                                                                                                                                                               acCAACTCTGGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCGGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgc                                  <  1:665322/251‑1 (MQ=255)                                                                                                                                                                                                 cAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctg                                >  1:497239/1‑251 (MQ=255)                                                                                                                                                                                                  aaCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctgc                               >  1:164416/1‑251 (MQ=255)                                                                                                                                                                                                      cctgtgTAGAGATTATGGATGGCACCACTCCTCTTGTGCTGGAGAACTCCGAAGTCGAGCGCACCACGCCTGCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGACGGCTAACCGTCATGCAGTGACGAACCTGCAAAACTCCCTGTTTGCGATTAAACGCCTGAGTAGCCGCCGATTCCAGGACGAAGAAGTACAGCGGGATGTTTCCTTCATGCCGTCCAAAATTATTGCTGCTGAt                           <  2:118647/249‑1 (MQ=255)                                                                                                                                                                                                        tgtgTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGa                            <  2:102822/248‑1 (MQ=255)                                                                                                                                                                                                         gtgtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTACGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAAAGCGTGATGTTTCCATCCTGCCGTTCAAAATTATTGCTGCTGATAAc                        >  2:557820/1‑251 (MQ=255)                                                                                                                                                                                                           gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                      >  2:183598/1‑203 (MQ=255)                                                                                                                                                                                                           gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                      <  1:183598/203‑1 (MQ=255)                                                                                                                                                                                                            tAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                     >  2:388563/1‑251 (MQ=255)                                                                                                                                                                                                             aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                     <  1:321035/250‑1 (MQ=255)                                                                                                                                                                                                                gATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCtata                                                                                                                                                                                                >  1:167116/1‑77 (MQ=255)                                                                                                                                                                                                                gATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCtata                                                                                                                                                                                                <  2:167116/77‑1 (MQ=255)                                                                                                                                                                                                                       gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                   >  1:218456/1‑178 (MQ=255)                                                                                                                                                                                                                       gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                   <  2:218456/178‑1 (MQ=255)                                                                                                                                                                                                                                acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa  >  1:176007/1‑250 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   TCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGG‑TCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTA  >  NC_000913/12006‑12478

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