Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,271,579 G→A 9.1% G161S (GGC→AGC)  PP_1112 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,271,5790GA 9.1% 123.0 / 2.6 44G161S (GGC→AGC) PP_1112membrane protein
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

ACTTCAACCGGTTCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGCCGGCCAGTGTCCGTGCCGA  >  NC_002947/1271445‑1271721
                                                                                                                                      |                                                                                                                                              
aCTTCAACCGGTTCGCGGCAACCGACAATCTTTCCCCGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGGTGCGTGACGTGGGTGGCTTCCCACTGttt                                                                                                                                  >  6:313355/1‑149 (MQ=255)
  ttCAACCGGTTCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTc                                                                                                                                >  3:149794/1‑149 (MQ=255)
       ccGGTTCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCtgt                                                                                                                           <  6:74813/149‑1 (MQ=255)
           ttCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTc                                                                                                                                <  1:103689/140‑1 (MQ=255)
           ttCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTc                                                                                                                                >  2:103689/1‑140 (MQ=255)
               cgGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGcc                                                                                                                   >  2:126777/1‑149 (MQ=255)
               cgGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGcc                                                                                                                   >  8:239313/1‑149 (MQ=255)
                   aaCCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGt                                                                                                               <  1:126777/149‑1 (MQ=255)
                        aCAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTc                                                                                                          >  4:206054/1‑149 (MQ=255)
                         cAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCg                                                                                                         <  7:239313/149‑1 (MQ=255)
                         cAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCg                                                                                                         <  3:206054/149‑1 (MQ=255)
                           aTCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGc                                                                                                       >  2:212827/1‑149 (MQ=255)
                             cTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCt                                                                                                     >  4:198024/1‑149 (MQ=255)
                                  ccAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGt                                                                                                >  7:148229/1‑149 (MQ=255)
                                  ccAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGt                                                                                                >  6:25064/1‑149 (MQ=255)
                                    aGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCcgc                                                                                                                 <  2:130963/130‑1 (MQ=255)
                                    aGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCcgc                                                                                                                 >  1:130963/1‑130 (MQ=255)
                                      gCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCt                                                                                            <  8:267251/149‑1 (MQ=255)
                                         ttGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCttt                                                                                          >  2:141132/1‑148 (MQ=255)
                                            aaGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAt                                                                                       <  3:283038/148‑1 (MQ=255)
                                               cccAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAttgtt                                                                                   <  4:149794/149‑1 (MQ=255)
                                               cccAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAttgtt                                                                                   >  2:220252/1‑149 (MQ=255)
                                               cccAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAttgtt                                                                                   <  8:12711/149‑1 (MQ=255)
                                                 cAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAt                                                                                       <  6:71176/143‑1 (MQ=255)
                                                 cAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCAt                                                                                       >  5:71176/1‑143 (MQ=255)
                                                      cgcTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTgcg                                                                            >  5:271680/1‑149 (MQ=255)
                                                           tCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCtgt                                                                        <  1:220252/148‑1 (MQ=255)
                                                             gCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGtgc                                                                      <  1:212827/148‑1 (MQ=255)
                                                                         gggcgggGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATg                                                         >  1:194206/1‑149 (MQ=255)
                                                                          ggcgggGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGc                                                        <  5:25064/149‑1 (MQ=255)
                                                                             ggggTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTAGccacccac                                                                                                                                       <  5:155360/67‑1 (MQ=255)
                                                                             ggggTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTAGccacccac                                                                                                                                       <  5:81356/67‑1 (MQ=255)
                                                                             ggggTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTAGccacccac                                                                                                                                       >  6:155360/1‑67 (MQ=255)
                                                                             ggggTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTAGccacccac                                                                                                                                       >  6:81356/1‑67 (MQ=255)
                                                                                        ggCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGATATGGCCagag                                          <  8:225080/149‑1 (MQ=255)
                                                                                        ggCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCATGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCagag                                          >  2:285683/1‑149 (MQ=255)
                                                                                                   ggCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCTCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAgc                               >  8:228830/1‑149 (MQ=255)
                                                                                                      cGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCgca                            <  1:141132/149‑1 (MQ=255)
                                                                                                            ttACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAg                      >  7:59167/1‑149 (MQ=255)
                                                                                                                 gAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGCCGGc                 <  3:198024/149‑1 (MQ=255)
                                                                                                                       cTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGCCGGCCAgtgt           >  5:124019/1‑149 (MQ=255)
                                                                                                                        tGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGc                     <  6:302986/138‑1 (MQ=255)
                                                                                                                        tGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGc                     >  5:302986/1‑138 (MQ=255)
                                                                                                                         gaaaacagtgggtagCCACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTCGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCCGTCAACCGCGCAGCAAAGCCGGCCAGTCTcc         >  4:57515/8‑149 (MQ=255)
                                                                                                                           gTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGc                     <  7:64770/135‑1 (MQ=255)
                                                                                                                                gtgggtggCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGCCGGCCAGTGTCCGTGCCGa  <  2:194206/149‑1 (MQ=255)
                                                                                                                                      ggCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCaa                                                              <  1:209466/83‑1 (MQ=255)
                                                                                                                                      ggCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCaa                                                              >  2:209466/1‑83 (MQ=255)
                                                                                                                                      |                                                                                                                                              
ACTTCAACCGGTTCGCGGCAACCGACAATCTTTCCCAGGCCTTGAAGCCCAAGGCGCTTTCGCTGGTCGTGGCGGGCGGGGTCATCGCGGCTGTCGTCGGGCCGACGCTTACCGAGTGGCTGCGTGACGTGGGTGGCTACCCACTGTTTTCGCTGTGTTATGCCGCGTTTGTCGGCCTGGCGTTGCTTTCATTGTTGATTGCGCTGTGCCTGCCCAACGATGCCGGTATGGCCAGAGCGGTCAACAGCGCAGCAAAGCCGGCCAGTGTCCGTGCCGA  >  NC_002947/1271445‑1271721

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