Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,879	T→C	Y239Y (TAT→TAC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,879	0	T	C	94.3%	197.1 / 1.8	70	Y239Y (TAT→TAC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (38/28);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

TCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGC  >  NC_000913/12641‑13117                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               ccGCTGGTCGGGAATTAAAACGTAACATCAACGAACCGACCGCAGCTGCGCTGGTTTACGGTCTGGAAAAAGGTACTGGCAACCGGACTATCGCGGTTTATGCCCTGGGTGGTGGTACTTTCTATATTTCTATTATCTAAATCGACGAAGTTGACGGCGAAAAAACCTTCGACGTTCTGTCAACCAACGGTGATACCCACCTGGTTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTtaagaa                                                                                                                                                                                                                                    <  2:449273/250‑1 (MQ=255)  cacgtgtctagaacgaaaacggttcatcaTCGAACCGACCGCAGCTGCGCTGTCTTACTGTCTGTAAACATGTACGTGCAACCTTACCATCGCGGTTTAGGACCTGGGTTGGGGTACATACGATATTTCTATTACCGAAATCGACGAAGTAGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGTTGACATCCCCCTGGGTGGTGAATACGTCGACAGTCGGCTTAGCAACTTCCTGGTTgaagaa                                                                                                                                                                                                                                    <  2:330456/227‑1 (MQ=255)      ggTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCa                                                                                                                                                                                                                                >  1:358995/1‑250 (MQ=255)       gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTCTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAACCCTTAGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGAGGTGAAGACTTCGACAACCGTCTGATCACCTACATGGTTGACGAAttaa                                                                                                                                                                                                                                >  2:442594/1‑247 (MQ=255)          tGGAAGTAAAACGTATCATCCACGAACCGACCGCAGCTGCGCTGGCATAAGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGCTATTTCTATTATCGAAATCGACGACGATGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGATGATACCAACCTGGGTGGTGAAAACTCCGACAGCCGGCTGATAAACAACCTGATTGAAGAATGCaagaa                                                                                                                                                                                                                            >  2:481796/1‑250 (MQ=255)                       tATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATAGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGCAGTTCTGGAAACAAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATAAACTACCAGGTGGAAGGATTAAAGAAAGATCAAGGCATTg                                                                                                                                                                                                              >  2:456975/1‑251 (MQ=255)                          catcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGAcc                                                                                                                                                                                                           <  2:111663/251‑1 (MQ=255)                           agcaACGAACCGACCGCAGCTGCGGTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTAGCGCGGTTTATGACCTGGGTGGTGGTACTTTAGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCGGGTTGAATAATTCAAGAAAGATCAGGGCATTGACCt                                                                                                                                                                                                          <  2:322370/249‑1 (MQ=255)                              aaCGAACCGTCCTCAGCTGCGCTGGTTTACGGTCTGGCCAAAGGTACTTGCAACCGTACTATCGCGGTTTATGGCCTGGGTGGTGGTATTTACGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACGTCGACAGCCTTCTGAGCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTTCATTGGCCTgcg                                                                                                                                                                                                       <  2:652953/251‑1 (MQ=255)                                cGAACCGACCGCAGCTGCGCTGGCTTACGGGCTGGACAAAAGTACTGGCACCCGGACTCTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGGCGGCGAAAAACCCTTCGAAGTTCTGGCACCCCACGGTGATACCCACCTGGGTGGGGACGACTTCGACCGCCATCTGATCCACTACCTGGTTGAAGTATTCaggaagtataagggacatgacctgcgc                                                                                                                                                                                                      >  2:174271/1‑228 (MQ=255)                                 gAACCGACCGCAGCCGCGCGGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTCCGACAGCCGTCTGATCAACTACCTGGTTGAAGAAt                                                                                                                                                                                                                                   <  1:101373/220‑1 (MQ=255)                                   accgaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACg                                                                                                                                                                                                  >  1:170583/1‑251 (MQ=255)                                                 cgcTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATg                                                                                                                                                                                    >  1:457894/1‑251 (MQ=255)                                                  gcTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGc                                                                                                                                                                                   >  1:426652/1‑251 (MQ=255)                                                     ggcgtgcGGTCTGGACAATGGTACTGGCAAACGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTTACGGCTAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAgc                                                                                                                                                                                <  2:678034/245‑1 (MQ=255)                                                     ggCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAg                                                                                                                                                                                 <  1:56840/250‑1 (MQ=255)                                                       cTTACGGTCTGGACAAAGGTACTGGCAAACGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTAACGACATCGACGAAGTTGGCGGCGAAAAAACCTTCGAAGTTCTGGCCACCAACGGTGATACACACCTGGGTGGTGAAGACTTAGACCGCCGTCTGATCAGATACCTGGTTGTAGAATGCAAGGAAGGTCAAGACATTGACCTGAGCAACGacccactgaccaagcaactc                                                                                                                                                                              >  2:312103/1‑232 (MQ=255)                                                                   aCAAAGGTACTGGCAACCGTAATATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGATTTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCg                                                                                                                                                                  >  1:101112/1‑251 (MQ=255)                                                                      aaGGTACTGGCAACCGTACTATCGCGGTTTATGCCCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCa                                                                                                                                                               <  1:28539/251‑1 (MQ=255)                                                                               gCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCg                                                                                                                                                      <  1:650643/251‑1 (MQ=255)                                                                                cAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGAACAAGGCATTGACATGCGCAACGAACAGCTGGCATTGCATCGCCTGACAGAAGCGGCAGaaaaatcg                                                                                                                                                      >  2:283169/1‑247 (MQ=255)                                                                                 aaCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGaa                                                                                                                                                    <  1:314256/251‑1 (MQ=255)                                                                                   ccGTACTATCGCGGTTTATGACCTGGGGGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTCCCTGGCTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCACCGCCTGAAAGAAGCGGCAGAAAAAGCGaaa                                                                                                                                                   >  2:91002/1‑250 (MQ=255)                                                                                          aTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGTTGATACCCACCTGGGTGTTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAAc                                                                                                                                            <  2:626074/250‑1 (MQ=255)                                                                                          aTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAAc                                                                                                                                            <  1:280995/250‑1 (MQ=255)                                                                                               ggTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCtt                                                                                                                                      <  2:170583/251‑1 (MQ=255)                                                                                                       aCCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGACAAAACCTTCGAAGTTCTGGCAACCAACGGCGATACCCACCTGGGTGGTGAAGACTTCGACAGCAGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAACAGCGAAAATCGAACTGTCTTCCGCTc                                                                                                                                >  1:641777/1‑249 (MQ=255)                                                                                                                 tggtACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAAAGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACg                                                                                                                    >  2:429219/1‑251 (MQ=255)                                                                                                                  tgtACTTTCGATATTTCTATTATCAAAATCGACGACGTTGCCGGCGAAAAAACCTTCGAAGTTCTGGCACCCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGACAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCCGACAAAGCGAAAATCGAACTATCTTCCCCCCAGCcgaccgacct                                                                                                                   >  2:344765/2‑249 (MQ=255)                                                                                                                     actatCGAGATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAATGCCCTCGAAGTTCTGACAACCAACGGTGATACCCACCTGGGCGGTGAAGACTTCGACAGCCTTCTGACCAACTACCTGGTTGAATAATTCAAGAAATATCAGTTAATTGACCGGCGGAACAATCCGCGCGCAATGCAGCGCCCCTAAGTCGCGGCACAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTaa                                                                                                                <  1:403627/247‑1 (MQ=255)                                                                                                                       tttCGATCTTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGACGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGCAGACTTCGACAGCCGTCTGATCAACTACATGGTTGAAGCATTCAAGAAAGCTCAGGGCATTGACCTGCGCAACGATCCGATGGCAATGCAGCGCCTTAAAGAAGCGGCAGAAAAAGCGAAAATCGAAATGGCTTCCGACCAGCAACCCGACGtacaac                                                                                                              >  2:8678/1‑246 (MQ=255)                                                                                                                        ttCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCt                                                                                                             >  2:567715/1‑251 (MQ=255)                                                                                                                          cGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCAACCTGGGAGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTAAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCCAAAATCGAACTGTCTTCCGCTCCGCAGACCGACGTTAAccggc                                                                                                           >  2:44820/1‑248 (MQ=255)                                                                                                                                   tattatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAATGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATc                                                                                                   >  1:211237/1‑250 (MQ=255)                                                                                                                                     ttatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCATCTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAg                                                                                                                                                              <  1:680384/189‑1 (MQ=255)                                                                                                                                     ttatCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                                                >  1:571416/1‑251 (MQ=255)                                                                                                                                      taacGAAATCGACGAAGGTGACGGCGAAAAAACCGACGAAGTTCTGGCAACCAACTGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACt                                                                                                <  2:130886/247‑1 (MQ=255)                                                                                                                                       atCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGc                                                                                              >  2:286867/1‑251 (MQ=255)                                                                                                                                                    aaGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCCACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGt                                                                                 >  1:639562/1‑251 (MQ=255)                                                                                                                                                        tGACGGCGAAAAAACCTTCGAAGTGCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCg                                                                              <  2:217711/250‑1 (MQ=255)                                                                                                                                                          ccctacacgacactctcccgactttGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGAAAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGTAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                           <  2:369292/227‑1 (MQ=255)                                                                                                                                                            ggCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGAt                                                                                                                                                                                                >  1:551996/1‑132 (MQ=255)                                                                                                                                                            ggCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGAt                                                                                                                                                                                                <  2:551996/132‑1 (MQ=255)                                                                                                                                                                   aaaCCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAc                                                                  >  1:667448/1‑251 (MQ=255)                                                                                                                                                                    aaCCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACa                                                                 <  1:456975/251‑1 (MQ=255)                                                                                                                                                                     aCCTTAGAAGTTCTGGCAACCAAAGGTAACACCAAACTGGGTGGTGACGACTCCGACCGCAGTCAGATCAACTACCTGGTTGAAGAATTCAAGAAAGCTCAGGGCATTGACCTGCGCAACGATACGCTGGCAATGCAGCGCCTGCAAGAAGAAGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGCCCGCAGTTAACCAGCCGTCAATCACTGCAGATGCGACCGCCACGAAACACATGaaat                                                                 >  2:659860/1‑248 (MQ=255)                                                                                                                                                                           gAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGt                                                          >  1:294601/1‑251 (MQ=255)                                                                                                                                                                                  tGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                           >  1:369292/1‑227 (MQ=255)                                                                                                                                                                                      aaccaacGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCg                                                >  1:326817/1‑250 (MQ=255)                                                                                                                                                                                         caacGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAAc                                            <  1:286867/251‑1 (MQ=255)                                                                                                                                                                                             ggTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATAAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCCGAAAAAGCGAAAATCGACCTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGaaccggg                                         >  2:306015/1‑245 (MQ=255)                                                                                                                                                                                              cctaTACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGACCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAA‑GCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                      <  1:342057/248‑1 (MQ=255)                                                                                                                                                                                               ggATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                      <  1:567715/250‑1 (MQ=255)                                                                                                                                                                                                gATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGACAGATCAGGGCATTGACCTGCGCAACGATACGCTGGAAAAGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGACAATCGCACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACAAATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                      >  2:399771/1‑250 (MQ=255)                                                                                                                                                                                                 aTACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGCAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAA‑GCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGaaa                                      >  2:342057/1‑248 (MQ=255)                                                                                                                                                                                                   aCCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAACATCGAACTGTCTTCCGCTCAGCAGACAGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGCAACACATAAACATCAAAGAGACTCGTGCGCAACTGGACAGCCt                                  >  2:284925/1‑251 (MQ=255)                                                                                                                                                                                                    cccACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTg                                 <  1:49258/251‑1 (MQ=255)                                                                                                                                                                                                      cACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGt                               >  1:427149/1‑251 (MQ=255)                                                                                                                                                                                                       cccTGGGTGGGGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAACTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGtt                              <  1:481796/250‑1 (MQ=255)                                                                                                                                                                                                       aCCTGGGTGGTGCAGACTTCGACAGCCGTCTAATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGAACCGCTGGCAATGCAGCGCCTGCAAGAAGCGGCAGAACAAGCGAAAATCGAACTGTCTTCCGCTCCGCAGGCCGCCGTTACCCTTCCGTACATCACTGCCGATGCGACCGGTCCGAAACAAATGACAATCAAAGTGACTCGTGCGAAACAGGAAAGCCTGGtt                              >  2:673158/1‑251 (MQ=255)                                                                                                                                                                                                        ccTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTg                             <  1:429219/251‑1 (MQ=255)                                                                                                                                                                                                             gtggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGa                            >  1:214669/3‑247 (MQ=255)                                                                                                                                                                                                             gtggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGAt                        <  1:35365/249‑1 (MQ=255)                                                                                                                                                                                                               ggTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTACGAAACACATTAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTAGAAGATCt                      >  2:324918/1‑251 (MQ=255)                                                                                                                                                                                                                     gACTTCGCCAGCCGTCTGATCAACTACCTGGTTGAAGAATTCACGACCGATCAGGGCATTGACCTGCGCAACGATCAGCTGGCAATGCAGCGCCTGAAAGAAGCGGAAGAAAAAGCGAAAATAGAACTGTCTTCCGCTCAGCAGACCGCCGTTAACCTGCCGTACATCACTGCAGCAGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGAGAAACTGGAAAGCCTGGTTGAAGATCTggaaa                 >  2:435838/1‑247 (MQ=255)                                                                                                                                                                                                                         tCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCg             <  1:399771/250‑1 (MQ=255)                                                                                                                                                                                                                                ccGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGAGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTg     >  1:487153/1‑251 (MQ=255)                                                                                                                                                                                                                                   tCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCGGGTTGAAGATCTGGTAAACCGTTCCATTgcgc  >  2:225634/1‑248 (MQ=255)                                                                                                                                                                                                                                     tGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATCGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCa                                                                                                  >  1:85603/1‑153 (MQ=255)                                                                                                                                                                                                                                     tGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATCGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCGTACATCa                                                                                                  <  2:85603/153‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               TCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACCTGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGC  >  NC_000913/12641‑13117

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