Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	812,329	A→G	9.6%	intergenic (‑78/‑150)	PP_0698 ← / → PP_0699	LysR family transcriptional regulator/LysE family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	812,329	0	A	G	9.6%	144.8 / 5.3	52	intergenic (‑78/‑150)	PP_0698/PP_0699	LysR family transcriptional regulator/LysE family transporter
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (2/3);  total (24/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 = 6.86e-01

CGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCT  >  NC_002947/812182‑812474                                                                                                                                                       |                                                                                                                                                  cgcgcACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAg                                                                                                                                                  <  7:292546/149‑1 (MQ=255)  gcgcACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAgg                                                                                                                                                 <  4:251998/149‑1 (MQ=255)     cACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGttt                                                                                                                                              <  2:133762/149‑1 (MQ=255)       cGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAc                                                                                                                                            <  8:39294/149‑1 (MQ=255)       cGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAc                                                                                                                                            <  4:163071/149‑1 (MQ=255)        gAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAcc                                                                                                                                           >  4:120400/1‑149 (MQ=255)          aGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCCTGCCTT‑‑GGCCTCTTCGCGGGTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTAccc                                                                                                                                          <  5:226497/149‑1 (MQ=14)           gCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCCTTGGCCTCTTCGCGGGTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTAccc                                                                                                                                          <  7:276700/149‑1 (MQ=17)            ctctCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCcgc                                                                                                                                        >  4:29174/1‑148 (MQ=255)                cGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAaga                                                                                                                                   <  8:226198/149‑1 (MQ=255)                  gTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGttt                                                                                                                                              >  3:117912/1‑136 (MQ=18)                  gTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGttt                                                                                                                                              <  4:117912/136‑1 (MQ=18)                  gTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAgg                                                                                                                                 >  3:225690/1‑149 (MQ=255)                         cggcAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATCATATCCCTGTTCAG‑CTGGGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAaggc                                                                                                                          <  5:52943/149‑1 (MQ=255)                           gcAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCGCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAaggcag                                                                                                                        <  7:249776/149‑1 (MQ=255)                              aTGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCt                                                                                                                     <  5:78098/149‑1 (MQ=255)                                gTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCCGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACACGCGAAGAGGCCAAGGCAGGCTg                                                                                                                    >  7:161645/1‑148 (MQ=255)                                     aGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACata                                                                                                              >  7:359390/1‑149 (MQ=255)                                       tAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATc                                                                                                            >  7:83132/1‑149 (MQ=255)                                          cGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCtgt                                                                                                         <  5:145253/149‑1 (MQ=255)                                          cGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCtgt                                                                                                         <  3:120400/149‑1 (MQ=255)                                                gCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCa                                                                                                    >  7:117529/1‑148 (MQ=255)                                                     gAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTa                                                                                              >  7:67997/1‑149 (MQ=255)                                                       gTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATc                                                                                            <  8:100418/149‑1 (MQ=255)                                                          gcggcAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTATTCGCGGGTAAACCCGCGCCTACATGGGCACATATTTACCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGGCGGGTGTACATATCTGTGAGCAAAACTAt                                                                                             >  4:379408/1‑145 (MQ=38)                                                                    cATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAGCCGCGAAGAGGCCAAGGCAGGCTGTAGATATCTGGGAGGAAAACTTTCAGATTCATCTCAg                                                                               >  6:200450/1‑149 (MQ=38)                                                                    cATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCa                                                                                <  8:83132/148‑1 (MQ=255)                                                                    cATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCa                                                                                <  4:225690/148‑1 (MQ=255)                                                                           tCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATaaa                                                                         <  3:306647/148‑1 (MQ=255)                                                                             ccTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCa                                                                                           >  5:354193/1‑128 (MQ=255)                                                                             ccTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCa                                                                                           <  6:354193/128‑1 (MQ=255)                                                                             ccTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCg                                                                      >  3:307118/1‑149 (MQ=255)                                                                                 ttCAG‑CTGTGCC‑‑CGCCGCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTACCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGAt                                                                    <  6:181026/147‑1 (MQ=255)                                                                                       cTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCg                                                             <  3:379408/149‑1 (MQ=255)                                                                                                  ccTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTaa                                                    <  4:44601/149‑1 (MQ=255)                                                                                                    tctTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAgc                                                  >  2:381956/1‑149 (MQ=255)                                                                                                        cgcgGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAgcgc                                                <  8:359390/147‑1 (MQ=255)                                                                                                        cgcgGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCtt                                              >  4:130243/1‑149 (MQ=255)                                                                                                         gggggTAAACCCTCTCCTACATGGCCTCATATTTCCCTGTAGGTGCAGGTTTACACGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTg                                             <  8:161645/147‑1 (MQ=255)                                                                                                              tAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGGTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAACCTATCAGATGCAGCTCAGCATAAATCGATTTTTCCGCATCTGGTAGCGCTTgatgat                                        >  6:317502/1‑149 (MQ=17)                                                                                                                 aCCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATcgg                                     <  3:29174/149‑1 (MQ=255)                                                                                                                   ccGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATcggcg                                   >  4:327284/1‑149 (MQ=255)                                                                                                                   ccGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATcggcg                                   >  6:53764/1‑149 (MQ=255)                                                                                                                     gCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGt                                 >  7:211209/1‑149 (MQ=255)                                                                                                                         cTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCtgt                             >  3:307645/1‑149 (MQ=255)                                                                                                                                 ccACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTc                      >  7:227190/1‑148 (MQ=255)                                                                                                                                     atatTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCg                 >  7:19671/1‑149 (MQ=255)                                                                                                                                        tttCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATc                  >  7:79313/1‑145 (MQ=255)                                                                                                                                        tttCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATc                  <  8:79313/145‑1 (MQ=255)                                                                                                                                             cTGTAGGCGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCCGTCTGTGCATTTCCATCGCCTATc           >  7:153993/1‑147 (MQ=255)                                                                                                                                                tAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGagag      <  4:307118/149‑1 (MQ=255)                                                                                                                                                    ggCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCt  <  3:327284/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAG‑CTGTGCC‑‑GGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCT  >  NC_002947/812182‑812474

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