Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,851,310 T→C 8.2% W18R (TGG→CGG)  PP_5530 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,851,3100TC 8.2% 131.1 / 2.9 49W18R (TGG→CGG) PP_5530hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base C (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

TGAGCTCACCCATATCGTGATCGAGGAGGTCGACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGCACAAGGCCACTGGCGTC  >  NC_002947/2851163‑2851457
                                                                                                                                                   |                                                                                                                                                   
tGAGCTCACCCATATCGTGATCGAGGAGGTCGACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAAATGCCCCGGCGGCGCGAATAAGCGCCTCCAGGGGGGGCTTgg                                                                                                                                                   >  2:552819/1‑150 (MQ=255)
       aCCCATATCGTGATCGAGGAGGTCGACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGaca                                                                                                                                            >  2:413504/1‑150 (MQ=255)
              tCGTGATCGAGGAGGTCGACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCACCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGc                                                                                                                                     >  1:739346/1‑150 (MQ=255)
                               gACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTg                                                                                                                    >  1:552427/1‑150 (MQ=255)
                                acacCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGa                                                                                                                   >  2:517478/1‑150 (MQ=255)
                                 cacCGATAACTGGGGCTACGCCGGGATGAGCCCCGGTGCGAGTGCCTGCCAATACTCAGAGGGATTAGTGCCTAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGaa                                                                                                                  <  2:188219/150‑1 (MQ=255)
                                     gATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGAcc                                                                                                              >  1:585194/1‑150 (MQ=255)
                                     gATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGAcc                                                                                                              >  1:161417/1‑150 (MQ=255)
                                        aaCTGGGGCTACGCCGGCATGCCCACCACCGAGTACAGCAGCAAAGAAACTCAGGGATGAGTTCCGATGTGCCCCGGCGGCGCGAATCAACGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGcc                                                                                                           <  1:271361/150‑1 (MQ=255)
                                         aCTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCa                                                                                                          <  1:523652/150‑1 (MQ=255)
                                            ggggCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCCGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGa                                                                                                       <  2:552427/150‑1 (MQ=255)
                                                cTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCaa                                                                                                   <  1:52207/150‑1 (MQ=255)
                                                 tACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCACCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAAt                                                                                                  <  2:739346/150‑1 (MQ=255)
                                                    gCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACa                                                                                               <  2:288539/150‑1 (MQ=255)
                                                     ccGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACAt                                                                                              >  1:218966/1‑150 (MQ=255)
                                                         cATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAg                                                                                          <  2:161417/150‑1 (MQ=255)
                                                         cATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAg                                                                                          <  2:585194/150‑1 (MQ=255)
                                                               ctccaccGAGTAGCGCACCCAAGACAATACGGGATGAGTTTCGAAGAACGCCGGCGGCGCGAATCAGGGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCGATACATCCAGGCCCGt                                                                                    <  2:437933/148‑1 (MQ=255)
                                                               caccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGt                                                                                    <  2:434868/150‑1 (MQ=255)
                                                               caccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTCGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGcc                                                                                                           <  2:168470/127‑1 (MQ=255)
                                                               caccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTCGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGcc                                                                                                           <  2:168554/127‑1 (MQ=255)
                                                               caccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTCGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGcc                                                                                                           >  1:168554/1‑127 (MQ=255)
                                                               caccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTCGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGcc                                                                                                           >  1:168470/1‑127 (MQ=255)
                                                                 ccaccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGa                                                                                  <  1:413504/150‑1 (MQ=255)
                                                                  caccGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGACTAAACTTGAAGACCGCCATGAGCGATACATCCAGGCCCGTGag                                                                                 >  1:437933/1‑149 (MQ=255)
                                                                         tACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCg                                                                          <  2:666841/150‑1 (MQ=255)
                                                                               agcaAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGAcc                                                                                                              <  1:685340/108‑1 (MQ=255)
                                                                               agcaAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGAcc                                                                                                              >  2:685340/1‑108 (MQ=255)
                                                                                              ggATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAgg                                                                            >  2:140299/1‑127 (MQ=255)
                                                                                              ggATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAgg                                                                            <  1:140299/127‑1 (MQ=255)
                                                                                              ggATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGgtgtgt                                                     >  2:429106/1‑150 (MQ=255)
                                                                                                 tGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGa                                                  >  1:289797/1‑150 (MQ=255)
                                                                                                    gTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGt                                               <  2:89225/150‑1 (MQ=255)
                                                                                                     ttCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTc                                              >  1:14442/1‑150 (MQ=255)
                                                                                                      tCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTcc                                             <  1:429106/150‑1 (MQ=255)
                                                                                                        cGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTAt                                                                      <  2:355637/123‑1 (MQ=255)
                                                                                                        cGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTAt                                                                      >  1:355637/1‑123 (MQ=255)
                                                                                                                cctggcggcGCTGATTCGCGCCGCCGGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTgg                                   >  2:849418/4‑150 (MQ=255)
                                                                                                                    gcggcgCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGaaaa                               <  2:14442/150‑1 (MQ=255)
                                                                                                                           gAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGAGGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTAcc                        >  1:778831/1‑150 (MQ=255)
                                                                                                                              tCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCaa                     <  2:289797/150‑1 (MQ=255)
                                                                                                                              tCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCaa                     >  1:800926/1‑150 (MQ=255)
                                                                                                                                 gcgcCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGca                  >  2:175952/1‑150 (MQ=255)
                                                                                                                                 gcgcCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGca                  >  1:112639/1‑150 (MQ=255)
                                                                                                                                      tCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCAATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGCACAAgg             <  1:552819/150‑1 (MQ=255)
                                                                                                                                                ccTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAgg                                                                                         <  1:781025/64‑1 (MQ=255)
                                                                                                                                                ccTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAgg                                                                                         >  2:781025/1‑64 (MQ=255)
                                                                                                                                                 cTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACTCAAGCACAAGGCCACTGGCgtc  >  1:845738/1‑150 (MQ=255)
                                                                                                                                                 cTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGCACAAGGCCACTGGCgtc  <  1:849418/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGAGCTCACCCATATCGTGATCGAGGAGGTCGACACCGATAACTGGGGCTACGCCGGCATGACCACCACCGAGTACCGCAGCAAAGAAACTCAGGGATGAGTTCCGAACTGCCCCGGCGGCGCGAATCAGCGCCTCCAGGGCGGCCTTGGCTCGACACGCTTGCGAAGGAGAATCCACTTGAAGACCGCCATGAGCAATACATCCAGGCCCGTGACGAAGGCGGTATGGGTGATGCTGGTGTGTTGAAGTCCGACGGATTGGAAAATGTTACCCAAGCACAAGGCCACTGGCGTC  >  NC_002947/2851163‑2851457

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