Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,036	G→T	S206Y (TCC→TAC)	folE ←	GTP cyclohydrolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,036	0	G	T	97.1%	121.9 / ‑2.9	35	S206Y (TCC→TAC)	folE	GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (24/10);  total (25/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACG  >  NC_000913/2242978‑2243097                                                           |                                                              ttttAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGa                                                         <  1:3357631/65‑1 (MQ=255)       tCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATc                                                   <  1:2212471/65‑1 (MQ=255)         aGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATcca                                                 <  1:928393/65‑1 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTCTTCTGACTGTATTTGAAGAATCCACCAagaga                                         >  1:2030345/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAagaga                                         >  1:45866/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAagaga                                         >  1:129611/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAagaga                                         >  1:3310292/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAagaga                                         >  1:2772739/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAagaga                                         >  1:2647469/1‑65 (MQ=255)                 tgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGCATTTGAACAATCCACCAagaga                                         >  1:2025868/1‑65 (MQ=255)                   aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAgg                                       <  1:1299908/65‑1 (MQ=255)                   aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAgg                                       <  1:3377812/65‑1 (MQ=255)                   aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAgg                                       <  1:3331756/65‑1 (MQ=255)                    cGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGt                                      >  1:2423719/1‑65 (MQ=255)                      cacaGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCg                                    >  1:946300/1‑65 (MQ=255)                       acaGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCCCCAAGAGAGGTCGt                                   <  1:119526/65‑1 (MQ=255)                         aGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTg                                 >  1:2484463/1‑65 (MQ=255)                         aGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTg                                 >  1:248502/1‑65 (MQ=255)                         aGCGCGCAGAAACTCGCGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTg                                 >  1:2827853/1‑65 (MQ=255)                                 gAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCa                         >  1:1195225/1‑65 (MQ=255)                                 gAAACTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCa                         >  1:1463272/1‑65 (MQ=255)                                    aCTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTg                      <  1:2851751/65‑1 (MQ=255)                                     cTCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTgg                     <  1:2142947/65‑1 (MQ=255)                                     cTCGTGGCGCGTATTCTGACTGTATTTGAACAACCTACCAAGAGAGGACGCTGTCGTGGCACTgg                     >  1:3278779/1‑65 (MQ=255)                                      tCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGt                    <  1:1570876/65‑1 (MQ=255)                                      tCGTGGCGCGTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGt                    >  1:896813/1‑65 (MQ=255)                                            cgcgTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATc              >  1:3232604/1‑65 (MQ=255)                                            cgcgTATTCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATc              >  1:1999107/1‑65 (MQ=255)                                                   tCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATg       >  1:419689/1‑65 (MQ=255)                                                   tCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATg       >  1:565707/1‑65 (MQ=255)                                                   tCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATg       >  1:2278207/1‑65 (MQ=255)                                                   tCTGACTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGCTGCATCGCGGATg       >  1:708391/1‑65 (MQ=255)                                                       aCTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAc   >  1:2647541/1‑65 (MQ=255)                                                       aCTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAc   >  1:918542/1‑65 (MQ=255)                                                       aCTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAc   >  1:2288269/1‑65 (MQ=255)                                                        cTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGcctgt  >  1:1196234/1‑62 (MQ=255)                                                        cTGTATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCAcg  >  1:1338523/1‑65 (MQ=255)                                                           |                                                              TTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACG  >  NC_000913/2242978‑2243097

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