Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,116,213 T→C 9.3% L303P (CTG→CCG)  PP_2732 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,116,2130TC 9.3% 114.1 / 2.3 43L303P (CTG→CCG) PP_2732hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (2/2);  total (22/21)
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.58e-01

GAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCATCGGG  >  NC_002947/3116071‑3116361
                                                                                                                                              |                                                                                                                                                    
gAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGAGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGAGGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTTGCCGCCTGAtgagc                                                                                                                                               >  2:443565/1‑146 (MQ=255)
    cGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGc                                                                                                                                           <  2:54486/150‑1 (MQ=255)
     gTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCg                                                                                                                                          >  1:584432/1‑150 (MQ=255)
      tACCTGTTTCATGCGCGCATCAGGGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGt                                                                                                                                         <  1:668382/150‑1 (MQ=255)
         cTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAg                                                                                                                                      <  1:217577/150‑1 (MQ=255)
           gTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCCGATCAGCTGGCgaatgaa                                                                                                                                    <  1:582815/150‑7 (MQ=255)
                      gcATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTa                                                                                                                                    <  1:650895/139‑1 (MQ=255)
                      gcATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTa                                                                                                                                    >  2:650895/1‑139 (MQ=255)
                        aTCAGTGCGGCGCACAACATCGGCCTTCTGGATTTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATc                                                                                                                       <  1:25134/150‑1 (MQ=255)
                            gTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTAtt                                                                                                                   <  1:200718/150‑1 (MQ=255)
                                     tCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCTCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCa                                                                                                          >  2:459504/1‑150 (MQ=255)
                                                   cTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGc                                                                                                                                           >  2:651355/1‑103 (MQ=255)
                                                   cTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGc                                                                                                                                           <  1:651355/103‑1 (MQ=255)
                                                     ggATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCgg                                                                                          <  1:630069/150‑1 (MQ=255)
                                                     ggATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCgg                                                                                          >  2:1183919/1‑150 (MQ=255)
                                                      gATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCggg                                                                                         >  1:585112/1‑150 (MQ=255)
                                                               cgtcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCAc                                                                                <  1:1183919/150‑1 (MQ=255)
                                                                gtcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAg                                                                                                              >  2:1125741/1‑119 (MQ=255)
                                                                gtcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAg                                                                                                              <  1:1125741/119‑1 (MQ=255)
                                                                 tcTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCt                                                                              <  2:585112/150‑1 (MQ=255)
                                                                    gtgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCTCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGaa                                                                           <  1:459504/150‑1 (MQ=255)
                                                                    gtgCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGaa                                                                           <  1:706836/150‑1 (MQ=255)
                                                                       cgacgaCGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGGCTGCCGGCATATGCCGGGATCAATCACCTGAACaa                                                                        >  2:31812/1‑150 (MQ=255)
                                                                        gacgacGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAAt                                                                       >  2:932320/1‑150 (MQ=255)
                                                                         acgacgTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCt                                                                                                                      <  2:1422109/102‑1 (MQ=255)
                                                                         acgacgTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCt                                                                                                                      >  1:1422109/1‑102 (MQ=255)
                                                                                 gACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATACCGGGCTCAATCACCTGAACAATCACCGGgcg                                                              >  1:633095/1‑150 (MQ=255)
                                                                                    aaaGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTg                                                           >  1:33743/1‑150 (MQ=255)
                                                                                                  ggggCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTg                                             >  1:1100954/1‑150 (MQ=255)
                                                                                                        ggTGCCGCTCGATGCAGCCGAAGGCTGCATTCGCCAGCTGATCGGCTGGCGTGCGTACGTGCGTGGCCGCTGTTGGCGGCGCATGCCGGCATGTGCCGGGCTCTATCACCTGAACAATCACCGGTAGCAGCAAAAATTCTAATGGACCgg                                       >  2:1369857/1‑150 (MQ=255)
                                                                                                          tGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGgc                                     >  2:884513/1‑150 (MQ=255)
                                                                                                            ccGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGgcgc                                   <  2:584432/150‑1 (MQ=255)
                                                                                                                      gaGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGTCGGCATATGCCGGGCTCAATCACCTGACCAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATgc                         <  1:893667/149‑1 (MQ=255)
                                                                                                                        gCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATgcgc                       >  1:831316/1‑150 (MQ=255)
                                                                                                                        gCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATgcgc                       <  1:932320/150‑1 (MQ=255)
                                                                                                                         ccGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCAGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGc                             >  1:1322466/1‑143 (MQ=255)
                                                                                                                         ccGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCAGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGc                             <  2:1322466/143‑1 (MQ=255)
                                                                                                                           gAAGGCTTCATTCGACAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGc                    >  2:235348/1‑150 (MQ=255)
                                                                                                                                  actcacGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGgcgc             >  1:801311/6‑150 (MQ=255)
                                                                                                                                     ttCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAg          >  2:490905/1‑150 (MQ=255)
                                                                                                                                      acGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCgg                                                                                                    >  1:768378/2‑59 (MQ=37)
                                                                                                                                      acGCCAGCCGATCAGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCgg                                                                                                    <  2:768378/58‑1 (MQ=37)
                                                                                                                                             cTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATACCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCATCggg  <  2:633095/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GAGCCGTACCTGTTTCATGCGCGCATCAGTGCGGCGCTCAACATCGGCCTTCTGGATGTGCGTCGTCTGTGCGACGACGTCGACAAAGCCTGGCGCGAGGGGCGGGTGCCGCTCAATGCAGCCGAAGGCTTCATTCGCCAGCTGATCGGCTGGCGTGAGTACGTGCGTGGCATCTATTGGCAGCGCATGCCGGCATATGCCGGGCTCAATCACCTGAACAATCACCGGGCGCTGCCAGAATTCTACTGGACCGGGCGCACCCGCATGCGCTGCATGGCGCAAGCCATCGGG  >  NC_002947/3116071‑3116361

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