Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,977	T→G	V366G (GTG→GGG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,977	0	T	G	97.7%	152.2 / ‑2.9	43	V366G (GTG→GGG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (20/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCT  >  NC_000913/1931832‑1932122                                                                                                                                                  |                                                                                                                                                  tATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCAc                                                                                                                                              >  1:86687/1‑151 (MQ=255)          gCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACt                                                                                                                                      <  2:125468/150‑1 (MQ=255)             gCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCag                                                                                                                                  <  1:777109/151‑1 (MQ=255)                 cTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCagag                                                                                                                                >  2:235081/1‑149 (MQ=255)                      gTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGtgttg                                                                                                                         >  2:36295/1‑151 (MQ=255)                             tGCCACAACTGACCAGTCATAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGc                                                                                                                  >  2:257369/1‑151 (MQ=255)                               ccACAACTGACCAGTCAGAATGTCTCGTTTGATAATATGCAGAATGCCGTAGGCGCTGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCATTGGCTACTGCAGTGAGTGTTGTTGACGc                                                                                                                  >  1:717457/1‑149 (MQ=255)                                  caACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTg                                                                                                             >  2:423320/1‑151 (MQ=255)                                   aaCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTg                                                                                                             >  1:719072/1‑150 (MQ=255)                                   aaCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGa                                                                                                            <  2:264018/151‑1 (MQ=255)                                           aGTCAGCATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcgc                                                                                                      <  2:748648/149‑1 (MQ=255)                                            gTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGaa                                                                                                   <  1:235118/151‑1 (MQ=255)                                               aGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGc                                                                                                 >  2:367530/1‑150 (MQ=255)                                                 aaTGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcg                                                                                              >  1:559104/1‑151 (MQ=255)                                                    gTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCAcgcc                                                                                            <  2:528292/150‑1 (MQ=255)                                                      cACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACgccgcc                                                                                         <  1:587958/151‑1 (MQ=255)                                                         gTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCgg                                                                                       >  1:321014/1‑150 (MQ=255)                                                             gATAATGTGCAGAGTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTAGGGGGGGCTCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCTCCGGACa                                                                                    <  2:539116/149‑1 (MQ=255)                                                                    tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGt                                                                           >  1:630734/1‑151 (MQ=255)                                                                     gCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGt                                                                           <  2:136400/150‑1 (MQ=255)                                                                          aTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAggg                                                                     <  1:699084/151‑1 (MQ=255)                                                                             ccGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTa                                                                  <  1:760974/151‑1 (MQ=255)                                                                                  ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAAccc                                                             <  2:36998/151‑1 (MQ=255)                                                                                   gcgcAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCg                                                            >  1:671807/1‑151 (MQ=255)                                                                                   gcgcAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCg                                                            >  1:760709/1‑151 (MQ=255)                                                                                   gcgcAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCg                                                            >  2:78776/1‑151 (MQ=255)                                                                                     gcAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGc                                                          <  2:552532/151‑1 (MQ=255)                                                                                       aGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCa                                                         >  1:293907/1‑150 (MQ=255)                                                                                                 aTTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGGAAAAGGACAGGGGTAAACCCGGCAAAAAAACg                                                 >  2:308617/1‑148 (MQ=255)                                                                                                    cggttAATTAGTAAGCCGGAAAATGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACAGCAGAGAGTGATGTTGACGCCATTGAACGCGCGAAACACGCCGTCGGACAGGTAAAAGTACAGGTTTAAACCCGGCAAAAAAACGCTACa                                            <  2:517217/147‑1 (MQ=255)                                                                                                    cGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACa                                            <  2:605932/150‑1 (MQ=255)                                                                                                     gTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaaa                                          <  2:241097/151‑1 (MQ=255)                                                                                                       ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaaaa                                         <  2:325813/150‑1 (MQ=255)                                                                                                       ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACaaaa                                         <  2:325822/150‑1 (MQ=255)                                                                                                            tGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCTCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGccc                                   >  2:309989/1‑151 (MQ=255)                                                                                                            tGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGccc                                   >  1:307289/1‑151 (MQ=255)                                                                                                             ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCg                                  <  2:320974/151‑1 (MQ=255)                                                                                                             ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCg                                  <  1:315340/151‑1 (MQ=255)                                                                                                               tAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGAt                                >  2:511472/1‑151 (MQ=255)                                                                                                                aaGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATc                               <  1:641961/151‑1 (MQ=255)                                                                                                                  gCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCt                             <  1:390432/151‑1 (MQ=255)                                                                                                                           tGATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                     <  2:130005/150‑1 (MQ=255)                                                                                                                            gATGGCAGCCGTCGTCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                    >  2:355228/1‑150 (MQ=255)                                                                                                                                             ggggggggCCCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCt  >  2:546044/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCT  >  NC_000913/1931832‑1932122

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