Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,934,036 A→T 18.2% L474Q (CTG→CAG)  PP_3467 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,934,0360AT18.2% 48.3 / 6.5 22L474Q (CTG→CAG) PP_3467sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base T (2/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATC  >  NC_002947/3933904‑3934179
                                                                                                                                    |                                                                                                                                               
atCGTCTCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGTTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTcc                                                                                                                                <  2:401372/150‑1 (MQ=255)
                   gTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCa                                                                                                             >  2:100205/1‑150 (MQ=255)
                    ttCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCaa                                                                                                            <  1:400406/150‑1 (MQ=255)
                                                         aCCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAgg                                                                       >  2:121051/1‑150 (MQ=255)
                                                            cGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                    >  2:530184/1‑150 (MQ=255)
                                                            cGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                    <  1:355506/150‑1 (MQ=255)
                                                                       gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGc                                                         <  1:65618/150‑1 (MQ=255)
                                                                          ttCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAg                                                      >  2:552834/1‑150 (MQ=255)
                                                                                        gctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGa                                        <  2:332658/150‑1 (MQ=255)
                                                                                            cTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTATGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtc                                    >  2:509480/1‑150 (MQ=255)
                                                                                                  cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                              >  1:389939/1‑150 (MQ=255)
                                                                                                   aGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAgg                             <  1:530184/150‑1 (MQ=255)
                                                                                                        cTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                             <  1:20813/129‑1 (MQ=255)
                                                                                                        cTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                             >  2:20813/1‑129 (MQ=255)
                                                                                                         ttGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCTCCAGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGc                       <  1:149394/150‑1 (MQ=255)
                                                                                                          tGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCg                      <  1:100205/150‑1 (MQ=255)
                                                                                                                  cgGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGccacca              <  1:552834/150‑1 (MQ=255)
                                                                                                                  atcgacgtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTAGTCGAAGAGAAGCGTGccacca              >  2:158279/7‑150 (MQ=255)
                                                                                                                        gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCAt                                                <  1:345968/110‑1 (MQ=255)
                                                                                                                        gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCAt                                                >  2:345968/1‑110 (MQ=255)
                                                                                                                            cgctgATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTCCa    <  1:158279/146‑1 (MQ=255)
                                                                                                                              cAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATc  >  1:330982/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                               
ATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATC  >  NC_002947/3933904‑3934179

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: