Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,609,896 G→A R224C (CGT→TGT)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,609,8960GA95.5% 62.4 / NA 22R224C (CGT→TGT) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/10);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATC  >  NC_002947/5609754‑5610041
                                                                                                                                              |                                                                                                                                                 
accaAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATAcc                                                                                                                                            >  2:144293/1‑150 (MQ=255)
       ccGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGc                                                                                                                                     <  1:506799/150‑1 (MQ=255)
              ccgccAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCg                                                                                                                                        <  1:104786/140‑1 (MQ=255)
              ccgccAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCg                                                                                                                                        >  2:104786/1‑140 (MQ=255)
                     tCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGAtcgccacaaag                                                                                                                       >  2:391106/1‑144 (MQ=255)
                        ttGTCCAGTTCGACCATGTGCGCGTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCACGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGccgccg                                                                                                                    <  1:573566/149‑1 (MQ=255)
                                      cATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGGAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTttg                                                                                                      >  1:180636/1‑150 (MQ=255)
                                          tgCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTtgttgt                                                                                                  <  1:299344/150‑1 (MQ=255)
                                                 ggTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGa                                                                                           >  2:545270/1‑150 (MQ=255)
                                                               tGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCCTGGCGTAGAAACTGTCGGCACACCTTACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCgg                                                                             >  1:399211/1‑150 (MQ=255)
                                                                     acaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCgg                                                                       >  2:326359/1‑150 (MQ=255)
                                                                     acaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCgg                                                                       <  2:270073/150‑1 (MQ=255)
                                                                             gggTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAAttg                                                               <  2:303396/150‑1 (MQ=255)
                                                                              ggTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAAttgt                                                              >  2:270580/1‑150 (MQ=255)
                                                                                  aCCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGc                                                          <  1:326359/150‑1 (MQ=255)
                                                                                                     gcgcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTcgcc                                       <  1:545270/150‑1 (MQ=255)
                                                                                                       gcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTcgcccg                                     >  2:171297/1‑150 (MQ=255)
                                                                                                              cagTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCgg                              <  1:391106/150‑1 (MQ=255)
                                                                                                                    ggCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGcc                        <  1:171297/150‑1 (MQ=255)
                                                                                                                       cATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGAc                     >  2:300652/1‑150 (MQ=255)
                                                                                                                              tAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGcac              >  1:491000/1‑150 (MQ=255)
                                                                                                                                          gCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTNACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCAtc  >  1:602371/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATC  >  NC_002947/5609754‑5610041

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