Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,076,041	A→G	8.7%	V183A (GTC→GCC)	phnN ←	ribose 1,5‑bisphosphate phosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,076,041	0	A	G	8.7%	124.8 / 2.9	46	V183A (GTC→GCC)	phnN	ribose 1,5‑bisphosphate phosphokinase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (2/2);  total (24/22)
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.97e-01

CGCTGCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAGCAAGCGCTGGCGCAGCACCTCA  >  NC_002947/5075897‑5076187                                                                                                                                                 |                                                                                                                                                   cgcTGCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGAccgcc                                                                                                                                               >  1:162279/1‑150 (MQ=255)  gcTGCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTATTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGAccgccg                                                                                                                                              >  1:268012/1‑150 (MQ=255)    tGCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATTCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCgccgcc                                                                                                                                            >  2:423346/1‑150 (MQ=255)     gCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCt                                                                                                                                           >  1:594883/1‑150 (MQ=255)        caGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGGTcgccgcatt                                                                                                                                          <  2:371861/148‑4 (MQ=255)        caGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGGTcgccgcatt                                                                                                                                          >  1:371861/1‑145 (MQ=255)                tGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATTCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGgaa                                                                                                                                <  1:423346/150‑1 (MQ=255)                        aTGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAgtgt                                                                                                                                    <  1:158449/138‑1 (MQ=255)                        aTGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAgtgt                                                                                                                                    >  2:158449/1‑138 (MQ=255)                                       cTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGCAACTAATGCGGCGACCGCCGCCTCCAgtgt                                                                                                                                    <  1:267632/123‑1 (MQ=255)                                       cTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAgtgt                                                                                                                                    >  2:267632/1‑123 (MQ=255)                                            gTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGt                                                                                                    >  2:146155/1‑150 (MQ=255)                                                   aTGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGc                                                                                             >  1:450958/1‑150 (MQ=255)                                                    tGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCa                                                                                            <  2:450958/150‑1 (MQ=255)                                                       cTGAACTGCGCTTCTGCGCGTAAGCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGAGAACCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCt                                                                                         <  1:39050/150‑1 (MQ=255)                                                        tGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCtt                                                                                        >  2:514995/1‑150 (MQ=255)                                                             tGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGt                                                                                                    <  2:561859/133‑1 (MQ=255)                                                             tGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGt                                                                                                    >  1:561859/1‑133 (MQ=255)                                                                gCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGt                                                                                <  2:162279/150‑1 (MQ=255)                                                                           tAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGgaagaa                                                                                                                             <  2:484602/94‑1 (MQ=255)                                                                           tAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGgaagaa                                                                                                                             >  1:484602/1‑94 (MQ=255)                                                                            aaTCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGcc                                                                    >  2:330403/1‑150 (MQ=255)                                                                                cTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGgc                                                                >  1:491673/1‑150 (MQ=255)                                                                                 ttttCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGTTCCCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGgcg                                                               <  1:340663/150‑1 (MQ=255)                                                                                 ttttCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGGTcgccgcatt                                                                                                                                          <  2:365828/75‑4 (MQ=255)                                                                                 ttttCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGGTcgccgcatt                                                                                                                                          >  1:365828/1‑72 (MQ=255)                                                                                 ttttCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGACCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGgcg                                                               <  1:340707/150‑1 (MQ=255)                                                                                     cTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGt                                                           <  1:330403/150‑1 (MQ=255)                                                                                            gctATCCCCGCGATAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATAAACCGACGGGTCGGCTGCAGCGTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCActt                                                    >  1:256445/1‑150 (MQ=255)                                                                                                 ccccGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTtcttct                                               <  1:514995/150‑1 (MQ=255)                                                                                                  cccGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTg                                              >  1:524930/1‑150 (MQ=255)                                                                                                  cccGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGATTGCAAGCATGCGCTGCGCGACGGGCTTTGGTCCACGTCTTCTg                                              >  2:133913/1‑150 (MQ=255)                                                                                                   ccGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCa                                                                                                                      >  2:297459/1‑77 (MQ=255)                                                                                                   ccGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCa                                                                                                                      <  1:297459/77‑1 (MQ=255)                                                                                                    cgcgGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTgcg                                            >  1:2671/1‑150 (MQ=255)                                                                                                           gtaCGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTccc                                     <  2:268012/150‑1 (MQ=255)                                                                                                           gtaCGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTccc                                     <  1:146155/150‑1 (MQ=255)                                                                                                              cGCCATGGTCGCGAAGCAGCGAAGCGAATGGCGCGACCGGCGCCTTCAGTGTGGAAGACTTATCCAGCACATGCACCGACGGGTCGGCTCCAGCTTGCAAGCCTGCGCTGCGCCCGGGGCGCTGTTCCACTTCTTCGGCGCTTTCCCGcc                                  <  2:186263/150‑1 (MQ=255)                                                                                                                        gcgAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGGCGGGTCGGCTGCAGCTTGCAAGCGTACGCTACGCGCGACGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAg                        >  2:479254/1‑150 (MQ=255)                                                                                                                                  gAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGACTGCAGCTTGCAAGCGTGCGCTGCGCGCCa                                                                   <  2:431977/97‑1 (MQ=255)                                                                                                                                  gAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGACTGCAGCTTGCAAGCGTGCGCTGCGCGCCa                                                                   >  1:431977/1‑97 (MQ=255)                                                                                                                                     aCGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAGCAAGCGCTGgcgc           >  2:183574/1‑150 (MQ=255)                                                                                                                                      cGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCATGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAGCAAGCGCTGGCgca          <  2:594883/150‑1 (MQ=255)                                                                                                                                       gAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCggg                            <  2:601678/131‑1 (MQ=255)                                                                                                                                       gAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCggg                            >  1:601678/1‑131 (MQ=255)                                                                                                                                              gcgACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAGCAAGCGCTGGCGCAGCACCTCa  <  2:491673/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CGCTGCACAGCGGGATGGCCAGCATGTTAACCAGCATGCTGGCGTTTGTCATGCCTGAACTGGGCTTTTTCGCGTAATCCTTTTCTGCGCTGCTATCCCCGCGGTAGTACGCCCTGATCGCGAAGCAGCGAAACGAATGCGGCGACCGCCGCCTCCAGTGTGGAAGAATTATCCAGCACATGCACCGACGGGTCGGCTGCAGCTTGCAAGCGTGCGCTGCGCGCCAGGCGCTGTTCCACTTCTTCTGCGCTTTCCCGCCCCCGGGCCAGCAAGCGCTGGCGCAGCACCTCA  >  NC_002947/5075897‑5076187

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