Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 104,876 G→C 12.2% P345P (CCG→CCC)  PP_0099 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947104,8760GC12.2% 114.8 / 8.8 41P345P (CCG→CCC) PP_0099hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base C (3/2);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GAGCTGTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCTTTTTC  >  NC_002947/104727‑105025
                                                                                                                                                     |                                                                                                                                                     
gAGCTGTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAGCGGCGGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCg                                                                                                                                                       <  2:705010/150‑1 (MQ=255)
   ttGTATGAAAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCTAAGGGCGCGAGGCGCAGGAGCTGGATTTCTCGGCGTTTTGCCCACCGGGt                                                                                                                                                    <  1:86209/149‑1 (MQ=255)
   cTGTAAGACAGCTGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGt                                                                                                                                                    <  2:1148782/150‑1 (MQ=255)
     gTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                                  >  2:937088/1‑150 (MQ=255)
     gTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTgg                                                                                                                                                  <  1:208126/150‑1 (MQ=255)
      tATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTggg                                                                                                                                                 <  1:641875/150‑1 (MQ=255)
        tGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCa                                                                                                                                               <  1:335861/150‑1 (MQ=255)
         gACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAg                                                                                                                                              <  1:79290/150‑1 (MQ=255)
             gCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATgg                                                                                                                                          <  2:573403/150‑1 (MQ=255)
              ccGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgc                                                                                                                                         >  1:1356638/1‑150 (MQ=255)
                gatctggCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGgcgc                                                                                                                                       <  2:27707/145‑1 (MQ=255)
                 gTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCg                                                                                                                                                       >  1:427172/1‑133 (MQ=255)
                 gTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCg                                                                                                                                                       <  2:427172/133‑1 (MQ=255)
                     ggCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTc                                                                                                                                  <  2:1356638/150‑1 (MQ=255)
                      gCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGTTGGCCAAGAGCTTGGCAAATGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCtagatc                                                                                                                                 >  1:27707/1‑145 (MQ=255)
                                    ggtgAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccccac                                                                                                                   >  1:1372687/1‑150 (MQ=255)
                                        aaGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGa                                                                                                                               >  1:1305697/1‑134 (MQ=255)
                                        aaGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGa                                                                                                                               <  2:1305697/134‑1 (MQ=255)
                                           gAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGACGGCGCTGTTCTGACATGAACCCCACCCCCcac                                                                                                            >  1:808959/1‑150 (MQ=255)
                                              tATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCgg                                                                                                                                                      <  2:1062352/105‑1 (MQ=255)
                                              tATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCgg                                                                                                                                                      >  1:1062352/1‑105 (MQ=255)
                                                     cGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGTGTTGTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTgc                                                                                                  <  2:697083/150‑1 (MQ=255)
                                                        gTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCa                                                                                               >  2:1313953/1‑150 (MQ=255)
                                                                       cTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCaa                                                                                >  1:868922/1‑150 (MQ=255)
                                                                             cAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCCGGTGGGCaaaa                                                                                                                                            >  1:1215714/1‑81 (MQ=255)
                                                                             cAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCCGGTGGGCaaaa                                                                                                                                            <  2:1215714/84‑4 (MQ=255)
                                                                              aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                         >  2:1232377/1‑150 (MQ=255)
                                                                              aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                         <  1:937088/150‑1 (MQ=255)
                                                                              aaGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACtaataa                                                                         <  1:986193/150‑1 (MQ=255)
                                                                                     tGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCATCTCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCaa                                                                  <  1:397142/150‑1 (MQ=255)
                                                                                     tGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCaa                                                                  <  1:1313953/150‑1 (MQ=255)
                                                                                        aaCGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCaataa                                                               >  2:272915/1‑150 (MQ=255)
                                                                                          cGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAAtt                                                             <  1:480358/150‑1 (MQ=255)
                                                                                            gCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTgc                                                           >  1:493109/1‑150 (MQ=255)
                                                                                                   ttCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATa                                                    <  2:493109/150‑1 (MQ=255)
                                                                                                       aaGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACata                                                >  2:181603/1‑150 (MQ=255)
                                                                                                       aaGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACata                                                >  2:852186/1‑150 (MQ=255)
                                                                                                           gcgcGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTgg                                            <  1:1232377/150‑1 (MQ=255)
                                                                                                                    gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccc                                                                                                                      <  2:1119982/67‑1 (MQ=255)
                                                                                                                    gcAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAaccc                                                                                                                      >  1:1119982/1‑67 (MQ=255)
                                                                                                                     cAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGt                                  >  1:565771/1‑150 (MQ=255)
                                                                                                                                         catctGCCCACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATGCAGtt              >  2:219894/5‑150 (MQ=255)
                                                                                                                                               cccACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACat                                                 <  1:1098120/109‑1 (MQ=255)
                                                                                                                                               cccACCCGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACat                                                 >  2:1098120/1‑109 (MQ=255)
                                                                                                                                                   ccGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCtttt    <  1:272915/150‑1 (MQ=255)
                                                                                                                                                   ccGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCtttt    >  2:229316/1‑150 (MQ=255)
                                                                                                                                                     gggTGGGCAGACGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCTTTTTc  <  2:808959/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCTGTATGACAGCCGGTTCGGCGTGCGCATGCGTGGTGAAGGAGTATTTGCCGAGTTGCTGGCGCAGCGCTTTGCCAAGGCAATGAAACGGCTGGATTTCGAAGGGCGCGAGGCGCAGGAGCTGGATTGCTCGGCGTTTTGCCCACCGGGTGGGCAGATGGCGCTGTTCTGACATGAACCCCACCCCCCACAGGGGGCTGCGCATTGCTCACAACGCAACTAATAACAACCAATAATTGCGCTAATACATATTGGCAACTGATGTTTTTTTGCTATTATCCAGTTCCCGCCTTTTTC  >  NC_002947/104727‑105025

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