Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,570	A→G	V28A (GTG→GCG)	folE ←	GTP cyclohydrolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,570	0	A	G	100.0%	153.7 / NA	43	V28A (GTG→GCG)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/18);  total (25/18)

AAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTG  >  NC_000913/2243426‑2243713                                                                                                                                                 |                                                                                                                                                aaaaTTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGcgggcg                                                                                                                                           >  1:167341/1‑151 (MQ=255)         aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGc                                                                                                                                   <  2:112652/151‑1 (MQ=255)         aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGc                                                                                                                                   <  1:340739/151‑1 (MQ=255)                          cgatgcgatgCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTcg                                                                                                                  >  2:535841/1‑151 (MQ=255)                            atgcgatgcgGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTcgcg                                                                                                                <  1:676330/151‑1 (MQ=255)                                   gcggcgTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTa                                                                                                          <  1:80165/150‑1 (MQ=255)                                     ggcgTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACg                                                                                                       >  2:150326/1‑151 (MQ=255)                                         tttCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTc                                                                                                   >  2:320614/1‑151 (MQ=255)                                          ttCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTc                                                                                                   >  1:630699/1‑150 (MQ=255)                                             cATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATg                                                                                                >  2:265613/1‑150 (MQ=255)                                                cAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAAcc                                                                                            >  2:138257/1‑151 (MQ=255)                                                  aaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAAcc                                                                                            >  1:326808/1‑149 (MQ=255)                                                     tGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGc                                                                                       >  1:733461/1‑151 (MQ=255)                                                        catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGAGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGc                                                                                    >  2:164142/1‑151 (MQ=255)                                                                aGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGc                                                                                    <  1:671004/143‑1 (MQ=255)                                                                aGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGc                                                                                    >  2:670964/1‑143 (MQ=255)                                                                    cGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGa                                                                         <  2:759655/150‑1 (MQ=255)                                                                     gagaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGt                                                                       >  1:584760/1‑151 (MQ=255)                                                                       gaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTg                                                                      <  2:639674/150‑1 (MQ=255)                                                                       gaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGa                                                                     <  1:628429/151‑1 (MQ=255)                                                                         ttCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                   >  1:478869/1‑151 (MQ=255)                                                                                      tGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAttt                                                      >  1:130882/1‑151 (MQ=255)                                                                                        aTTTCGGTCATATGACAAGCAATAAGGCTTTTGCGCGTTTAGTTATCCATTTCATGCGCGGGCGGGTGCAGCGGTGTTTCCAGTCCTCGTGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtct                                                    >  2:766071/1‑151 (MQ=255)                                                                                           tCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTcc                                                  <  2:608785/150‑1 (MQ=255)                                                                                            cGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTg                                                >  2:57289/1‑151 (MQ=255)                                                                                             ggTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGc                                               <  1:435168/151‑1 (MQ=255)                                                                                                cATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAg                                             <  2:167334/150‑1 (MQ=255)                                                                                                cATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAg                                             <  1:272451/150‑1 (MQ=255)                                                                                                    tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgtgt                                         >  1:473775/1‑150 (MQ=255)                                                                                                           cAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTc                                 >  2:34993/1‑151 (MQ=255)                                                                                                            aaTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                <  2:474861/151‑1 (MQ=255)                                                                                                               aaGGCTTTTGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGcc                             <  2:630662/151‑1 (MQ=255)                                                                                                                    ttttGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCg                        >  1:590258/1‑151 (MQ=255)                                                                                                                     tttGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                <  2:472730/142‑1 (MQ=255)                                                                                                                     tttGCGCGTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                >  1:472752/1‑142 (MQ=255)                                                                                                                          gcgTTTCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggggg                  >  1:142542/1‑151 (MQ=255)                                                                                                                               tCGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCaaa              <  1:55871/150‑1 (MQ=255)                                                                                                                                cGTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt            <  2:664701/151‑1 (MQ=255)                                                                                                                                 gTTATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGtt           <  2:650759/151‑1 (MQ=255)                                                                                                                                  ttATCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGtt           >  2:360211/1‑150 (MQ=255)                                                                                                                                    aTCCATTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGttatt        >  2:59677/1‑151 (MQ=255)                                                                                                                                        aTTTCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGc    >  1:44486/1‑151 (MQ=255)                                                                                                                                          ttCATGCGCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTg  >  2:3817/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                AAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTG  >  NC_000913/2243426‑2243713

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