Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,373	C→T	Q165* (CAG→TAG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,373	0	C	T	100.0%	115.2 / NA	36	Q165* (CAG→TAG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/19);  total (17/19)

GGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGT  >  NC_000913/1931240‑1931506                                                                                                                                      |                                                                                                                                      ggcggcAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACg                                                                                                                                  <  1:19339/139‑1 (MQ=255)                  agctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCaga                                                                                                                 <  1:990371/139‑1 (MQ=255)                  agctgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCaga                                                                                                                 <  1:721591/139‑1 (MQ=255)                            ttccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCt                                                                                                       <  2:710397/139‑1 (MQ=255)                             tccacttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTc                                                                                                      <  2:766854/139‑1 (MQ=255)                                  ttATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCa                                                                                                 <  2:1315021/139‑1 (MQ=255)                                          ttGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTa                                                                                         >  2:472343/1‑139 (MQ=255)                                           tGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCt                                                                                                       >  2:634647/1‑124 (MQ=255)                                           tGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCt                                                                                                       <  1:634647/124‑1 (MQ=255)                                                          ccTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGt                                                                         <  2:71585/139‑1 (MQ=255)                                                              ttCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTcgcg                                                                     >  1:656033/1‑139 (MQ=255)                                                                   cgAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGa                                                                <  1:1324184/139‑1 (MQ=255)                                                                        cGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCAGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                           <  1:202438/139‑1 (MQ=255)                                                                             gCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgcg                                                      <  1:1142582/139‑1 (MQ=255)                                                                             gCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgcg                                                      <  1:510269/139‑1 (MQ=255)                                                                                       gCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGa                                            >  2:833680/1‑139 (MQ=255)                                                                                           tCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAg                                                                              >  2:1330478/1‑101 (MQ=255)                                                                                           tCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAg                                                                              <  1:1330478/101‑1 (MQ=255)                                                                                                 cATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGtt                                  >  1:1080485/1‑139 (MQ=255)                                                                                                     gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGt                              >  1:322414/1‑139 (MQ=255)                                                                                                     gTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGt                              >  2:653094/1‑139 (MQ=255)                                                                                                      tAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgcg                                                      <  1:1313178/114‑1 (MQ=255)                                                                                                      tAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCcgcg                                                      >  2:1313178/1‑114 (MQ=255)                                                                                                      tAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGtt                             >  1:543852/1‑139 (MQ=255)                                                                                                      tAAAACCGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGtt                             >  2:1049377/1‑139 (MQ=255)                                                                                                           ccGGTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACt                        >  1:261231/1‑139 (MQ=255)                                                                                                              gTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTc                                                                                                      <  1:245958/58‑1 (MQ=255)                                                                                                              gTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTc                                                                                                      >  2:245958/1‑58 (MQ=255)                                                                                                              gTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgt                                      <  1:713192/122‑1 (MQ=255)                                                                                                              gTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGcgt                                      >  2:713192/1‑122 (MQ=255)                                                                                                              gTGATGAGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCg                     >  2:256216/1‑139 (MQ=255)                                                                                                                    aGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                           >  2:655337/1‑95 (MQ=255)                                                                                                                    aGCTCTTCCGGCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCggg                                                           <  1:655337/95‑1 (MQ=255)                                                                                                                              gCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATg                                                                                               >  2:741750/1‑49 (MQ=255)                                                                                                                              gCAAGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATg                                                                                               <  1:741750/49‑1 (MQ=255)                                                                                                                                 aGGGGTAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGt  <  1:256216/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      GGCGGCAGAAGAGCTGCAGCTGCCCACTTCCACTTATCGTTTTGCCGATAGCGAAAGCCTTTTCCGCGAGGCGGTTGCTGACATTGGCTATCCCTGCATTGTAAAACCGGTGATGAGCTCTTCCGGCAAGGGGCAGACGTTTATTCGTTCTGCAGAGCAACTTGCTCAGGCATGGAAGTACGCTCAGCAAGGCGGTCGCGCCGGAGCGGGCCGCGTAATTGTTGAAGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGT  >  NC_000913/1931240‑1931506

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