Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,101,367 C→G 8.1% G74A (GGC→GCC)  PP_3609 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,101,3670CG 8.1% 151.7 / 2.7 49G74A (GGC→GCC) PP_3609membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base G (2/2);  total (25/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 = 9.73e-01

CCGTGATGATCTGCCCCGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCG  >  NC_002947/4101227‑4101514
                                                                                                                                             |                                                                                                                                                   
ccGTGATGATCTGCCCCGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGt                                                                                                                                               <  6:141159/147‑1 (MQ=255)
   tgatgaTCTGCCCCGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGgc                                                                                                                                          >  6:351849/1‑149 (MQ=255)
             ccccGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCCTCGCCGGCTACATACATTGCCCCGAGCCCGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTAccc                                                                                                                                >  1:278137/1‑149 (MQ=255)
             ccccGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCATCGCCGGCTACATACATTGCCCCGCGCACGCCACCAATCCCCAACCACCAAGGCCCCTGCC‑GGGCCTTGCCCATGTCTGGCGCCGGTcccc                                                                                                                                >  1:282328/1‑149 (MQ=255)
                cGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGa                                                                                                                             <  1:145127/149‑1 (MQ=255)
                 gCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGAt                                                                                                                            <  8:235132/149‑1 (MQ=255)
                              ccAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCa                                                                                                               <  7:287600/149‑1 (MQ=255)
                                 ggAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGc                                                                                                             >  2:136185/1‑148 (MQ=255)
                                    aCCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                         <  1:136185/149‑1 (MQ=255)
                                    aCCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                         >  5:70212/1‑149 (MQ=255)
                                          atgcGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCACGAGCACGCCACCAAACCACAACCACCAAGTCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATccacac                                                                                                    >  2:164113/2‑143 (MQ=255)
                                             cgcCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCAccacca                                                                                                  <  5:351849/147‑1 (MQ=255)
                                                    tttAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCa                                                                                         >  6:95955/1‑149 (MQ=255)
                                                          ggTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCAGCAATAGACAACCACCAAGGCCGCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGAGCACCACCAGCATGCTCACCGTCa                                                                                   <  1:164113/149‑1 (MQ=255)
                                                           gTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCaa                                                                                  >  4:336708/1‑149 (MQ=255)
                                                                          ccGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGcccc                                                                   <  6:70212/149‑1 (MQ=255)
                                                                              cTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCca                                                                >  1:98849/1‑149 (MQ=255)
                                                                                  atacatTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGcc                                                           >  8:283010/1‑149 (MQ=255)
                                                                                          ccccGAGCACGCCACCACTCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCCCCAGCCTGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGccc                                                    <  2:282328/148‑1 (MQ=255)
                                                                                          ccccGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGcac                                                    <  2:278137/148‑3 (MQ=255)
                                                                                             cGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGt                                                <  2:98849/149‑1 (MQ=255)
                                                                                                    gCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                         <  7:89988/85‑1 (MQ=255)
                                                                                                    gCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                         >  8:89988/1‑85 (MQ=255)
                                                                                                      caccaATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGccc                                                                                                                                                   <  2:193301/41‑1 (MQ=255)
                                                                                                      caccaATCCACAACCACCAAGGCCCCGGCA‑GGGCCTTGccc                                                                                                                                                   >  1:193301/1‑41 (MQ=255)
                                                                                                         caATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGcc                                    >  2:156205/1‑149 (MQ=255)
                                                                                                                caACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAcc                                        <  4:47436/138‑1 (MQ=255)
                                                                                                                caACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAcc                                        >  3:47436/1‑138 (MQ=255)
                                                                                                                caACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCt                                       >  8:81318/1‑139 (MQ=255)
                                                                                                                caACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCt                                       <  7:81318/139‑1 (MQ=255)
                                                                                                                   ccaccaAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTAc                          <  5:95955/149‑1 (MQ=255)
                                                                                                                   acaccaAGGCCCCTTCA‑GGGCCTTGCCCATGTCTGGCGCCTGTACCCCGATAATCAGCATGGCAGCGATCACCTCCCTCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTAc                          <  8:313749/148‑1 (MQ=255)
                                                                                                                     ggcaAGG‑CCCTGCAGGGGCCTTGGCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAc                                         <  2:336527/130‑1 (MQ=255)
                                                                                                                     ggcaAGG‑CCCTGCAGGGGCCTTGGCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAc                                         <  2:331313/130‑1 (MQ=255)
                                                                                                                     ggcaAGG‑CCCTGCAGGGGCCTTGGCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAc                                         >  1:331313/3‑132 (MQ=255)
                                                                                                                     ggcaAGG‑CCCTGCAGGGGCCTTGGCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGAc                                         >  1:336527/3‑132 (MQ=255)
                                                                                                                     accaAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCAGGCTCCCCGTCAACGAGGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCt                        >  7:132521/1‑149 (MQ=255)
                                                                                                                      gcaAGG‑CCCTGCAGGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTg                       >  4:35690/2‑149 (MQ=255)
                                                                                                                       caAGG‑CCCTGCAGGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCgcagca                              >  2:244057/1‑141 (MQ=255)
                                                                                                                       caAGG‑CCCTGCAGGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCgcagca                              <  1:244057/141‑1 (MQ=255)
                                                                                                                         aaggcccTGCAGGGGCCTTGCCCATGTCTGTCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCAACCCACAGCCAGTGCCCAAGTGCCCGAcc                                        >  2:202632/5‑130 (MQ=255)
                                                                                                                         aaggcccTGCAGGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCGCGAcc                                        <  1:202632/126‑1 (MQ=255)
                                                                                                                          ggCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGc                   <  7:283010/149‑1 (MQ=255)
                                                                                                                            ccccTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCt                 >  7:240994/1‑149 (MQ=255)
                                                                                                                            ccccTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCt                 <  8:10006/149‑1 (MQ=255)
                                                                                                                                 gCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCATCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCgg                    <  6:262334/141‑1 (MQ=255)
                                                                                                                                 gCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCATCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCgg                    >  5:262334/1‑141 (MQ=255)
                                                                                                                                      ggCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCCTGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGc         >  5:327710/1‑148 (MQ=255)
                                                                                                                                        ccTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCAtt                            <  6:351929/127‑1 (MQ=255)
                                                                                                                                        ccTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCAtt                            >  5:351929/1‑127 (MQ=255)
                                                                                                                                            gCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAAcaccg  >  8:215468/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CCGTGATGATCTGCCCCGCCACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCA‑GGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCG  >  NC_002947/4101227‑4101514

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