Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,548 G→C 8.5% G11G (GGC→GGG)  PP_4850 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,518,5480GC 8.5% 140.1 / 4.7 47G11G (GGC→GGG) PP_4850membrane protein
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5518400‑5518692
                                                                                                                                                    |                                                                                                                                                
gTTCGAACACCAGCAGACCGGCCAGCACGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATCACGATGATCAGGCTGGCCAGGGGGATCTGACGGGGGGGAAAGGTCATACAGCCTAGCAGGCAGACAGCGGGAATg                                                                                                                                                  >  8:127566/1‑149 (MQ=255)
     aaCACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATATGACGGGGGGTAAAGGTCAGACAGCCTCGCCGGCAGCCAGCGGGAATGCCg                                                                                                                                               >  8:263269/1‑147 (MQ=255)
          cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTg                                                                                                                                        >  6:34207/1‑149 (MQ=255)
           agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                       >  4:261626/1‑149 (MQ=255)
             cagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                                     <  1:176949/149‑1 (MQ=255)
                  cGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                       <  8:218365/142‑1 (MQ=255)
                  cGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                       >  7:218365/1‑142 (MQ=255)
                    gCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                              <  8:23119/149‑1 (MQ=255)
                     ccAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggc                                                                                                                             <  5:221963/149‑1 (MQ=255)
                           ggCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCa                                                                                                                       <  5:34207/149‑1 (MQ=255)
                                        cTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAg                                                                                                          <  3:261626/149‑1 (MQ=255)
                                           cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGt                                                                                                        >  1:322189/1‑148 (MQ=255)
                                                gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGcc                                                                                                                                      <  2:34562/113‑1 (MQ=39)
                                                gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGcc                                                                                                                                      >  1:34562/1‑113 (MQ=39)
                                                        agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAg                                                                                          <  3:86920/149‑1 (MQ=255)
                                                           agcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                               >  8:206098/1‑125 (MQ=255)
                                                           agcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                               <  7:206098/125‑1 (MQ=255)
                                                            gcagcCCGGCCATGACGAAGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGACATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGc                                                                                        >  8:221309/1‑147 (MQ=255)
                                                             cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAAt                                                                                     <  2:322189/149‑1 (MQ=255)
                                                               gcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCg                                                                                   >  3:32556/1‑149 (MQ=255)
                                                               gcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCg                                                                                   >  3:37084/1‑149 (MQ=255)
                                                                cccGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGc                                                                                  >  7:277498/1‑149 (MQ=255)
                                                                  cGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTc                                                                                <  7:221309/149‑1 (MQ=255)
                                                                              tgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGACAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTg                                                                    >  1:137697/1‑149 (MQ=255)
                                                                                   aGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                               >  3:252515/1‑117 (MQ=255)
                                                                                   aGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                               <  4:252515/117‑1 (MQ=255)
                                                                                       tGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAAt                                                                                     <  5:358949/123‑1 (MQ=255)
                                                                                       tGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAAt                                                                                     >  6:358949/1‑123 (MQ=255)
                                                                                       tGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATcc                                                           >  8:171939/1‑149 (MQ=255)
                                                                                              ggTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                                  >  1:287026/1‑103 (MQ=255)
                                                                                              ggTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                                  <  2:287026/103‑1 (MQ=255)
                                                                                                tGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAgg                                                  <  5:53180/149‑1 (MQ=255)
                                                                                                    cTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGAAGGCATAGGCGCAGGTCATTTCTTCTTCCGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCt                                               <  7:127566/148‑1 (MQ=255)
                                                                                                           ggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGGAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                       <  4:37084/149‑1 (MQ=255)
                                                                                                           ggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                       <  4:32556/149‑1 (MQ=255)
                                                                                                                    gTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCt                                                                      >  7:138459/1‑109 (MQ=255)
                                                                                                                    gTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCt                                                                      <  8:138459/109‑1 (MQ=255)
                                                                                                                        tACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                          >  7:322636/1‑149 (MQ=255)
                                                                                                                           aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtg                                      <  6:328135/134‑1 (MQ=255)
                                                                                                                           aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtg                                      >  5:328135/1‑134 (MQ=255)
                                                                                                                           aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAg                       >  5:254817/1‑149 (MQ=255)
                                                                                                                               tCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCt                   >  3:42451/1‑149 (MQ=255)
                                                                                                                                       aGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtgg                                        <  3:153269/120‑1 (MQ=255)
                                                                                                                                       aGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtgg                                        >  4:153269/1‑120 (MQ=255)
                                                                                                                                              gggAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCt    <  8:322636/149‑1 (MQ=255)
                                                                                                                                              gggAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCt    <  7:171939/149‑1 (MQ=255)
                                                                                                                                               ggAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTg   <  8:277498/149‑1 (MQ=255)
                                                                                                                                                gAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTATTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCttg  >  1:144498/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5518400‑5518692

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