Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,920,763 G→A 9.3% G399S (GGC→AGC)  PP_4328 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,920,7630GA 9.3% 128.0 / 3.6 43G399S (GGC→AGC) PP_4328hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/22);  new base A (2/2);  total (19/24)
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.51e-01

AACAGCTGCTGCGCCTGGCTGCCGCTGCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCG  >  NC_002947/4920629‑4920897
                                                                                                                                      |                                                                                                                                      
aaCAGCTGCTGCGCCTGGCTGCCGCTTCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACgg                                                                                                                                       <  2:23968/136‑1 (MQ=255)
aaCAGCTGCTGCGCCTGGCTGCCGCTGCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACgg                                                                                                                                       <  2:366001/136‑1 (MQ=255)
 aCAGCTGCTGCGCCTGGCTGCCGCTGCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGc                                                                                                                                      <  1:554095/136‑1 (MQ=255)
                ggCTGCCGCTGCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCg                                                                                                                       >  1:496471/1‑136 (MQ=255)
                             gTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCACACCACCTGGCAAACTGAAGAGACGATCCGCCACGGCCAGGACTTCAACCCGCTGCAAGTaaaaa                                                                                                          >  2:123859/1‑135 (MQ=255)
                               ttCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTg                                                                                                        >  2:317544/1‑136 (MQ=255)
                                 cACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCa                                                                                                      <  1:11936/136‑1 (MQ=255)
                                       tGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCga                                                                                                <  1:442888/136‑1 (MQ=255)
                                          agagCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCgaaga                                                                                             <  1:76020/136‑1 (MQ=255)
                                          agagCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCgaaga                                                                                             >  1:319611/1‑136 (MQ=255)
                                           gagCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCgaagaa                                                                                            >  1:531305/1‑136 (MQ=255)
                                           gagCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCgaagaa                                                                                            >  1:415064/1‑136 (MQ=255)
                                               cACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAAccc                                                                                        <  1:197843/136‑1 (MQ=255)
                                                   ccATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACAGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCcaca                                                                                    >  1:559629/1‑136 (MQ=255)
                                                   ccATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACAGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCcaca                                                                                    >  1:562124/1‑136 (MQ=255)
                                                                gCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGAc                                                                       <  2:220781/136‑1 (MQ=255)
                                                                     aaaCCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGccc                                                                  >  1:73086/1‑135 (MQ=255)
                                                                     aaaCCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCa                                                                  >  2:100471/1‑136 (MQ=255)
                                                                       aCCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAAc                                                                <  1:197574/136‑1 (MQ=255)
                                                                               gCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGcc                                                        <  2:478380/136‑1 (MQ=255)
                                                                                   gAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAg                                                    <  1:317544/136‑1 (MQ=255)
                                                                                         ggCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTc                                              <  2:366136/136‑1 (MQ=255)
                                                                                         ggCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTc                                              <  1:249692/136‑1 (MQ=255)
                                                                                            aaCCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGc                                                         <  1:124052/122‑1 (MQ=255)
                                                                                            aaCCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGc                                                         >  2:124052/1‑122 (MQ=255)
                                                                                              ccTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACAGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACaa                                         <  2:562124/136‑1 (MQ=255)
                                                                                              ccTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACAGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACaa                                         <  2:559629/136‑1 (MQ=255)
                                                                                                        accTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGgc                               >  2:398992/1‑136 (MQ=255)
                                                                                                        accTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGgc                               <  1:110398/136‑1 (MQ=255)
                                                                                                        accTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGgc                               <  1:53054/136‑1 (MQ=255)
                                                                                                            ggCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCAt                           >  1:395874/1‑136 (MQ=255)
                                                                                                            ggCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGGCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCAt                           >  2:96313/1‑136 (MQ=255)
                                                                                                              cAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCg                         <  1:465803/136‑1 (MQ=255)
                                                                                                               aaaCTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGa                        <  1:282595/136‑1 (MQ=255)
                                                                                                               aaaCTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGa                        <  1:48615/136‑1 (MQ=255)
                                                                                                                       gTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCt                >  1:181063/1‑136 (MQ=255)
                                                                                                                       gTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCt                >  1:528553/1‑136 (MQ=255)
                                                                                                                         gCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTc              <  1:163303/136‑1 (MQ=255)
                                                                                                                         gCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTc              <  2:367277/136‑1 (MQ=255)
                                                                                                                         cccGATCCGCCACGGCCAGGACTTCCTCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTc              <  2:206348/135‑1 (MQ=255)
                                                                                                                             aTCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGAcc          >  2:497699/1‑136 (MQ=255)
                                                                                                                               ccGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTg        <  2:415064/136‑1 (MQ=255)
                                                                                                                               ccGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTg        <  2:531305/136‑1 (MQ=255)
                                                                                                                                     cGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCg  >  2:536818/1‑136 (MQ=255)
                                                                                                                                     cGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCATGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCg  >  1:443344/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AACAGCTGCTGCGCCTGGCTGCCGCTGCCGTTTCACGCTTGCAGAGCCACGCCATGAGCAGCCTGCAACAAACCGGCACGCTGGAAAACGGCAACCTGCAGACCACCTGGCAAACTGAAGTGCCGATCCGCCACGGCCAGGACTTCATCCCGCTGCAAGTCAAACTGCAGCGCGAAGAAACCCCACAACAACAGGCCGACCGCCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCG  >  NC_002947/4920629‑4920897

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑