Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,092,790 C→A 11.3% A294S (GCC→TCC)  PP_4482 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,092,7900CA11.3% 121.9 / 6.0 44A294S (GCC→TCC) PP_4482AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base A (2/3);  total (22/22)
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.64e-01

AACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTC  >  NC_002947/5092647‑5092923
                                                                                                                                               |                                                                                                                                     
aaCACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCCCTCCAAACCGCCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTggg                                                                                                                                 >  1:688903/1‑149 (MQ=255)
       aCCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGc                                                                                                                          <  1:1034086/150‑1 (MQ=255)
        ccccAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCg                                                                                                                         <  2:1125976/150‑1 (MQ=255)
                     acCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTTTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGaa                                                                                                            <  1:305953/150‑1 (MQ=255)
                       cccGACCCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                                                                                          <  2:1073267/150‑1 (MQ=255)
                       cccGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                                                                                          >  1:367046/1‑150 (MQ=255)
                       cccGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                                                                                          <  1:723922/150‑1 (MQ=255)
                          gAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCcac                                                                                                       >  1:734699/1‑150 (MQ=255)
                           aaCCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCcaca                                                                                                      <  1:737545/150‑1 (MQ=255)
                           aaCCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCcaca                                                                                                      <  1:645698/150‑1 (MQ=255)
                              ccacccGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAAGAGAAGCAACAgag                                                                                                   >  2:283957/1‑148 (MQ=255)
                                acccGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACa                                                                                                 >  1:769348/1‑150 (MQ=255)
                                  ccGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATAGCGGTCTTCCCGAGGCGTTGCACCAAAGAAATTACAAAAAAAAAAAAAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGc                                                                                               <  2:244033/150‑1 (MQ=255)
                                   cGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                                                                                          >  1:457550/1‑138 (MQ=255)
                                   cGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGc                                                                                                          <  2:457550/138‑1 (MQ=255)
                                              cccTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGgatg                                                                                   >  2:768235/1‑150 (MQ=255)
                                                    tCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGAc                                                                             <  2:249443/150‑1 (MQ=255)
                                                       gCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTg                                                                          >  1:881661/1‑150 (MQ=255)
                                                            cgcCGAGCTCAATTCAACCGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGt                                                                     <  2:96692/150‑1 (MQ=255)
                                                                         tCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACt                                                                            >  1:111614/1‑130 (MQ=255)
                                                                         tCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACt                                                                            <  2:111614/130‑1 (MQ=255)
                                                                                  ctgctgctgCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGcc                                               <  1:768235/150‑1 (MQ=255)
                                                                                        ctgcGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGAACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGt                                                                     <  1:440113/122‑1 (MQ=255)
                                                                                        ctgcGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGAACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGt                                                                     <  1:1015738/122‑1 (MQ=255)
                                                                                        ctgcGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGAACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGt                                                                     >  2:440113/1‑122 (MQ=255)
                                                                                        ctgcGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGAACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGt                                                                     >  2:1015738/1‑122 (MQ=255)
                                                                                          gcGACGCTGATTGCGGTCTTCCCGCGGCGCTACACCCAGCAAGTTGTGACAGTCTCTGTAAAAATGCGGCCTCGAGGACAAGTGATATTATTGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCAAGCGCCCGTTGTTCaa                                       >  1:809908/1‑150 (MQ=255)
                                                                                           cGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAAc                                      <  2:717983/150‑1 (MQ=255)
                                                                                                 tGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAAc                                <  2:634301/150‑1 (MQ=255)
                                                                                                      gCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAg                           <  1:902983/150‑1 (MQ=255)
                                                                                                      gCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAg                           >  1:601118/1‑150 (MQ=255)
                                                                                                      gCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAg                           >  1:601210/1‑150 (MQ=255)
                                                                                                            ttCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGgtagta                     >  2:559146/1‑150 (MQ=255)
                                                                                                                cgcgGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTgg                                                                                      <  1:903371/81‑1 (MQ=255)
                                                                                                                cgcgGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTgg                                                                                      >  2:903371/1‑81 (MQ=255)
                                                                                                                    gcgTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTActggctgg             <  2:388667/150‑1 (MQ=255)
                                                                                                                      gTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGc           >  2:82927/1‑150 (MQ=255)
                                                                                                                         gCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCTCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCAcg        <  2:688903/150‑1 (MQ=255)
                                                                                                                           aCCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCa                                                                <  2:728404/92‑1 (MQ=255)
                                                                                                                           aCCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCa                                                                >  1:728404/1‑92 (MQ=255)
                                                                                                                            ccAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATAAGCATCTGCCGCGCCTTGTTCAACCTCAACTCCAGGTAGTACTGGCTGGGCACGCgg     >  1:478202/1‑150 (MQ=255)
                                                                                                                             cAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACTGCCCGATCTAACTGACCGCAGCAGAGGGCTGCATCAGCATCTCCCGCGCCTTGTTCAACCGCGACTCCAGGTAGCACTGGCGGGGCACGCggg    >  2:587699/1‑149 (MQ=255)
                                                                                                                              aaaGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCAACAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGtt   >  2:898586/1‑150 (MQ=255)
                                                                                                                               aaGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTc  >  2:735805/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                     
AACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAAAGAAATTACGGTAGGCACTGGAAAAATGCGGCCCCGAGGAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTC  >  NC_002947/5092647‑5092923

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