Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,336,743	Δ1 bp	10.5%	coding (193/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,743	0	A	.	10.5%	112.6 / 6.6	38	coding (193/1443 nt)	PP_4696	nitrogen regulation protein NR(I)
Reads supporting (aligned to +/- strand):  ref base A (16/18);  new base . (2/2);  total (18/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 = 1.00e+00

GCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTAGCGCGTATTCTGC  >  NC_002947/5336596‑5336890                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                                 >  4:229359/1‑92 (MQ=255)                                                             gAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCt                                                                                                                                                 <  3:229359/92‑1 (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                                                                                >  8:119614/1‑110 (MQ=255)                                                                                                            cGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGa                                                                                <  7:119614/110‑1 (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                                          >  2:203512/1‑125 (MQ=255)                                                                                                                                   ccGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTc                                          <  1:203512/125‑1 (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)                                                                                                                                                tGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTAGCGCGTAt       >  4:57141/1‑148 (MQ=255)                                                                                                                                                  aTCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTa               >  1:316357/1‑138 (MQ=255)                                                                                                                                                  aTCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTa               <  2:316357/138‑1 (MQ=255)                                                                                                                                                     aaGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACACCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTAGCGCGTATTCTGc  <  1:82552/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCCTGCGTCGCCTGCTTGAGCGGAACCAGTACCAGGTCAGCGAGGCCGGCTCGGTGCAGGAAGCCCAGGAACGCTTCAGCATTGCCACTTTCGACCTCATCGTCAGCGACCTGCGCCTGCCCGGCGCGCCCGGCACCGAGC‑TGATCAAGCTCGGCCAGGGCACCCCGGTGCTGATCATGACCAGTTACGCCAGCCTGCGCTCAGCGGTGGACTCGATGAAAATGGGCGCGGTGGACTACATCGCCAAGCCCTTCGACCACGACGAGATGCTGCAGGCAGTAGCGCGTATTCTGC  >  NC_002947/5336596‑5336890

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