Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,892,528	G→A	15.0%	intergenic (+38/‑57)	PP_1697 → / → PP_1698	transcriptional regulator, GntR family/MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,892,528	0	G	A	15.0%	94.9 / 8.3	40	intergenic (+38/‑57)	PP_1697/PP_1698	transcriptional regulator, GntR family/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base A (2/4);  total (17/23)
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.69e-01

GTCGTGATCGAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTCGGCGGCTGGTCGCTCGATGCGCTGGAAGTG  >  NC_002947/1892386‑1892674                                                                                                                                               |                                                                                                                                                   gtCGTGATCGAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACgg                                                                                                                                              <  6:261992/149‑1 (MQ=255)    gTGATCGAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGGCAGCACGgcgc                                                                                                                                           <  3:9697/149‑1 (MQ=255)          gAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGGCAGCACGGCGCAgaca                                                                                                                                      <  1:256309/148‑4 (MQ=255)          gAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGc                                                                                                                                     >  2:113883/1‑149 (MQ=255)                  ggCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACg                                                                                                                                               >  4:209170/1‑131 (MQ=255)                  ggCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACg                                                                                                                                               <  3:209170/131‑1 (MQ=255)                   gCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCATTGGCTTCGCCTGc                                                                                                                            <  8:95079/149‑1 (MQ=255)                   gCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGc                                                                                                                            <  4:77630/149‑1 (MQ=255)                         gATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTca                                                                                                                      <  8:141670/149‑1 (MQ=255)                                  ggCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGa                                                                                                             <  1:113883/149‑1 (MQ=255)                                        ccTCTGGAATGGCTCGGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGaa                                                                                                       <  4:190197/149‑1 (MQ=255)                                         ctctGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAg                                                                                                      >  1:198040/1‑149 (MQ=255)                                          tctGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGtgctgc                                                                                                                                                   <  1:154469/103‑1 (MQ=255)                                          tctGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGtgctgc                                                                                                                                                   >  2:154469/1‑103 (MQ=255)                                          tctGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCacac                                                                                                                   <  1:182229/135‑1 (MQ=255)                                          tctGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCacac                                                                                                                   >  2:182229/1‑135 (MQ=255)                                                   gcgcCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACACACGACTAGAAGGATAACAAt                                                                                             >  8:188439/1‑148 (MQ=255)                                                    cgccgcTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCa                                                                                           <  5:77183/149‑1 (MQ=255)                                                     gccgcTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGc                                                                                                                                <  4:242834/111‑1 (MQ=255)                                                     gccgcTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGc                                                                                                                                >  3:242834/1‑111 (MQ=255)                                                      ccgcTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATg                                                                                         >  2:61965/1‑149 (MQ=255)                                                       cgcTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGt                                                                                        <  8:156634/149‑1 (MQ=255)                                                          tCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCAACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGtt                                                                                       >  8:198285/1‑147 (MQ=255)                                                              gcgGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTg                                                                                 >  1:65203/1‑149 (MQ=255)                                                               cgGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTgg                                                                                <  3:97281/149‑1 (MQ=255)                                                                        gcagcaCTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATcgc                                                                          >  8:174093/1‑146 (MQ=255)                                                                        gcagcaCTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCgag                                                                       <  7:198285/149‑1 (MQ=255)                                                                         cagcaCTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCCTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAg                                                                                                                                         <  2:19535/82‑1 (MQ=255)                                                                         cagcaCTTCCGCTACCGCATCGACATTCGATAAACCGGTCGCCTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAg                                                                                                                                         >  1:19535/1‑82 (MQ=255)                                                                           gcaCTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATc                                                                    <  5:144694/149‑1 (MQ=255)                                                                            caCTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCa                                                                   <  2:198040/149‑1 (MQ=255)                                                                               ttCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCt                                                                >  2:239418/1‑149 (MQ=255)                                                                                   gCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACg                                                            >  2:251450/1‑149 (MQ=255)                                                                                    cTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTGGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACgg                                                           >  1:85192/1‑149 (MQ=255)                                                                                       ccGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCgaga                                                                      >  5:118447/1‑135 (MQ=255)                                                                                       ccGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCgaga                                                                      <  6:118447/135‑1 (MQ=255)                                                                                              gagaTTCGTTAATCCGGTCGCTTCACCCTTCACCTGTCTGCGACGTGCTGCAACGGCGCAGTGGCTTCGCCTGCCCGTAACACCCACGAATAGCAGGATAACAATCATGTTCAGCTGGTCGCTCGAGATCACCTCACGGGAACGCAAg                                                  <  5:187837/145‑1 (MQ=255)                                                                                                         aaCCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGc                                       <  4:257847/148‑1 (MQ=255)                                                                                                                                   cTGCGCCGTGGCAGCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTcgg                              <  1:20784/131‑1 (MQ=255)                                                                                                                                   cTGCGCCGTGGCAGCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTcgg                              <  4:103679/131‑1 (MQ=255)                                                                                                                                   cTGCGCCGTGGCAGCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTcgg                              >  3:103679/1‑131 (MQ=255)                                                                                                                                   cTGCGCCGTGGCAGCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTcgg                              >  2:20784/1‑131 (MQ=255)                                                                                                                                              tgcCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTCGGCGGCTGGTCGCTCGATGCGCTGGAAGTg  >  3:205644/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                   GTCGTGATCGAGCGCCTGGCCCGCGATTACGCCGGCAAACCTCTGGAATGGCGCCGCTCGCGCGGCCATGCGCAGCACTTCCGCTACCGCATCGACATTCGTTAAACCGGTCGCTTCACCCTTCACCTGTCTGCGCCGTGCTGCCACGGCGCAGTGGCTTCGCCTGCCCGTCACACCCACGACTAGAAGGATAACAATCATGTTCAGCTGGTATCGCGAGATCACCTCACGGGAACGCAAGACATTCTGGGCCTGTTTCGGCGGCTGGTCGCTCGATGCGCTGGAAGTG  >  NC_002947/1892386‑1892674

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