Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,593,944 C→T 13.7% R10C (CGC→TGC)  PP_5510 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,593,9440CT13.7% 66.6 / 3.7 30R10C (CGC→TGC) PP_5510hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (2/2);  total (14/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 = 7.80e-01

AGCGACCACCGTGGCTCCGCGTCCGGCCTATGTGGCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCGCTACGAACTGAGCAACTACCTCAACGCTT  >  NC_002947/2593796‑2594088
                                                                                                                                                     |                                                                                                                                                
accGACCCCCGTGCCTCCGCGGCCGGCCTATGTGGCCCCAAGCCCAGCCTACGTTCCACCGCGAGCCACCTACCGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCcgc                                                                                                                                                  <  8:248202/147‑1 (MQ=255)
aGCGACCACCGTGGCGCCGCGTCCTGCCTATGTGGCTCCACGCCCAGCCTACGTTTCACCGCTAGCCAGCTACAGCTATCCCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAACCCATGACGACTCCAGCCCAACGGCCTCcgc                                                                                                                                                  <  8:248204/149‑1 (MQ=255)
          gTGGCTCCGCGTCCGGCCTATGTGGCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAg                                                                                                                                         >  4:255136/1‑148 (MQ=255)
                         gCCTATGTGGCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCAt                                                                                                                         >  4:190055/1‑149 (MQ=255)
                                  gCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTgc                                                                                                                             >  2:49851/1‑136 (MQ=255)
                                  gCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTgc                                                                                                                             <  1:49851/136‑1 (MQ=255)
                                      cACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCTAACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCATCCCAACGGCATCAGAGCCTGAAAGAACTGTACTTGCGAATTGAGCGTGTGGa                                                                                                             >  3:142795/1‑148 (MQ=255)
                                            cAGCCTACGTTGCACCGCGAGCCAGCTACAGTCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACgg                                                                                                       >  1:28245/1‑149 (MQ=255)
                                            cAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGa                                                                                                      <  3:255136/149‑1 (MQ=255)
                                            cAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGa                                                                                                      <  1:226131/149‑1 (MQ=255)
                                                        cACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCgg                                                                                          >  8:182913/1‑149 (MQ=255)
                                                         aCCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGTGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGt                                                                                         >  1:252385/1‑149 (MQ=255)
                                                         aCCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGTGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGt                                                                                         >  1:248348/1‑149 (MQ=255)
                                                                 ccAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGac                                                                                 <  2:61837/149‑1 (MQ=255)
                                                                 ccAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCACCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTTAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCgaga                                                                                  <  2:45396/148‑1 (MQ=255)
                                                                            ctaTCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTgg                                                                       <  7:182913/149‑1 (MQ=255)
                                                                                           cTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAg                                                        <  1:17825/149‑1 (MQ=255)
                                                                                              ctcCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGt                                                      >  2:188986/1‑148 (MQ=255)
                                                                                              ctcCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGt                                                      >  2:208641/1‑148 (MQ=255)
                                                                                                 cAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACt                                                  >  4:112526/1‑149 (MQ=255)
                                                                                                 cAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACt                                                  >  2:75112/1‑149 (MQ=255)
                                                                                                     ttCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCt                                                                         >  3:83143/1‑122 (MQ=255)
                                                                                                     ttCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCt                                                                         <  4:83143/122‑1 (MQ=255)
                                                                                                             cagTATGTAAGCCATTACGACTCCATCCCAACTTCCTCCGCGCCTGAAAGAACTTTCCTTGCGCATTTAGCGTGTGTACAACGGAGTCATCGTCCTGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGTTGACTGCAACGTGTTCa                                      <  4:142795/147‑1 (MQ=255)
                                                                                                              cGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCAc                                     <  1:208641/149‑1 (MQ=255)
                                                                                                              cGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCAc                                     <  1:188986/149‑1 (MQ=255)
                                                                                                                     aaGCCATGACGACTCCAGCCCAACGGCCTCCGTGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCg                               <  2:252385/148‑1 (MQ=255)
                                                                                                                     aaGCCATGACGACTCCAGCCCAACGGCCTCCGTGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCg                               <  2:248348/148‑1 (MQ=255)
                                                                                                                        ccATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCCTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTTACTGCAACGTGTTCACCAATCGCTAc                           <  3:185562/149‑1 (MQ=255)
                                                                                                                                   ccAGCCCAACGGCCTCCGCGCCTGAAAGAACTGACCTAGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGCTGACTGAAACGTGTTCACCAATCCCTACGAACTGAGCa                 >  7:65888/1‑148 (MQ=255)
                                                                                                                                         cAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCGCTACGAACTGAGCAACTACCt          <  1:75112/149‑1 (MQ=255)
                                                                                                                                          aaCGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGac                                                                                 <  8:211979/77‑1 (MQ=255)
                                                                                                                                                   cgcgCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAAGATCCTGGGACTGGCAAACGAAGGTCACTGCTTCGTGTTCTCCATTCGCTACGAACTTAGCAACTACCTCAACGCtt  <  2:28245/147‑1 (MQ=255)
                                                                                                                                                   cgcgCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCGCTACGAACTGAGCAACTACCTCAACGCtt  <  3:190055/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCGACCACCGTGGCTCCGCGTCCGGCCTATGTGGCTCCACGCCCAGCCTACGTTGCACCGCGAGCCAGCTACAGCTATCGCACCAGCAGCTTCTCCAGCTTCCGCTCCCGTAGGTAAGCCATGACGACTCCAGCCCAACGGCCTCCGCGCCTGAAAGAACTGTCCTTGCGCATTGAGCGTGTGGACAACGGAGTCATCGTCCGGTCCCGAGACAACATCCTGGGCCTGGCAAACGAAGGTGACTGCAACGTGTTCACCAATCGCTACGAACTGAGCAACTACCTCAACGCTT  >  NC_002947/2593796‑2594088

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: