Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,277,210 G→C 11.1% A52G (GCC→GGC)  PP_2006 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,277,2100GC11.1% 104.8 / 5.6 36A52G (GCC→GGC) PP_2006membrane protein
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base C (2/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CGACCCCGGCGTCGATGTTGGCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGCATGACCTTGCTTTCATT  >  NC_002947/2277064‑2277358
                                                                                                                                                  |                                                                                                                                                    
cGACCCCGGCGTCGATGTTGGCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTgcgcggcgc                                                                                                                                                   <  1:486437/150‑1 (MQ=255)
      cGGCGTCGATGTTGGCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAgg                                                                                                                                             <  1:536907/150‑1 (MQ=255)
           tCGATGTTGGCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCt                                                                                                                                        >  2:448690/1‑150 (MQ=255)
                   ggCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACt                                                                                                                                <  2:733734/150‑1 (MQ=255)
                             gTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAAc                                                                                                                      >  2:661048/1‑150 (MQ=255)
                                      aGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTg                                                                                                             >  1:750478/1‑150 (MQ=255)
                                                                 aCCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTc                                                                                  >  1:701596/1‑150 (MQ=255)
                                                                   ccGGCAGATGCTCACGGTATCCCAGCGTTCCGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTg                                                                                >  1:157361/1‑150 (MQ=255)
                                                                   ccGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTg                                                                                >  2:686504/1‑150 (MQ=255)
                                                                   ccGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTg                                                                                >  1:410251/1‑150 (MQ=255)
                                                                    cGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTgg                                                                               >  2:123462/1‑150 (MQ=255)
                                                                              tCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGttt                                                                     >  2:104361/1‑150 (MQ=255)
                                                                                  ggTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCAtt                                                                 <  1:448690/150‑1 (MQ=255)
                                                                                  ggTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCAtt                                                                 <  2:701596/150‑1 (MQ=255)
                                                                                    tATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGTGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAggg                                                                                                   >  2:430547/1‑114 (MQ=255)
                                                                                    tATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGTGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAggg                                                                                                   <  1:430547/114‑1 (MQ=255)
                                                                                      tCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAgg                                                             >  2:30372/1‑150 (MQ=255)
                                                                                          aGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACt                                                                                                                                >  1:471847/1‑79 (MQ=255)
                                                                                          aGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACt                                                                                                                                <  2:471847/79‑1 (MQ=255)
                                                                                           gCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATc                                                        <  1:453042/150‑1 (MQ=255)
                                                                                           gCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATc                                                        <  1:661048/150‑1 (MQ=255)
                                                                                                     gTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGc                                              <  1:104361/150‑1 (MQ=255)
                                                                                                           gTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATg                                        <  2:439218/150‑1 (MQ=255)
                                                                                                                ggCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTg                                   >  2:193273/1‑150 (MQ=255)
                                                                                                                  ccAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTggg                                 >  1:339173/1‑150 (MQ=255)
                                                                                                                   cAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGc                                >  1:103109/1‑150 (MQ=255)
                                                                                                                       gtTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAAc                            <  2:189385/150‑1 (MQ=255)
                                                                                                                          ccGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCg                         <  2:157361/150‑1 (MQ=255)
                                                                                                                           cGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCgg                        >  2:720700/1‑150 (MQ=255)
                                                                                                                            ggCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGa                       >  1:114952/1‑150 (MQ=255)
                                                                                                                                   gTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGCATg                >  1:90704/1‑150 (MQ=255)
                                                                                                                                       gTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGc                   >  2:271707/1‑143 (MQ=255)
                                                                                                                                       gTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGc                   <  1:271707/143‑1 (MQ=255)
                                                                                                                                        accttGCGCGCCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTggg                                 >  1:63216/5‑128 (MQ=255)
                                                                                                                                        accttGCGCGCCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTggg                                 <  2:63216/124‑1 (MQ=255)
                                                                                                                                          cttGCGCGCCGCGCAAGGTGTCTTCGGAACTGTGTACAAAc                                                                                                                      <  1:489574/39‑1 (MQ=255)
                                                                                                                                          cttGCGCGCCGCGCAAGGTGTCTTCGGAACTGTGTACAAAc                                                                                                                      >  2:489574/3‑41 (MQ=255)
                                                                                                                                               gcggcgcgCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGCATGACCTTGCTTTCa    <  1:193273/150‑1 (MQ=255)
                                                                                                                                                 ggcgcgCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGCATGACCTTGCTTTCAtt  <  2:103109/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGACCCCGGCGTCGATGTTGGCGGGGTAGGTTTTTTCCAGGGCGATGCGCTGGTTGAGCGCAGAGACCCGGCAGATGCTCACGGTATCCCAGCCATTGGTCGTGACCGTACTGGCCAGTGTTCCGGCATTGGTCAGTAGCTGCGCGGCGCGCAAGGTGTCTTCGGAACTGTGTACAAACGCCGCGTTGCGGTACAGGGGCGATTCGCGCAGGGTCTGGCGCGCGGTTTCATTGAGGTGATCGAGGCGCAGCGAGATGGGTTGGGCAAACCCGGACAGCATGACCTTGCTTTCATT  >  NC_002947/2277064‑2277358

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: