Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,632,046	G→T	P164T (CCG→ACG)	guaA ←	GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,632,046	0	G	T	95.7%	154.7 / ‑1.0	46	P164T (CCG→ACG)	guaA	GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (22/22);  total (23/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 = 9.69e-01

GTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTC  >  NC_000913/2631910‑2632184                                                                                                                                         |                                                                                                                                           gTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGt                                                                                                                                >  6:134932/1‑149 (MQ=255)         ggCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAg                                                                                                                                        >  8:97149/1‑133 (MQ=255)         ggCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAg                                                                                                                                        <  7:97149/133‑1 (MQ=255)          gCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTAc                                                                                                                       <  6:135581/149‑1 (MQ=255)            ttCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGc                                                                                                                              >  1:126552/1‑140 (MQ=255)            ttCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGc                                                                                                                              <  2:126552/140‑1 (MQ=255)               acacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtga                                                                                                                  >  4:97370/1‑149 (MQ=255)                              aCGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGaa                                                                                                   >  3:134957/1‑149 (MQ=255)                                 cacaAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAg                                                                                                <  6:92510/149‑1 (MQ=255)                                            ccAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTa                                                                                     <  7:123876/149‑1 (MQ=255)                                             cAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacgga                                                                                                    >  7:71969/1‑133 (MQ=255)                                             cAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacgga                                                                                                    <  8:71969/133‑1 (MQ=255)                                                cATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGaa                                                                                                              <  2:19660/120‑1 (MQ=255)                                                cATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGaa                                                                                                              >  1:19660/1‑120 (MQ=255)                                                       aTACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCt                                                                          >  8:144098/1‑149 (MQ=255)                                                        tACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTgg                                                                                                                           >  2:113778/1‑99 (MQ=255)                                                        tACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTgg                                                                                                                           <  1:113778/99‑1 (MQ=255)                                                           ccTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATc                                                                      <  3:97370/149‑1 (MQ=255)                                                                 gggTATGAGTCACTTCCGGGTGGAACTGTACGACATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACa                                                                >  2:49149/1‑149 (MQ=255)                                                                  ggTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACAt                                                               >  1:162714/1‑149 (MQ=255)                                                                   gTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                            >  5:154518/1‑135 (MQ=255)                                                                   gTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                            <  6:154518/135‑1 (MQ=255)                                                                      tGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGag                                                           <  7:144098/149‑1 (MQ=255)                                                                          tCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGagcagc                                                       >  2:42489/1‑149 (MQ=255)                                                                            aCTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTc                                                                                                           >  1:113780/1‑95 (MQ=255)                                                                            aCTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTc                                                                                                           >  1:116551/1‑95 (MQ=255)                                                                            aCTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTc                                                                                                           <  2:113780/95‑1 (MQ=255)                                                                            aCTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTc                                                                                                           <  2:116551/95‑1 (MQ=255)                                                                             cTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGaa                                                                                                   <  4:23027/102‑1 (MQ=255)                                                                             cTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGaa                                                                                                   >  3:23027/1‑102 (MQ=255)                                                                              ttCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGtt                                                   >  6:136743/1‑149 (MQ=255)                                                                                       gAACTGTACGTCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCt                                           <  1:49149/148‑1 (MQ=255)                                                                                       gAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTg                                          <  5:134932/149‑1 (MQ=255)                                                                                         aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCg                                        <  1:42489/149‑1 (MQ=255)                                                                                         aCTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGAGCTGGCTACGGTGATGAAGTCGGACGGAATAGCTGTAACTTTATCGCCGTGGCTCATCCAGAAATCGAGGAGCGGTTTACCGTCTGCg                                        <  4:134957/149‑1 (MQ=255)                                                                                               cGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAg                                   <  5:136743/148‑1 (MQ=255)                                                                                                  cATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                             <  7:135311/55‑1 (MQ=255)                                                                                                  cATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCt                                                                                                                             >  8:135311/1‑55 (MQ=255)                                                                                                      gAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCtt                           <  2:162714/149‑1 (MQ=255)                                                                                                       aaGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTc                          >  2:84926/1‑149 (MQ=255)                                                                                                                   ccGTTTGCCATTATGGCTAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAAc              >  5:179295/2‑149 (MQ=255)                                                                                                                         gCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTgc        >  7:96843/1‑149 (MQ=255)                                                                                                                             tAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGGGCATCTTCGATACCGCGAACCAGTGCGCTg    >  5:140348/1‑149 (MQ=255)                                                                                                                                tGGCTAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTc  <  6:179295/148‑1 (MQ=255)                                                                                                                                tGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGa                        <  7:113719/126‑1 (MQ=255)                                                                                                                                tGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGa                        >  8:113719/1‑126 (MQ=255)                                                                                                                                         |                                                                                                                                           GTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTGCGCTGTC  >  NC_000913/2631910‑2632184

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