Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,093,200 C→T 12.1% S141N (AGC→AAC)  PP_4482 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,093,2000CT12.1% 82.7 / 3.0 34S141N (AGC→AAC) PP_4482AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base T (2/2);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

CCACCAGCTCTTCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACACCCCGGCCGACAGGCCACCGATCAGGC  >  NC_002947/5093069‑5093332
                                                                                                                                   |                                                                                                                                    
ccaccaGCTCTTCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCtt                                                                                                                      <  3:224548/148‑1 (MQ=255)
 caccaGCTCTTCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTgg                                                                                                                    <  7:290135/149‑1 (MQ=255)
        tctTCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCg                                                                                                             >  4:244647/1‑149 (MQ=255)
          ttCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGtt                                                                                                           >  4:200530/1‑149 (MQ=255)
                 acTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGGAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCaaa                                                                                                    <  5:331317/149‑1 (MQ=255)
                 acTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCaaa                                                                                                    <  7:268850/149‑1 (MQ=255)
                   tGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCtggtgg                                                                                                                                <  3:195727/119‑1 (MQ=255)
                   tGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCtggtgg                                                                                                                                >  4:195727/1‑119 (MQ=255)
                     ccccGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTc                                                                                                 <  3:200530/148‑1 (MQ=255)
                               cTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAgc                                                                                      <  6:137396/149‑1 (MQ=255)
                                     cccGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGt                                                                                >  5:324162/1‑149 (MQ=255)
                                          gATCGCGCGCCAGCACAGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGAGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACg                                                                           >  5:41289/1‑149 (MQ=255)
                                             cgcgcgCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCt                                                                        >  7:207083/1‑149 (MQ=255)
                                               cgcgcCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCa                                                                                  >  3:303472/1‑137 (MQ=255)
                                               cgcgcCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCa                                                                                  <  4:303472/137‑1 (MQ=255)
                                                    cAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGc                                                                                          >  8:200034/1‑124 (MQ=255)
                                                    cAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGc                                                                                          <  7:200034/124‑1 (MQ=255)
                                                       cACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTAcc                                                              >  2:213754/1‑149 (MQ=255)
                                                         ccGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCa                                                            <  3:108420/149‑1 (MQ=255)
                                                                  gcagGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAgg                                                   <  1:213754/149‑1 (MQ=255)
                                                                   cagGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGc                                                  <  3:244647/149‑1 (MQ=255)
                                                                         ggtcACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGc                                             <  8:207083/148‑1 (MQ=255)
                                                                              cGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGc                                        >  4:15269/1‑148 (MQ=255)
                                                                               ggCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCa                                      <  4:103390/149‑1 (MQ=255)
                                                                                          cgcAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACAcc                           >  6:242640/1‑149 (MQ=255)
                                                                                          cgcAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACAcc                           <  3:15269/149‑1 (MQ=255)
                                                                                                       tcggtcgCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACACCCCGGCCGACAGGc              >  2:101552/1‑149 (MQ=255)
                                                                                                              cgATCCCAGTCGAACAGGTGGTTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGc                                                                         >  4:30407/1‑83 (MQ=255)
                                                                                                              cgATCCCAGTCGAACAGGTGGTTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGc                                                                         >  4:142403/1‑83 (MQ=255)
                                                                                                              cgATCCCAGTCGAACAGGTGGTTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGc                                                                         <  3:30407/83‑1 (MQ=255)
                                                                                                              cgATCCCAGTCGAACAGGTGGTTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGc                                                                         <  3:142403/83‑1 (MQ=255)
                                                                                                                 tCCCCGTGGAACAGGTGGGTGGTGGCTATCACCTTGGTGAAGCGTTCGGCAATGTCTTCGTGCCAGCGCCTGTGAACGGCTGCCCGGTACCCATCGAGCAGGGCCTGCATTGCCAGCGGATACACCCCGGCCGACAGGCCACCGATCAg    <  6:41289/149‑1 (MQ=255)
                                                                                                                 tCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACACCCCGGCCGACAGGCCACCGATCAg    <  6:221690/149‑1 (MQ=255)
                                                                                                                   ccAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACACCCCGGCCGACAGGCCACCGATCAGGc  <  7:73429/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CCACCAGCTCTTCCGACACTGCCCCGGCCAGCTCCGCCCCGTGATCGCGCGCCAGCACCGCCAGCAGCAGGTCGGTCACGGCCATACCGCCGCAGGCGGTCAGTCGGTCGCGATCCCAGTCGAACAGGTGGCTGGTGGCAATCACCTTGGGGAAGCGTTCGGCAAAGTCATCCTGCCAGCGCCAGTGCACGGCTGCCCGGTACCCATCGAGCAGGCCCAGCATTGCCAGCGGATACACCCCGGCCGACAGGCCACCGATCAGGC  >  NC_002947/5093069‑5093332

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