Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,132,360 G→A 8.5% R40H (CGC→CAC)  PP_4518 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,132,3600GA 8.5% 129.9 / 2.4 47R40H (CGC→CAC) PP_4518hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

ACGCCACCACTGAGCAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGCCTGTCGT  >  NC_002947/5132222‑5132503
                                                                                                                                          |                                                                                                                                               
acgccacCACTGAGCAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTcgccg                                                                                                                                       >  8:112053/1‑149 (MQ=255)
acgccacCACTGAGCAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTcgccg                                                                                                                                       >  6:298214/1‑149 (MQ=255)
   ccaccaCTGAGCAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCgcgc                                                                                                                                    <  6:69856/149‑1 (MQ=255)
              cAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActg                                                                                                                         >  6:334035/1‑149 (MQ=255)
               aaGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATGActgg                                                                                                                        >  6:309440/1‑149 (MQ=255)
                aGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActggc                                                                                                                       <  7:112053/149‑1 (MQ=255)
                    gtgCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGt                                                                                                                   >  1:22445/1‑149 (MQ=255)
                     tgCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTc                                                                                                                  <  2:22445/149‑1 (MQ=255)
                               gctgctgGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCatcatc                                                                                                           >  4:142083/1‑146 (MQ=255)
                               gctgctgGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCatcatc                                                                                                           <  3:142083/146‑1 (MQ=255)
                                ctgctgGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATccc                                                                                                         <  5:298214/147‑1 (MQ=255)
                                      gCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTg                                                                                                   <  5:334035/147‑1 (MQ=255)
                                      gCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCACAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTggg                                                                                                 >  3:250331/1‑149 (MQ=255)
                                      gCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCACAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTggg                                                                                                 >  3:245561/1‑149 (MQ=255)
                                       ccACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCACAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTgggg                                                                                                <  4:250331/149‑1 (MQ=255)
                                       ccACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCACAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTgggg                                                                                                <  4:245561/149‑1 (MQ=255)
                                           cagcaccctCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCTCTGGCTGGTCCATCATCCCAAGCTGGGGCCAt                                                                                            <  1:38323/142‑1 (MQ=255)
                                             gctggctCAGCGTCGCGCTGGGGGTGCTGGGGATCTTCCTGCCGGTATTGCCCACCACCCCATTCCTGGTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTACATCATCCCAAGCTGGGGTCATg                                                                                           >  3:37687/1‑148 (MQ=255)
                                             gctggctCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTc                                                                                                                                 <  8:88503/110‑1 (MQ=255)
                                             gctggctCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTc                                                                                                                                 >  7:88503/1‑110 (MQ=255)
                                                   tCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTACCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGt                                                                                    >  1:251142/1‑149 (MQ=255)
                                                                 ggggTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAgtgg                                                                      >  2:162991/1‑149 (MQ=255)
                                                                 ggggTGATGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGTCATGGATTCGTGACTACCTCAgtg                                                                       <  4:37687/148‑1 (MQ=255)
                                                                           ggATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATc                                                            <  6:242317/149‑1 (MQ=255)
                                                                               atTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATcccc                                                         >  2:57305/1‑148 (MQ=255)
                                                                                 tCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCt                                                       >  3:323511/1‑148 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActggc                                                                                                                       <  4:116376/81‑1 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActggc                                                                                                                       >  3:116376/1‑81 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActgg                                                                                                                        >  5:191821/1‑80 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCActgg                                                                                                                        <  6:191821/80‑1 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATcc                                                                                                          >  8:130665/1‑94 (MQ=255)
                                                                                    tGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATcc                                                                                                          <  7:130665/94‑1 (MQ=255)
                                                                                      ccGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAggg                                                 <  4:323511/149‑1 (MQ=255)
                                                                                               cccACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCACAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTa                                        <  1:57305/149‑1 (MQ=255)
                                                                                                 caccacCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACg                                      >  5:254030/1‑149 (MQ=255)
                                                                                                  accaccCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGc                                     >  1:202128/1‑149 (MQ=255)
                                                                                                       cccATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGCGGCCATGGATTCGTGACTACCTCAGAGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCg                                >  3:140449/1‑149 (MQ=255)
                                                                                                          aTTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGcc                             <  2:251142/149‑1 (MQ=255)
                                                                                                               ttcTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCTCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTAc                                       <  2:229354/132‑1 (MQ=255)
                                                                                                               tGCTTTTGGCGGCGGCCTGCTTTGACCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATACCCCTCAAGGGCAAGGTCTAc                                       >  1:229354/1‑134 (MQ=255)
                                                                                                                         ggcggcCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAgg                                                                >  3:328346/1‑99 (MQ=255)
                                                                                                                         ggcggcCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAgg                                                                <  4:328346/99‑1 (MQ=255)
                                                                                                                             gcCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTgg          <  2:202128/149‑1 (MQ=255)
                                                                                                                              ccTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGc         <  6:254030/149‑1 (MQ=255)
                                                                                                                               cTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGcc        >  6:222092/1‑149 (MQ=255)
                                                                                                                                     tGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAACCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGCCTGTCGt  >  4:188703/1‑149 (MQ=255)
                                                                                                                                      gCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGCCTGTCGt  <  5:222092/148‑1 (MQ=255)
                                                                                                                                          gCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGt                                                                                    <  6:245506/62‑1 (MQ=255)
                                                                                                                                          gCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGt                                                                                    >  5:245506/1‑62 (MQ=255)
                                                                                                                                          |                                                                                                                                               
ACGCCACCACTGAGCAAGCTGTGCGCTACCTGCTGCTGGCCACTGGCTGGCTCAGCGTCGCGCTGGGGGTGCTGGGGATATTCCTGCCGGTATTGCCCACCACCCCATTCCTGCTTTTGGCGGCGGCCTGCTTTGCCCGCAGCTCGCCGCGCTTCCATCACTGGCTGGTCCATCATCCCAAGCTGGGGCCATGGATTCGTGACTACCTCAGTGGTGAAGGCATCCCCCTCAAGGGCAAGGTCTACGCCATCGGCCTGATGTGGGCCAGCATTGGCCTGTCGT  >  NC_002947/5132222‑5132503

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