Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,273,070	C→G	15.6%	intergenic (+35/+256)	PP_2003 → / ← PP_2004	hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,273,070	0	C	G	15.6%	101.5 / 14.2	45	intergenic (+35/+256)	PP_2003/PP_2004	hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (17/20);  new base G (3/4);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

CGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCA  >  NC_002947/2272923‑2273217                                                                                                                                                     |                                                                                                                                                    cGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcg                                                                                                                                                    <  7:90195/149‑1 (MQ=255)           tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggtc                                                                                                                                          <  2:312662/149‑3 (MQ=255)             aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCgg                                                                                                                                        <  3:37208/149‑1 (MQ=255)             aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCgg                                                                                                                                        <  5:316043/149‑1 (MQ=255)             aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTAGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                                         <  5:373162/148‑1 (MQ=38)               aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCATTAGCCGGCAGCCCGCGCAACGGCCGGt                                                                                                                                       >  2:167159/1‑148 (MQ=255)                   aGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAg                                                                                                                                   >  4:311054/1‑148 (MQ=38)                   aGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAg                                                                                                                                   <  3:311054/148‑1 (MQ=38)                   aGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAgg                                                                                                                                  <  8:279367/149‑1 (MQ=255)                         cGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTa                                                                                                                             >  8:217323/1‑148 (MQ=255)                            aGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACa                                                                                                                         >  6:93679/1‑149 (MQ=255)                                  tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATcccc                                                                                                                    <  5:26573/148‑1 (MQ=255)                                                    cTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTGCAGACGCCCCGCGa                                                                                                 >  8:250594/1‑149 (MQ=255)                                                     tGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGa                                                                                                 >  1:260084/1‑149 (MQ=255)                                                                aGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTGCAGACGCCCCGCGATAACTTTCTGaa                                                                                     <  7:250594/149‑1 (MQ=255)                                                                   tGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAacac                                                                                  <  4:294189/149‑1 (MQ=255)                                                                       gcagGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGc                                                                              <  7:217323/149‑1 (MQ=255)                                                                                ccAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGATAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                     >  4:71446/1‑149 (MQ=255)                                                                                ccAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                     <  1:162672/149‑1 (MQ=255)                                                                                ccAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                     <  1:335865/149‑1 (MQ=255)                                                                                ccAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCATTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACCCCCCGCGATAACTTTCTGAACACCAGCGACGTTTAc                                                                     >  4:314231/1‑149 (MQ=255)                                                                                            ggCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCt                                                         <  2:260084/149‑1 (MQ=255)                                                                                             gcc‑gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                                        <  1:167159/149‑1 (MQ=255)                                                                                               c‑gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCGTTAGCCGGCAGCCCGCGAAAGGGCCGGGTCAGGGTGTAAACATCCCCAGGTTCAGACGCCCCGCGAGACCTTTGTGAACACCTGGGGCGGTTTCCACTCAGAAGCGTg                                                       >  2:298101/1‑148 (MQ=255)                                                                                                 gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCAc                                                                  >  5:305943/1‑136 (MQ=255)                                                                                                 gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCAc                                                                  <  6:305943/136‑1 (MQ=255)                                                                                                  ccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtca                                                    <  5:93679/149‑1 (MQ=255)                                                                                                     cccccTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGTAAGGGCCGGTTCAGGCTGTAAAGATCCCCTGTTTCAGTCGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatca                                                 <  3:188302/149‑1 (MQ=255)                                                                                                         ctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATccc                                             >  5:178261/1‑149 (MQ=255)                                                                                                           ccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                           >  1:196235/1‑149 (MQ=255)                                                                                                             ctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                                        <  7:314932/102‑1 (MQ=255)                                                                                                             ctTTGATGGGTCGACTGAACCGACCCTTTCGCCGGCAGCAAGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTTCAGACGCCCCGCGATAACTTTCAGAACCACGGATACGTTTACCACTCATAAGATTCTCATCATCCCTTc                                          >  2:235664/1‑148 (MQ=255)                                                                                                              tttGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCAt                                        >  7:95371/1‑149 (MQ=255)                                                                                                                      gTCGCCTGAACCGCCCCATTAGCCGGCCGCCCGCGAAAGGGCCGGTTCAGGCTGTTAACATACCCAGTTAAAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTTCCACGCATAAGCTTGtaa                                                    >  4:77563/1‑127 (MQ=255)                                                                                                                      gTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtca                                                    <  3:77563/129‑1 (MQ=255)                                                                                                                          ccTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCt                              >  7:339620/1‑147 (MQ=255)                                                                                                                          ccTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa                            <  6:178261/149‑1 (MQ=255)                                                                                                                           cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg                           <  3:338313/149‑1 (MQ=255)                                                                                                                             gAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGc                         >  7:12886/1‑149 (MQ=255)                                                                                                                                 cggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTg                       >  8:298952/4‑147 (MQ=39)                                                                                                                                   aCCCTTTCGCCGGCAGCCCGCGATAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAgg                   <  3:71446/149‑1 (MQ=255)                                                                                                                                      cTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa                <  8:53816/149‑1 (MQ=255)                                                                                                                                      cTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa                >  7:53816/1‑149 (MQ=255)                                                                                                                                      cTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATAATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa                >  8:124391/1‑149 (MQ=255)                                                                                                                                           gCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCt           <  2:196235/149‑1 (MQ=255)                                                                                                                                             cGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACACCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACGCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTa          >  3:291601/1‑148 (MQ=255)                                                                                                                                                    ccGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGc   >  3:149661/1‑148 (MQ=255)                                                                                                                                                     cgcgAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCa  >  5:1488/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCC‑GCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGCA  >  NC_002947/2272923‑2273217

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