Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,794 T→C 17.7% F389L (TTC→CTC)  PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,7940TC17.7% 72.9 / 9.6 36F389L (TTC→CTC) PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base C (3/3);  total (20/16)
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.14e-01

CTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACG  >  NC_002947/4614652‑4614941
                                                                                                                                              |                                                                                                                                                   
cTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGc                                                                                                                                               <  8:9561/149‑1 (MQ=255)
                gCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCGTCCTCCAGCTGGAAGAAGGCCaag                                                                                                                                <  1:19294/148‑3 (MQ=255)
                  ccGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCg                                                                                                                             >  7:37443/1‑149 (MQ=255)
                           acaTCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCg                                                                                                                       >  2:47706/1‑146 (MQ=255)
                                  cAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAg                                                                                                                                                <  6:247474/114‑1 (MQ=255)
                                  cAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAg                                                                                                                                                >  5:247474/1‑114 (MQ=255)
                                   aaGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTg                                                                                                                                             <  8:261045/116‑1 (MQ=255)
                                   aaGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTg                                                                                                                                             >  7:261045/1‑116 (MQ=255)
                                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                           >  6:22944/1‑119 (MQ=255)
                                                  ggATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGtt                                                                                                                           <  5:22944/119‑1 (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                  <  5:24794/140‑1 (MQ=255)
                                                                                                                                    tctGGCCGTCCTCCAGCTGGAAGAAGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTc                  >  6:24794/3‑142 (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)
                                                                                                                                              |                                                                                                                                                   
CTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACGACG  >  NC_002947/4614652‑4614941

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: