Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,614,791 T→G 16.7% F388V (TTC→GTC)  PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,614,7910TG16.7% 81.8 / 11.2 36F388V (TTC→GTC) PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (17/13);  new base G (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 = 9.05e-01

CTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACG  >  NC_002947/4614652‑4614938
                                                                                                                                           |                                                                                                                                                   
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                                                                                               >  2:132433/1‑137 (MQ=255)
                                                         cATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCttt                                                                                               <  1:132433/137‑1 (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)
                                                                                                                                           |                                                                                                                                                   
CTATGAAGTGGAGCTGGCCCGCGTCCCACATCCTCAAGGCGGTACGCTGTGGATCTACATCGACAACGACGGTAACCGACTGGAACTGGGCCGCCTCAGTGCCGGTGATGCCTTCATCAGTGCCGTGGACGGCTTGGCCTTCTTCCAGCTGGAGGACGGCCAGACCGTTGTCGAAGATATCAATGAAGGGCTTTATCGGGTCTTCGATACCGACCCGCTGAATCCGAAACGCTCACGGCTGATCAGGCTCGGCGACCGCAATCTCAATTGCCTCGACCTGCTGTACG  >  NC_002947/4614652‑4614938

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