Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,336,737:1	+T	11.1%	coding (187/1443 nt)	PP_4696 →	nitrogen regulation protein NR(I)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,336,737	1	.	T	11.1%	124.5 / 9.4	36	coding (187/1443 nt)	PP_4696	nitrogen regulation protein NR(I)
Reads supporting (aligned to +/- strand):  ref base . (14/18);  new base T (2/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CGCTCGGCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGC  >  NC_002947/5336590‑5336873                                                                                                                                                     |                                                                                                                                         cGCTCGGCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑t                                                                                                                                         >  3:61551/1‑149 (MQ=255)       gCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCaa                                                                                                                                   <  3:168408/149‑1 (MQ=255)          cTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCt                                                                                                                                >  6:309950/1‑149 (MQ=255)                  cagggccgAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAgg                                                                                                                        <  1:205385/142‑1 (MQ=255)                             ggAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTg                                                                                                             <  6:56205/149‑1 (MQ=255)                                 ccAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGa                                                                                                         >  3:265160/1‑149 (MQ=255)                                        cAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCa                                                                                                                          >  2:269855/1‑125 (MQ=255)                                           gTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCGCCGAGCTTGATC‑AGCTCggtgccg                                                                                                                        >  1:199424/1‑119 (MQ=255)                                           gTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCGCCGAGCTTGATC‑AGCTCggtgccg                                                                                                                        <  2:199424/124‑6 (MQ=255)                                                              tGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCa                                                                            <  2:272794/149‑1 (MQ=255)                                                                   gAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCt                                                                                                                                <  3:229359/92‑1 (MQ=255)                                                                   gAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCt                                                                                                                                >  4:229359/1‑92 (MQ=255)                                                                    aaGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTg                                                                      <  3:276870/149‑1 (MQ=255)                                                                        ccAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACt                                                                  <  5:309950/149‑1 (MQ=255)                                                                                  ttCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAAt                                                        >  8:121412/1‑149 (MQ=255)                                                                                       cATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGgcg                                                   <  4:265160/149‑1 (MQ=255)                                                                                             cACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTg                                              >  2:82552/1‑148 (MQ=255)                                                                                                  tCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCAGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTa                                         >  8:122717/1‑148 (MQ=255)                                                                                                    gACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACa                                       <  1:5817/148‑1 (MQ=255)                                                                                                          aTCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCaa                                <  7:121412/149‑1 (MQ=255)                                                                                                               cAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCGCCGCGGTGCTGATc                                                                                                       >  7:111608/1‑73 (MQ=255)                                                                                                               cAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCGCCCCGGTGCTGATc                                                                                                       <  8:111608/73‑1 (MQ=255)                                                                                                                aGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCAAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCtt                          >  2:14246/1‑149 (MQ=255)                                                                                                                  cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGc                                                                                     <  2:240737/88‑1 (MQ=255)                                                                                                                  cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGc                                                                                     >  1:240737/1‑88 (MQ=255)                                                                                                                  cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGa                                                               <  7:119614/110‑1 (MQ=255)                                                                                                                  cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGa                                                               >  8:119614/1‑110 (MQ=255)                                                                                                                  cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCg                        >  3:187820/1‑149 (MQ=255)                                                                                                                         cgcCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCacg                  <  1:271928/148‑1 (MQ=255)                                                                                                                               cccGGCGCGCCAGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGAt           <  7:122717/149‑1 (MQ=255)                                                                                                                                  ggCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGAtgct        <  1:14246/149‑1 (MQ=255)                                                                                                                                         ccGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTc                         <  1:203512/125‑1 (MQ=255)                                                                                                                                         ccGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTc                         >  2:203512/1‑125 (MQ=255)                                                                                                                                          cctggccGAGCTTGATC‑AGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTgcaggc  >  2:286522/6‑147 (MQ=255)                                                                                                                                          cctggccGAGCTTGATC‑AGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTgcaggc  <  1:286522/142‑1 (MQ=255)                                                                                                                                             cACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACgaga            <  4:58117/134‑1 (MQ=255)                                                                                                                                             cACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACgaga            >  3:58117/1‑134 (MQ=255)                                                                                                                                                     |                                                                                                                                         CGCTCGGCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGC  >  NC_002947/5336590‑5336873

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: ‑