Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,352,492	A→G	8.0%	intergenic (‑135/+5)	PP_1178 ← / ← nrdA	hypothetical protein/ribonucleoside‑diphosphate reductase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,352,492	0	A	G	8.0%	137.5 / 2.2	51	intergenic (‑135/+5)	PP_1178/nrdA	hypothetical protein/ribonucleoside‑diphosphate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base G (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CGCGCCATGGGCGTACATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACTTCTCGG  >  NC_002947/1352350‑1352636                                                                                                                                                |                                                                                                                                                          cgcgCCATGGGCGTACATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGATCAGTGGGCAAAGGGCAAGTTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGACACCCGTAGGCACCCCGCACGGCATTCAGGCCt                                                                                                                                                      <  6:21397/148‑1 (MQ=255)  gagcCATGGGCGTACATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGAAAAGGGCAAGGTTTCTCTGAGGCGCGTA‑TACAAAAATGCCGCGCGGGTCACCCGTATGCACCCCGCACCGCATTCAGTCCtt                                                                                                                                                     >  2:314677/3‑148 (MQ=255)               aCATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCg                                                                                                                                       >  3:114684/1‑149 (MQ=255)                cATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGATTCCAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTc                                                                                                                                        <  1:271407/148‑1 (MQ=255)                 aTCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCg                                                                                                                                       >  3:99301/1‑147 (MQ=255)                 aTCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCa                                                                                                                                     <  4:114684/149‑1 (MQ=255)                  tCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACACCGCACGGCATTCAGGCCTTACTGGAAGGCTTCGCAg                                                                                                                                    >  1:224296/1‑149 (MQ=255)                     ttttttGAGAGGTGCTGTTCTTTTACCGAGCTGCGTTAACAGTGGGCAAAGGGCAAGGGTTCTCTGAGTCGAGGA‑TACAAAAATGCCTCGCGGGTCACCCGTAGGCACCCCGCACGGCATTCAGTCCTTACTGGCATGCTTCGCAG‑‑‑TC‑a                                                                                                                             <  1:314677/149‑1 (MQ=255)                       ttttGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCTCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑a                                                                                                                             <  6:286071/147‑1 (MQ=255)                       ttttGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑a                                                                                                                             <  5:234056/147‑1 (MQ=255)                         ttGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCCCCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCt                                                                                                                         >  6:185234/1‑149 (MQ=255)                                    tGTACTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTATCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATcgcg                                                                                                         >  4:135074/1‑149 (MQ=255)                                      ttCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTgg                                                                                                                                                >  2:239735/1‑117 (MQ=255)                                      ttCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTgg                                                                                                                                                <  1:239735/117‑1 (MQ=255)                                      ttCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCa                                                                                                       <  7:276864/149‑1 (MQ=255)                                       tCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCa                                                                                                       <  4:99301/148‑1 (MQ=255)                                        cTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCa                                                                                                                                     >  7:261178/1‑126 (MQ=255)                                        cTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCa                                                                                                                                     <  8:261178/126‑1 (MQ=255)                                         ttttACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCa                                                                                                                                     <  3:4047/125‑1 (MQ=255)                                         ttttACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGc                                                                                                    >  8:227983/1‑149 (MQ=255)                                         ttttACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACAGGCAGGCTTCGCa                                                                                                                                     >  4:4047/1‑125 (MQ=255)                                            tACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGc                                                                                                                          <  8:248369/129‑1 (MQ=255)                                              cctaGCTGCGTGAAGAGTGGTCAAAGGGCATGTTTTGTCTGAGGCGCGGA‑TAAAAAAATGCGGGGCGGGGCACCCGTAGCCACCCCGCACGGAATTCAGGCCTTACTGGCAGGTTTCCCTG‑‑‑TC‑AGGCT‑‑‑‑GGTC‑CTTCGCGCAGGCCttt                                                                                                <  2:224296/145‑1 (MQ=255)                                              ccGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTg                                                                                               <  1:24565/149‑1 (MQ=255)                                               cGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCGGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTgg                                                                                              >  7:195318/1‑149 (MQ=255)                                               cGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCGGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTgg                                                                                              >  7:200563/1‑149 (MQ=255)                                                  gCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCACGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑ATGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCAc                                                                                           <  1:218367/149‑1 (MQ=255)                                                       gTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGCCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGgcg                                                                                      >  7:305161/1‑149 (MQ=255)                                                                           aaGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACt                                                                  >  3:311618/1‑149 (MQ=255)                                                                                     tGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGt                                                        >  4:200727/1‑149 (MQ=255)                                                                                       aGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGt                                                        <  4:311618/147‑1 (MQ=255)                                                                                                  aCAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGgc                                                             >  7:290687/1‑132 (MQ=255)                                                                                                  aCAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGgc                                                             <  8:290687/132‑1 (MQ=255)                                                                                                    aaaaaTGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCAGGACTGGGGCGCTGTCGCCACCGCTGGa                                           >  4:230978/1‑148 (MQ=255)                                                                                                          gCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCg                                    >  7:298618/1‑149 (MQ=255)                                                                                                              cgcgGGGCACCCGTAGGCACCCCGCACGGCATTCGGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGa                                <  8:200563/149‑1 (MQ=255)                                                                                                              cgcgGGGCACCCGTAGGCACCCCGCACGGCATTCGGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGa                                <  8:195318/149‑1 (MQ=255)                                                                                                              cgcgGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGa                                <  7:227983/149‑1 (MQ=255)                                                                                                               gcgGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAg                               <  1:326404/149‑1 (MQ=255)                                                                                                                     cACCCGTAGGCACCCCGCACGGCATTCTGGCCTTACTGGGAGGCTTCGCTGTTTTCTAAGCTCCCCCTTCTTCTCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGtt                                  <  1:279553/149‑1 (MQ=255)                                                                                                                              gCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGgtgt                    >  2:59287/1‑145 (MQ=255)                                                                                                                              gCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGgtgt                    <  1:59287/145‑1 (MQ=255)                                                                                                                               cACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGa                 <  1:284845/147‑1 (MQ=255)                                                                                                                               cACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGa                 >  2:284845/1‑147 (MQ=255)                                                                                                                               cACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATg               <  3:135074/149‑1 (MQ=255)                                                                                                                                      cACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACt        <  5:185234/149‑1 (MQ=255)                                                                                                                                         ggCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACTtct     >  7:312074/1‑149 (MQ=255)                                                                                                                                           cATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACTTCTCg   <  3:93709/149‑1 (MQ=255)                                                                                                                                            aTTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACTTCTCgg  >  5:346280/1‑149 (MQ=255)                                                                                                                                               cAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAAc                                        >  6:175937/1‑108 (MQ=255)                                                                                                                                               cAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAAc                                        <  5:175937/108‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                          CGCGCCATGGGCGTACATCCTTTTTTGAGAGGTCCTGTTCTTTTACCGAGCTGCGTGAACAGTGGGCAAAGGGCAAGGTTTCTCTGAGGCGCGGA‑TACAAAAATGCCGCGCGGGGCACCCGTAGGCACCCCGCACGGCATTCAGGCCTTACTGGCAGGCTTCGCAG‑‑‑TC‑AGGCT‑‑‑‑CGTC‑GATCGCGCAGGCCTTTGGCACTGGCGCAGGGCCGGCTGCCTGGACTGGGGCGCTGTCGCCACCGCTGGAAACGGCGTTGAGCTTGCCGGTGTTGATGGTCGACTTCTCGG  >  NC_002947/1352350‑1352636

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