Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,846	G→T	G228G (GGG→GGT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,846	0	G	T	94.7%	157.3 / ‑6.5	59	G228G (GGG→GGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (2/0);  major base T (24/30);  minor base A (3/0);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

AACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTG  >  NC_000913/12601‑13090                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      aaCGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGgtggtg                                                                                                                                                                                                                                                   >  2:482590/1‑250 (MQ=255)     gtgcacAGCGTCAGGCCACCAAAGATGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACTGTTTGGACAAAGGTACTGTCAACCGTACTAGCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑Ac                                                                                                                                                                                                                                             <  2:23396/246‑1 (MQ=255)       gCTCAGAGTCAGGCAACCAAAGACGCAGACCGTATCGCTGGTCTGGAAGTAAAACGTATCATAAACGAACCGACCGCAGATGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGGACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATCCCCACCTGGGTGGGGAAG‑ACtt                                                                                                                                                                                                                                           >  2:550981/1‑251 (MQ=255)          cAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCg                                                                                                                                                                                                                                         >  2:543804/1‑250 (MQ=255)              gTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGc                                                                                                                                                                                                                                    >  2:41494/1‑251 (MQ=255)                      aCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGAt                                                                                                                                                                                                                            >  1:390946/1‑251 (MQ=255)                      aCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGAAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGAACGCCGTCTGAt                                                                                                                                                                                                                            >  1:250862/1‑251 (MQ=255)                           aGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACt                                                                                                                                                                                                                       >  1:265496/1‑251 (MQ=255)                                  ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGt                                                                                                                                                                                                                <  1:543804/251‑1 (MQ=255)                                      gTATCGCTGGTCTGGAAGTTAAACGTATCATCAACGAACCGACAGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTga                                                                                                                                                                                                             >  1:485689/1‑250 (MQ=255)                                               gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTCCCTGGTTGAAGAATTCaag                                                                                                                                                                                                   <  2:18440/251‑1 (MQ=255)                                                   ggAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTTCTGTCTGGACAAAGGTTCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGACATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑TCTTCGACAGCCCTCTTATCAACTACCTGGTTGAACAATTCaagaa                                                                                                                                                                                                 <  1:378525/249‑1 (MQ=255)                                                    gAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTTAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCaagaaaga                                                                                                                                                                                              <  1:47119/251‑1 (MQ=255)                                                      aGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATc                                                                                                                                                                                            <  1:198382/251‑1 (MQ=255)                                                       gTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATc                                                                                                                                                                                            >  1:625801/1‑250 (MQ=255)                                                         aaaaCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAgg                                                                                                                                                                                         <  2:55209/251‑1 (MQ=255)                                                                  catcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGAc                                                                                                                                                                                 <  2:390946/250‑1 (MQ=255)                                                                     caACGAACCGACCGCAGCTGCGCTGGGTTACGGTCTTGACAAAGGTACTGGCAACCGTACTATCGCGGTTGATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTGCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTg                                                                                                                                                                              <  1:84621/250‑1 (MQ=255)                                                                      aaCGACCCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCCACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATAGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAACCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTgcg                                                                                                                                                                            >  2:206573/1‑251 (MQ=255)                                                                         gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCACCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGCGAAG‑ACTTCGACAGCCGTCCGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCa                                                                                                                                                                          >  2:443930/1‑250 (MQ=255)                                                                                  gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCt                                                                                                                                                                <  1:612287/251‑1 (MQ=255)                                                                                  gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGGATTGACCTGCGCAACGATCCGCt                                                                                                                                                                <  2:102943/251‑1 (MQ=255)                                                                                               cTTACGGTCTCGACAAAGGTACTGGCAACCGTACTATCGCGGTTTGTGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGGTACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGCTGAAGAATTCAAGAAAGCTCAGGGGATTGACCTGCGCAACGCTCCGCTGGCAATGCAgcgc                                                                                                                                                   >  2:274912/1‑251 (MQ=255)                                                                                                      tCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGCCGAAGTTGACGGCGAACAAACCTTCGAAGTTCTGGCCACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGTATTGACCTGAGCAACGATCCGCTGGCAATGCAGCGCCTgca                                                                                                                                              >  2:543454/1‑247 (MQ=255)                                                                                                          gACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGc                                                                                                                                        <  2:567615/251‑1 (MQ=255)                                                                                                              aaGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCa                                                                                                                                    <  2:557460/251‑1 (MQ=255)                                                                                                                   aCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGaaaa                                                                                                                               <  2:625801/251‑1 (MQ=255)                                                                                                                      ggCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGaaaaa                                                                                                                              <  1:456812/249‑1 (MQ=255)                                                                                                                              tACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCg                                                                                                                    >  1:62058/1‑251 (MQ=255)                                                                                                                                 tATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGGATTGACCAGCGCAACGATCCGCTGGCAATGCAGCGCCTGCAAGAAGCGGCAGAAAAAGCGAAAATCGAAc                                                                                                                 >  2:137802/1‑251 (MQ=255)                                                                                                                                      cgGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCt                                                                                                            >  2:6800/1‑251 (MQ=255)                                                                                                                                            aTGACCTGGGTGGTGGTACTTTCGATAGTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCt                                                                                                      <  2:249344/251‑1 (MQ=255)                                                                                                                                                ccTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGAAAAAACCTTCGAAGATCTGGCAACCAACGGTGATACCCACCTGGGAGGTGAAG‑ACTTCCACAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGACGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                  >  1:275980/1‑251 (MQ=255)                                                                                                                                                ccTGGGTGGTGGTACTTTCGATATTTCTATTAGCGAAATCGACGAAGTGGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagc                                                                                                  <  2:281324/251‑1 (MQ=255)                                                                                                                                                  ggggTGGGGGTTCTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACTGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTcagcag                                                                                                <  2:336696/250‑1 (MQ=255)                                                                                                                                                    ggtggtggtACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGAc                                                                                              >  1:257232/1‑251 (MQ=255)                                                                                                                                                            tACTTTCGATATTAATATTATCGAAATCGACGAAGTTGCCGGAGACACAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGAGGTGCAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGtt                                                                                       >  2:419170/1‑250 (MQ=255)                                                                                                                                                             aCTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTa                                                                                      <  2:490508/250‑1 (MQ=255)                                                                                                                                                              atttCGATATTTCTATTATCGAAATCGACGAAGTAGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAAc                                                                                    <  1:294944/250‑1 (MQ=255)                                                                                                                                                                           tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAAt                                                                                                                                                          <  1:500646/168‑1 (MQ=255)                                                                                                                                                                           tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAAt                                                                                                                                                          >  2:500646/1‑168 (MQ=255)                                                                                                                                                                             ttatCGACATCCGCGAAGTTGACGACGCAAATACCTTCGAAGTTCTGAAAACCAAAGATGATACCCAGATGGGAGGTAACG‑CCTTCGACAGCCGTCTGCTCAACTCCATGGTTGAAGAATTCAAGAAAGCTCAGGGCATTGACCTGCGCAACGATCCGCTGGCACTGCAGCGCCTGAAAGTAGCGGCAGCAAAAACGAAAATCGACCTGTCTTCCGCTCAGCAGCCCGACGTTCACCTTCCGAACATCACt                                                                     >  2:309762/1‑251 (MQ=255)                                                                                                                                                                             ttatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                     >  1:484449/1‑251 (MQ=255)                                                                                                                                                                                 cGAAATCGCCGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGACGTACATCACTGc                                                                   >  2:394471/1‑249 (MQ=255)                                                                                                                                                                                    aaTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGatg                                                              >  2:286925/1‑249 (MQ=255)                                                                                                                                                                                     aTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTGTATTTACCTGCGCAACGATCCGCTGTCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGatgc                                                             <  1:591604/251‑4 (MQ=255)                                                                                                                                                                                        gacgaAGTTGACGGCGAAAAAACCTTCGACGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGAc                                                          >  2:185390/1‑251 (MQ=255)                                                                                                                                                                                          cgaATTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTTCCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTACCCTGCCGTACATCGCTGCAGATGCGACCg                                                        <  1:482719/251‑1 (MQ=255)                                                                                                                                                                                               ttGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGATGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGACTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCg                                                   >  1:375631/1‑251 (MQ=255)                                                                                                                                                                                                    ggCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGACGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCGCTGCAGATGCGACCGGTCCGAAaca                                              <  2:349468/251‑1 (MQ=255)                                                                                                                                                                                                     gCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTACCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAacac                                             <  2:162937/251‑1 (MQ=255)                                                                                                                                                                                                                  cGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAACGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGCTGCGACCGGTCCGAAACACATGAACATCaa                                  >  1:359086/1‑249 (MQ=255)                                                                                                                                                                                                                      gTTCTGTCAACCAACGGTGAGACCCACTTGGGTGGTGAAG‑ACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATGGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACGGTCTTCCGCGCAGCAGACCGACGTTAACCTGCCGTACATCACTGGAGATGCGACCGGTCCGAAACACAGGAACATCAAAGttac                            <  2:62050/251‑4 (MQ=255)                                                                                                                                                                                                                      gTTCTGGCAACCAACGGTGATACCCCCCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGGTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGAc                            <  2:53247/251‑1 (MQ=255)                                                                                                                                                                                                                           ggCAACCAACGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCg                         <  1:6800/249‑1 (MQ=255)                                                                                                                                                                                                                                  aacGGTGATACCCACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACt                >  1:79777/1‑251 (MQ=255)                                                                                                                                                                                                                                        gatccccACCTGGGTGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTACCGGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACAGGAACATCAAAGTGACTCGTGCGAAACTGGAAAg          <  2:359086/247‑1 (MQ=255)                                                                                                                                                                                                                                        aaaaCCCACCTGGGTGGTGAAG‑ACTTCGAATGCCGTCTTAGCAACTTCCTGGGTGAAAAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAATTGACTCGTGCGAAACTGGAAAg          <  2:80329/248‑1 (MQ=255)                                                                                                                                                                                                                                            cccACCTGGGGGGAGAAGCACT‑CGCCACCAGTCTGATCACCTACCTGGTTGAAGAATTCAGGAACGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGACAAAGCGAAAATCGAACTGTCTTCCGCTAAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACCTGAACATCAAAGTGCCTCGTGCGAAACTGGGAAGCCTg      >  2:621188/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cTGGGTGGTGAAG‑ACTTCGCCAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATAAGGGTCTTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTg  >  2:62082/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      AACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAG‑ACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTG  >  NC_000913/12601‑13090

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