Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,315,892 G→A 9.5% V297I (GTA→ATA)  PP_3787 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,315,8920GA 9.5% 111.9 / 2.9 42V297I (GTA→ATA) PP_3787hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base A (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGCGCTGCGACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAC  >  NC_002947/4315751‑4316040
                                                                                                                                             |                                                                                                                                                    
gcgGGCGCTGCGACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCca                                                                                                                                               <  5:205971/149‑1 (MQ=255)
    gcgcTGCGACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCcacaca                                                                                                                                           <  1:209562/149‑1 (MQ=255)
        tGCGACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGc                                                                                                                                         >  5:301696/1‑147 (MQ=255)
           gACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGa                                                                                                                                    <  3:15660/149‑1 (MQ=255)
               cAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAgcgc                                                                                                                                <  8:133208/149‑1 (MQ=255)
                     gCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCATACGCCACACAGCAAGGAGCGCAGCAt                                                                                                                           >  5:317455/1‑148 (MQ=255)
                     gCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCATACGCCACACAGCAAGGAGCGCAGCAt                                                                                                                           >  3:159761/1‑148 (MQ=255)
                         aCTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAgg                                                                                                                                     <  1:112384/134‑1 (MQ=255)
                         aCTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAgg                                                                                                                                     >  2:112384/1‑134 (MQ=255)
                         aCTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAAc                                                                                                                       <  1:180845/148‑1 (MQ=255)
                         aCTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAAc                                                                                                                       <  1:229170/148‑1 (MQ=255)
                         aCTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTAAGCGGATGTTCATCAGGGAATCAGTCTCGTACGCCACCAAGCAAGGAGCGCAGCATGAAc                                                                                                                       >  5:76197/1‑148 (MQ=255)
                          cTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCa                                                                                                                     <  5:144352/149‑1 (MQ=255)
                             cTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAgg                                                                                                                   >  5:77853/1‑148 (MQ=255)
                                    aGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACt                                                                                                            >  2:18493/1‑148 (MQ=255)
                                            ccTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCAc                                                                                                   <  8:25645/149‑1 (MQ=255)
                                             cTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGGATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACg                                                                                                  >  7:176955/1‑149 (MQ=255)
                                              tGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACgg                                                                                                 >  5:197075/1‑149 (MQ=255)
                                                aTGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACgg                                                                                                 <  5:269337/147‑1 (MQ=255)
                                                aTGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACgg                                                                                                 <  7:298913/147‑1 (MQ=255)
                                                aTGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGa                                                                                               <  1:17509/149‑1 (MQ=255)
                                                                          gCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAACGAGCGCAGCTTGAACCAGGCAGAACTCACGATCACGGGAGGGCCTGCCTCCCCCCTGCACTtcac                                                                     >  5:257681/1‑147 (MQ=255)
                                                                              aTAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTg                                                                 <  6:301696/149‑1 (MQ=255)
                                                                                 gctgCCGTCGAGCACTACCTTCTGATCTGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCAGACAGCAAGGCGCGCAGCATGAACCAGGCATTACTGACGATCACGGGTGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTc                                                              <  6:76197/146‑1 (MQ=255)
                                                                                 gCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGt                                                               <  8:176955/148‑1 (MQ=255)
                                                                                 gCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTc                                                              >  4:103754/1‑149 (MQ=255)
                                                                                          gaGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATc                                                     >  4:227957/1‑149 (MQ=255)
                                                                                              aCTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCATACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGc                                                 <  4:159761/149‑1 (MQ=255)
                                                                                              aCTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCATACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGc                                                 <  6:317455/149‑1 (MQ=255)
                                                                                                tACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCg                                                >  1:22795/1‑148 (MQ=255)
                                                                                                      cTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTatca                                           <  6:197075/147‑1 (MQ=255)
                                                                                                        gATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCa                                       >  1:333848/1‑149 (MQ=255)
                                                                                                             gCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGCTCACGGGAGGGCCTGCCTACCCCGTGCAATTCTAGCTGGTCGAATACATCGAGCGAATCATACAACGc                                   >  7:156707/1‑148 (MQ=255)
                                                                                                                 tGCGGATCTTCATCAGGGAATCAGTCTCGTACTCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTTGTCGAATACATCGAGCGTATCATCCACCGCTAc                                <  3:227957/147‑1 (MQ=255)
                                                                                                                    ggATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATc                           >  5:112603/1‑149 (MQ=255)
                                                                                                                           catcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAAt                       <  2:333848/146‑1 (MQ=255)
                                                                                                                             tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGc                  >  5:188632/1‑149 (MQ=255)
                                                                                                                               aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCg                 <  6:188632/148‑1 (MQ=255)
                                                                                                                                 ggAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCa              >  5:73764/1‑149 (MQ=255)
                                                                                                                                     tCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTc           >  6:241050/1‑148 (MQ=255)
                                                                                                                                       aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac        <  6:77853/149‑1 (MQ=255)
                                                                                                                                         tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac        >  7:268518/1‑147 (MQ=255)
                                                                                                                                             gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc  <  3:103754/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCGGGCGCTGCGACGCAAATGGCTCACTGCTTCGCCAGCGATTACCTGATGGCGCTCGCCGAGCTTTCGCAGCGGCAAATAGCGGCCGACGAGCACTACCTTCTGATCGGCTATGCGGATCTTCATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAC  >  NC_002947/4315751‑4316040

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