Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,808 C→A 14.6% D393E (GAC→GAA)  PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,8080CA14.6% 98.5 / 9.8 41D393E (GAC→GAA) PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base A (3/3);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

GCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAGGGCGCCT  >  NC_002947/4614668‑4614954
                                                                                                                                            |                                                                                                                                                  
gCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCGTCCTCCAGCTGGAAGAAGGCCaag                                                                                                                                             <  1:19294/148‑3 (MQ=255)
  ccGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCg                                                                                                                                          >  7:37443/1‑149 (MQ=255)
           acaTCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCg                                                                                                                                    >  2:47706/1‑146 (MQ=255)
                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                                        <  5:22944/119‑1 (MQ=255)
                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                                        >  6:22944/1‑119 (MQ=255)
                                        aCATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCTTGGACGACTTGGCCTTATTCCAGCTGGAGGACGGCCAGACCGTTTTAGAAg                                                                                                                                 >  3:26735/1‑120 (MQ=255)
                                        aCATCGAAAACGACGGTAACGGACTGGAACTGGGCCGCCTAAGTGCCGGTGATGCCTTCATCAGTGCAGTGGACTGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAg                                                                                                                                 <  4:26735/120‑1 (MQ=255)
                                         cATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCttt                                                                                                               <  1:132433/137‑1 (MQ=255)
                                         cATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCttt                                                                                                               >  2:132433/1‑137 (MQ=255)
                                               cAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGat                                                                                                                               >  8:78519/1‑115 (MQ=255)
                                               cAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGat                                                                                                                               <  7:78519/115‑1 (MQ=255)
                                                                 ggAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAAc                                                                           >  2:210607/1‑149 (MQ=255)
                                                                  gAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGAGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTATTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAAACCGAAACg                                                                          >  4:238370/1‑149 (MQ=255)
                                                                      tGGGCCGCCTCAGTGCCGGTGATGCCTTCATGAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCa                                                                      >  1:91743/1‑149 (MQ=255)
                                                                            gccTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCt                                                                 >  5:293363/1‑148 (MQ=255)
                                                                                     ccGGTGATGCCTTCATCAGTGCCGTGGACGGCTGGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTCTCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGcgc                                                       <  6:159849/149‑3 (MQ=255)
                                                                                      cGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCg                                                      >  2:191653/1‑149 (MQ=255)
                                                                                         tGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGc                                                    >  2:270046/1‑148 (MQ=255)
                                                                                           aTGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGAc                                                 <  1:191653/149‑1 (MQ=255)
                                                                                               cTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGc                                              >  4:273775/1‑148 (MQ=255)
                                                                                                  catcaGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAAt                                           <  1:177293/148‑1 (MQ=255)
                                                                                                                aCGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGAc                            <  3:273775/149‑1 (MQ=255)
                                                                                                                  ggCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACct                          >  8:92121/1‑149 (MQ=255)
                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                               >  6:24794/3‑142 (MQ=255)
                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                               <  5:24794/140‑1 (MQ=255)
                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc                        >  1:224613/3‑149 (MQ=255)
                                                                                                                    cTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgc                        >  6:2071/1‑149 (MQ=255)
                                                                                                                     ttGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACctgct                       <  1:47706/149‑1 (MQ=255)
                                                                                                                        gCCTTCATCCAGCTGGACGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGt                     <  4:189957/148‑1 (MQ=255)
                                                                                                                        gCCGTCCTCGAGCTGTAAGAAGGCCAGACCGTTGTCGTAGATATCTATGTAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAATCGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGcgt                                <  4:46862/137‑3 (MQ=255)
                                                                                                                        gCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTAGGCGACCGCAATCTAAATTGCCt                                >  3:46862/1‑137 (MQ=255)
                                                                                                                          cttcttCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacg                  >  3:249166/1‑149 (MQ=255)
                                                                                                                              ttcCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTacgacg               >  3:182364/1‑148 (MQ=255)
                                                                                                                                  aGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACaa          <  7:92121/149‑1 (MQ=255)
                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCg                                               >  5:45471/1‑111 (MQ=255)
                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCg                                               <  6:45471/111‑1 (MQ=255)
                                                                                                                                   gCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAg         <  8:37443/149‑1 (MQ=255)
                                                                                                                                       gaggaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAAt                                     <  3:81426/117‑1 (MQ=255)
                                                                                                                                       gaggaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGAAACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAAt                                     >  4:81426/1‑117 (MQ=255)
                                                                                                                                          gaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAGGGCGCCt  <  5:2071/149‑1 (MQ=255)
                                                                                                                                          gaCGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAGGGCGCCt  <  1:210607/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACGAACAAGGGCGCCT  >  NC_002947/4614668‑4614954

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: ‑