Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,544,963	G→T	15.9%	G281G (GGC→GGA)	hyuC ←	bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,544,963	0	G	T	15.9%	115.8 / 13.9	44	G281G (GGC→GGA)	hyuC	bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base T (5/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATGGG  >  NC_002947/4544818‑4545105                                                                                                                                                   |                                                                                                                                               gCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGccc                                                                                                                                             <  3:135504/149‑1 (MQ=255)   aCCAACTCATTGCTCAGGCCTTTCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCATGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCgg                                                                                                                                           <  6:188371/149‑1 (MQ=255)   aCCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCgg                                                                                                                                           <  8:52700/149‑1 (MQ=255)     cAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGt                                                                                                                                         <  5:169521/149‑1 (MQ=255)         tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgggcgaag                                                                                                                                     >  8:44318/1‑144 (MQ=255)         tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCACCCGGGC‑TCGCGTAACGTGATCCCgggc                                                                                                                                         <  8:222003/145‑2 (MQ=255)         tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCACCCGGGC‑TCGCGTAACGTGATCCCgggc                                                                                                                                         >  7:222003/1‑144 (MQ=255)               cTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGTGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAg                                                                                                                               <  5:114030/149‑1 (MQ=255)                tCAGGCCATGCTTGGCACAGGACGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCCCGTTACGCGAGCCCGGGTACGCCTGAAgg                                                                                                                              >  8:131682/1‑149 (MQ=255)                     ccATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAAcc                                                                                                                         <  6:250107/149‑1 (MQ=255)                             ggCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGc                                                                                                                 <  8:171354/149‑1 (MQ=255)                                       gCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTg                                                                                                       <  1:35373/149‑1 (MQ=255)                                        ccTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTgg                                                                                                      <  1:253228/149‑1 (MQ=255)                                         cTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTggg                                                                                                     >  2:267003/1‑149 (MQ=255)                                                aCCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCtggtg                                                                                              <  2:92768/149‑1 (MQ=255)                                                aCCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCtggtg                                                                                              <  1:267003/149‑1 (MQ=255)                                                              ggCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGc                                                                                >  3:129254/1‑149 (MQ=255)                                                                  atcatcGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGtt                                                                            <  2:152628/149‑1 (MQ=255)                                                                          aTTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCct                                                                    >  1:249136/1‑149 (MQ=255)                                                                            tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgg                                                                                                                                           <  8:71002/76‑1 (MQ=37)                                                                            tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgg                                                                                                                                           >  7:71002/1‑76 (MQ=37)                                                                            tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCctc                                                                   >  8:104245/1‑148 (MQ=255)                                                                              aaCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                                >  8:122119/1‑149 (MQ=255)                                                                              aaCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                                >  4:248842/1‑149 (MQ=255)                                                                               aCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                                >  7:106647/1‑148 (MQ=255)                                                                                cTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACca                                                              <  3:248842/149‑1 (MQ=255)                                                                                   cTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGCTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAg                                                           >  7:43284/1‑149 (MQ=255)                                                                                                        gTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGt                                      >  8:193794/1‑149 (MQ=255)                                                                                                         tCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTc                                     <  1:30501/149‑1 (MQ=255)                                                                                                           gAGGGTCATGCGCACTTCGCCCGGGAGTCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTcc                                    >  1:2580/1‑149 (MQ=255)                                                                                                              ggTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTgc                                <  8:106647/149‑1 (MQ=255)                                                                                                                     cgcACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGt                          <  3:115156/148‑1 (MQ=255)                                                                                                                        aCTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGccccc                                                >  4:33748/1‑123 (MQ=255)                                                                                                                        aCTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGccccc                                                <  3:33748/123‑1 (MQ=255)                                                                                                                          ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAg                            >  8:13262/1‑141 (MQ=255)                                                                                                                          ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAg                            <  7:13262/141‑1 (MQ=255)                                                                                                                          ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACAGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCg                    >  6:112081/1‑149 (MQ=255)                                                                                                                           gtacccGGGC‑TCGCGTAACGTGATCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCgg                   >  1:76/4‑149 (MQ=255)                                                                                                                             gCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCg                 >  5:65151/1‑149 (MQ=255)                                                                                                                                        gcGTAACGTGATCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCgg        >  2:8360/2‑148 (MQ=255)                                                                                                                                           ttACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATg    <  2:249136/149‑1 (MQ=255)                                                                                                                                              cgcgAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACAGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATggg  <  5:112081/148‑1 (MQ=255)                                                                                                                                               gcgAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGc                                        >  2:241688/1‑109 (MQ=255)                                                                                                                                               gcgAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGc                                        <  1:241688/109‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATGGG  >  NC_002947/4544818‑4545105

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