Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,329 A→G 16.7% 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_002947812,3290AG16.7% 79.7 / 9.1 36intergenic (‑78/‑150)PP_0698/PP_0699LysR family transcriptional regulator/LysE family transporter
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base G (4/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

GCGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCTCA  >  NC_002947/812181‑812476
                                                                                                                                                    |                                                                                                                                                   
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:52517/1‑148 (MQ=21)
                                                                                                        gggTAAACCTGCCCCTACAGGGAAATATGTGGCCATGTAGGAGCGGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGa                                              >  6:148687/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)
                                                                                                                                                   cAGGTTTACCCTCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCTCa  <  3:42651/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCGCACGAAGCTCTCGGTCTGGGCGGCAATGTCGAAGTAGCGGAAGGCGTTGAGTGGCGGCAGTTTCATGATATCCCTGTTCAGCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGCCACATATTTCCCTGTAGGGGCAGGTTTACCCGCGAAGAGGCCAAGGCAGGCTGTACATATCTGTGAGCAAAACTATCAGATTCATCTCAGCATAAATCGATTTTTCCGCATCTGTAAGCGCTTGATAATCGGCGGTCTGTGCATTTCCATCGCCTATCGAGAGCCCTCA  >  NC_002947/812181‑812476

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: