Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,765 A→T 12.5% H215Q (CAT→CAA)  PP_2004 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,7650AT12.5% 111.9 / 8.2 40H215Q (CAT→CAA) PP_2004AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base T (3/2);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCCAGGTCCGGTTCACGGCCGCCATTGAGCAGCGGGCCCAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTCCAGGCGTT  >  NC_002947/2273628‑2273904
                                                                                                                                          |                                                                                                                                           
tcCAGGTCCGGTTCACGGCCGCCATTGAGCAGCGGGCCCAGCAGGTGGGCGATGCGT‑TTATCCAGACTACTCGCCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTg                                                                                                                                 <  1:492441/150‑1 (MQ=255)
 ccAGGTCCGGTTCACGGCCGCCATTGAGCAGCGGGCCCAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCATTGCCCgagtacct                                                                                                                                <  1:29871/150‑8 (MQ=255)
     gTCCGGTTCACGGCCGCCATTGAGCAGCGGGCCCAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAg                                                                                                                            >  1:322211/1‑150 (MQ=255)
                     caaTTGAGCAGCGGGACCAGCAGGTGGGATATGCGC‑TCACCCATACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCCACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCgg                                                                                                            >  2:772067/3‑150 (MQ=255)
                         tGAGCAGCGGGCCCAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCa                                                                                                        >  2:240318/1‑150 (MQ=255)
                                cGGGCCCAGCAGGTGGGTGATGCGC‑TCCCCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGccg                                                                                                 <  1:114789/150‑3 (MQ=255)
                                    cccAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGa                                                                                             <  1:251871/150‑1 (MQ=255)
                                          gtctgggtgA‑GCGCATCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCg                                                                                       >  1:380405/4‑150 (MQ=255)
                                           tctgggtgA‑GCGCATCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGt                                                                                      >  2:77030/3‑150 (MQ=255)
                                                gtgATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTcggcg                                                                                 >  1:355703/1‑150 (MQ=255)
                                                            aCCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGcac                                                                      >  1:356943/1‑150 (MQ=255)
                                                                 tgcagcgcACGCGGGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGtt                                                                                                                          >  1:140429/5‑93 (MQ=255)
                                                                 tgcagcgcACGCGGGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGtt                                                                                                                          <  2:140429/89‑1 (MQ=255)
                                                                       gcACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGtgctgc                                                           >  2:155477/1‑150 (MQ=255)
                                                                                tGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGtt                                                  <  2:123389/150‑1 (MQ=255)
                                                                                tGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGtt                                                  <  1:797980/150‑1 (MQ=255)
                                                                                     tGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAg                                             <  1:240318/150‑1 (MQ=255)
                                                                                        agcagcTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCATCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGgatgatgccgc                                          >  2:260014/1‑145 (MQ=255)
                                                                                                tGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATa                                  >  2:597942/1‑150 (MQ=255)
                                                                                                  cTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATagg                                <  1:155477/150‑1 (MQ=255)
                                                                                                   ttCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATagga                               <  1:182267/150‑1 (MQ=255)
                                                                                                   ttCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATagga                               >  1:185958/1‑150 (MQ=255)
                                                                                                        acaACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATagga                               <  2:539179/145‑1 (MQ=255)
                                                                                                        acaACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATagga                               >  1:539179/1‑145 (MQ=255)
                                                                                                            cTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGc                      >  2:51551/1‑150 (MQ=255)
                                                                                                             tGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGgcgc                         <  1:751751/146‑1 (MQ=255)
                                                                                                             tGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGgcgc                         >  2:751751/1‑146 (MQ=255)
                                                                                                                aaCAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAAc                  >  1:718545/1‑150 (MQ=255)
                                                                                                                  cAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTc                <  1:772067/150‑1 (MQ=255)
                                                                                                                     gTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGAt             <  2:380405/150‑1 (MQ=255)
                                                                                                                         tgcCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTcc         >  2:507620/1‑150 (MQ=255)
                                                                                                                         tgcCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTcc         <  2:195217/150‑1 (MQ=255)
                                                                                                                          gcCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATa                                              >  2:71196/1‑112 (MQ=255)
                                                                                                                          gcCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATa                                              <  1:71196/112‑1 (MQ=255)
                                                                                                                             aGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTCCAGGc     >  2:388655/1‑150 (MQ=255)
                                                                                                                              caaccccGGGCATTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCgg                       >  2:275709/7‑131 (MQ=255)
                                                                                                                              caaccccGGGCATTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCgg                       <  1:275709/125‑1 (MQ=255)
                                                                                                                               aaccccGGGCATTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGAGCTCCAGGCGt   >  1:362393/6‑150 (MQ=255)
                                                                                                                               aaccccGGGCATTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCGCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTCCAGGCGt   >  2:639042/6‑150 (MQ=255)
                                                                                                                                tACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTCCAGGCGtt  <  2:718545/150‑1 (MQ=255)
                                                                                                                                         aaTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGcgcg                        <  2:53922/119‑1 (MQ=255)
                                                                                                                                         aaTGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGcgcg                        >  1:53922/1‑119 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TCCAGGTCCGGTTCACGGCCGCCATTGAGCAGCGGGCCCAGCAGGTGGGTGATGCGC‑TCACCCAGACTGCGCACCCGCGTGCGCTGCAGCAGCTCTGCTTCGCACAACTGCAACAGGTGCTGCCAGGTACTCGGGCAATGCCCGGGGTTGGCCAGTTGCAAGGTAGCGCGGCTCAGACGCAGCTGATTGCCGTCGGCGGCGAACTGCACAGGCGTGCTGCAAAGCGGTTGATAGCGCGCGGCATAGGAGGGCGCGGCAAACTCGATCTCCAGGCGTT  >  NC_002947/2273628‑2273904

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