Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,326 G→A 17.6% intergenic (‑75/‑153) 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_002947812,3260GA17.6% 74.5 / 9.0 34intergenic (‑75/‑153)PP_0698/PP_0699LysR family transcriptional regulator/LysE family transporter
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base A (4/2);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GCGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCT  >  NC_002947/812181‑812474
                                                                                                                                                 |                                                                                                                                                    
gcgcgcACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGATTGGCGGCATTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCa                                                                                                                                                   <  5:250744/149‑1 (MQ=255)
        gAAGCTCTCGCTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTATTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAcc                                                                                                                                           <  5:228260/149‑1 (MQ=255)
         aaGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTAccc                                                                                                                                          <  6:27207/149‑1 (MQ=255)
         aaGCGCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTAccc                                                                                                                                          <  4:140910/149‑1 (MQ=18)
                         cggcAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAgg                                                                                                                                 <  5:23028/142‑1 (MQ=255)
                                 tCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCAAATATTTCCCTGTAGGGGCAGGTTTAACCGCGAAGAGGCCAAGGCAGGCTGt                                                                                                                   >  7:71961/1‑148 (MQ=255)
                                   gAAGTAGCGGAAGGCGTTGAGTGTCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTAc                                                                                                                 <  4:125413/148‑1 (MQ=255)
                                          cGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTACCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCtgt                                                                                                         >  4:147810/1‑149 (MQ=255)
                                                                           tCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAAt                                                                        <  3:147810/149‑1 (MQ=255)
                                                                            cccTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATc                                                                       >  5:200951/1‑149 (MQ=255)
                                                                                    aGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTc                                                               >  3:200004/1‑149 (MQ=255)
                                                                                     gCTGTGCCGGCCTCTTCGCGGGTTAACAGGCTCCTACTTGGCCACATATTTCCCTGTTGGGGCAGGTTTACCCGCTAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTcc                                                              <  4:93981/149‑1 (MQ=255)
                                                                                      cTGTGCCGGCCTGTTCGCGGGTATACCCTCTCCTACATGGCCTCATTTTTCCCTGTAGGGGCAGGTTTACCAGCGAAGAGGCCAAGGCAGGCTGTAGATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCg                                                             <  3:40456/149‑1 (MQ=255)
                                                                                         tgCCGGCCTGTTCGCGGGTAACCCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCAt                                                          <  7:90020/149‑1 (MQ=255)
                                                                                           ccGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGAGTTTTCCGCATCt                                                        >  7:172164/1‑149 (MQ=255)
                                                                                           ccGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTACCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTAAGGATAAATCGATTTTTCCGCATCt                                                        >  7:216302/1‑149 (MQ=255)
                                                                                              gCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGt                                                      >  7:165658/1‑148 (MQ=255)
                                                                                                   ttCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAgcgc                                                >  6:69895/1‑149 (MQ=255)
                                                                                                   ttCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAgc                                                  <  8:71961/147‑1 (MQ=255)
                                                                                                     cgcgGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCt                                               <  8:165658/148‑1 (MQ=255)
                                                                                                        gggTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                            >  5:168970/1‑148 (MQ=21)
                                                                                                        gggTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                            >  6:148687/1‑148 (MQ=21)
                                                                                                        gggTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                            >  6:52517/1‑148 (MQ=21)
                                                                                                               cccGCTGGTACATGGCCAGATAGTTCCGTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAGCTATCAGATTCATCTCAGCATTAATCGATTTGTCGGCATCTGTATGCGCTTGATATTcggg                                    >  2:168412/1‑148 (MQ=255)
                                                                                                                ccGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGAATCTGTAAGCGCTTGATAATcggcg                                   >  8:89412/1‑149 (MQ=255)
                                                                                                                      cTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCtgt                             >  7:115296/1‑149 (MQ=255)
                                                                                                                      cTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCt                               <  5:69895/147‑1 (MQ=255)
                                                                                                                        aCATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATGTGTATGCGCTTGATAATCGGCGGTCTGTGc                           >  2:72332/1‑149 (MQ=255)
                                                                                                                              ccACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTc                      >  6:627/1‑148 (MQ=255)
                                                                                                                                       gccatGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                            <  2:207564/113‑1 (MQ=255)
                                                                                                                                       gccatGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                            >  1:207564/5‑117 (MQ=255)
                                                                                                                                            gTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTc                                                               <  4:24288/93‑1 (MQ=255)
                                                                                                                                            gTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTc                                                               >  3:24288/1‑93 (MQ=255)
                                                                                                                                             tAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCgaga       <  8:115296/148‑1 (MQ=255)
                                                                                                                                                gggCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGccc   >  5:265580/1‑149 (MQ=255)
                                                                                                                                                 ggCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATAGCCTATCGAGAGCCCt  >  5:103881/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCT  >  NC_002947/812181‑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: