Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,507	T→C	I115I (ATT→ATC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,507	0	T	C	88.9%	116.9 / 1.0	45	I115I (ATT→ATC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (26/14);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGT  >  NC_000913/12260‑12753                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       cacgacccacCCCTTCTAACATTGCCTATACCCGGGATGGTGAACCCCTGGTGGGTCAGCCGGCTAAACGTCAGGCGGTGACGAACCCGCAAACCACCCTGTTTGCGATAAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCAAGAAAATGGCACCGCCGCAGATctct                                                                                                                                                                                                                                                     <  2:106687/245‑1 (MQ=255)      accacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTctgctg                                                                                                                                                                                                                                                 >  1:251452/1‑250 (MQ=255)               tCTATCATTCCCTATCCCCAGGATGGTGAACCTCTAGTGGGTCAGCCGGCAAAACGTCAGGCAGTGACGACCCCGCAAAACCCTCTGTTTGCGATTAACCGCCTGTTTGGCCGCCGTTTCCAGGACGAAGAAGTACAGCGTGATTTTTCCA‑TCATGCCGTTAAAAATTATTTCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGa                                                                                                                                                                                                                                        <  1:234150/250‑1 (MQ=255)                          cTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATga                                                                                                                                                                                                                               >  1:187465/1‑248 (MQ=255)                                       ggTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGCTCTCTGCTGGAGTGCTGAAAAAAATGAAGAAAACCGCTGa                                                                                                                                                                                                                  >  2:117177/1‑248 (MQ=255)                                           aaaCTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTTATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGACGTTAAAGGCCAGAAACTGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAc                                                                                                                                                                                                           >  1:185339/1‑251 (MQ=255)                                                  ggtgggTCAGCCGGCTAACCGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGCCCGCCGCTTCCAGGACGAAAAATTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGt                                                                                                                                                                                                     <  1:218885/246‑1 (MQ=255)                                                  ggTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGAAGTTAAAGGCCAGAAAATGGAACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGACGATTACCTGGGTg                                                                                                                                                                                                    >  2:436727/2‑251 (MQ=255)                                                     tGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                                                                                                                                                 >  2:79492/1‑251 (MQ=255)                                                                             aGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACATGGGTGAACCGGTAACTGAAGCTGTTATCACCg                                                                                                                                                                         >  2:22173/1‑251 (MQ=255)                                                                                      cccGAAAAACCCCCTGTTTGCGATTAACCGCCTTATTGGCCGCTGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCGTTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAATTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAt                                                                                                                                                                <  1:308310/251‑1 (MQ=255)                                                                                       ccGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTCCAGCGAGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGACAAACATGACGAAAACCGCTGAAGATTACCTGGGTAAACCGGTAACTGAAGATGTTATCCCCGTACCgacca                                                                                                                                                                >  2:159682/1‑246 (MQ=255)                                                                                              accccctGTTTGCGATTAAACGCCTGATTGCCCGCCGCTCCCAGGACGAAGAAGTACAGCGTGATTTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAATTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAc                                                                                                                                                        <  1:502702/246‑1 (MQ=255)                                                                                                 ccctGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCGGGACGAAGCAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCACAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGCAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATCCTTTAACGAt                                                                                                                                                     >  2:427324/3‑251 (MQ=255)                                                                                                     gTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTc                                                                                                                                                 >  1:497469/1‑251 (MQ=255)                                                                                                               aaaCGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGc                                                                                                                                       >  1:488404/1‑251 (MQ=255)                                                                                                                    cgtgatgggccgccgcTTCCAGGATAAAGAAGTACAGGGTTATGTTTACA‑AAAGGCGGATAAAAATTATTGGTGTGGAAAAGGGCGACGAAGGGGGCGAAGTAAAAGGCAAGAAAATGGAACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCGGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCa                                                                                                                                  <  2:488404/241‑1 (MQ=255)                                                                                                                     cTGATTGGCCGCCGCTTCCAGGACGAAGAAGTCCAGCGAGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTCCTGAAAAAAATGACGAAACCCGCTGAAAATTACCTGGGGGAACCGGTAACTGAAGCTGTTATCACCATACCGGCATACTTTCACGATGCTCAGCGGCAGGCAACCa                                                                                                                                  >  2:321609/1‑250 (MQ=255)                                                                                                                      tGATTGGCCGCAGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCAA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAAcaaa                                                                                                                                 >  2:443461/1‑247 (MQ=255)                                                                                                                      gtaTTGGCCGCCGCTTCCAGGACGAAAAAGTACGGCTGGATGTTTCAA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGCCGCATGGGTCAAATTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaaa                                                                                                                                <  1:253051/249‑1 (MQ=255)                                                                                                                         ttggccgcacctctcAGGACGAGGAAGTCCAGTGTGATGTTCCCA‑TAATGCAGTACAAAATTTTCTCTGCTGATAACGGCGACGCATGGGTAGAATTTAAAGGCCAGAAAATGGCCCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATCCTTTAACGATGCTCAGCGTCAGGCAACCAAAGa                                                                                                                              <  1:414235/236‑1 (MQ=255)                                                                                                                             ccgccgcTTCCAGGACGAAGAAGTACAGCGTGATGTTTCAA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAg                                                                                                                         <  1:118848/250‑1 (MQ=255)                                                                                                                                        aGGAAGAAGACGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGACGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAACTGGGAGAACCGGAAACAGAACCTGTTATCACCGTACCGGAATACTTTAACGATGCTCAGCGTCAGGCAACCAAAAACGCAGGACGCATCGCt                                                                                                              >  1:37520/1‑251 (MQ=255)                                                                                                                                           gcaaagaagTACAGCGGGATTTTCCCA‑TCATGCCTTTCAAAATAATTGCGGTTAAAAACGGCGACGCAGGGGTAAAATTAAAAGCCCAGAAAATGGCCCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGt                                                                                                           <  1:49330/248‑1 (MQ=255)                                                                                                                                               aaaagtacgggtggttttttCCA‑TCAGGCCGTGAAAAATTATTGCTCCGAATAACGGCGACGCAGGGGTAGAAGTAAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACATGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTgg                                                                                                       <  1:22173/234‑1 (MQ=255)                                                                                                                                                             aTGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGCCGCAAGCCGTATCGCTGGTCTGGAAGTAAAACGTatc                                                                                         >  1:126536/1‑251 (MQ=255)                                                                                                                                                                tttCCAATCATGCCGTTAAAAGTAATGGCTGCTGATAACGGGGAAGAATGGGTGAAAGTAAAAGGCCAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTatca                                                                                        <  2:50247/250‑1 (MQ=255)                                                                                                                                                                   ccA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                                                                                                                                                   >  2:2224/1‑171 (MQ=255)                                                                                                                                                                   ccA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAAc                                                                                   >  1:14311/1‑251 (MQ=255)                                                                                                                                                                   ccA‑TCAGGCTGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCAAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                                                                                                                                                   <  1:2224/171‑1 (MQ=255)                                                                                                                                                                     a‑tcatGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAACAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTCTCACCGTACCGGCATACTTTAACGATGCTCAGCGTCCGGCAAcc                                                                                                                                   >  1:275579/1‑201 (MQ=255)                                                                                                                                                                         atGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGGATCATCAACGAacc                                                                              >  1:118273/1‑251 (MQ=255)                                                                                                                                                                              gTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGCTTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCACAGCGTCAGGCAACCAAAGACGCAGGCCATATCGCTGGTCTAGAAGGAAAACATATCATCAACGAaacgac                                                                           >  2:271454/1‑249 (MQ=255)                                                                                                                                                                                          tgctgcTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAg                                                                                                                               >  1:20478/1‑185 (MQ=255)                                                                                                                                                                                                 gATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGAAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACgg                                                      >  1:355893/1‑251 (MQ=255)                                                                                                                                                                                                                   gTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAACGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCCGCTGCGCTGGCTTAAGGTCTGGACAAAggtcatgg                                    >  2:306924/1‑245 (MQ=255)                                                                                                                                                                                                                         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)                                                                                                                                                                                                                                                 ccgcAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTggg      <  1:313007/251‑1 (MQ=255)                                                                                                                                                                                                                                                     aGATCTCTGCTGAAGTGCTGAACAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGCAGTAACACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAATGTACTGGCAACCGTACTATCGCGGTTTATGACCTGggtggg  >  2:183988/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       ATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCA‑TCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGT  >  NC_000913/12260‑12753

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