Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,645,727 A→G 7.8% A343A (GCA→GCG)  PP_4954 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,645,7270AG 7.8% 143.5 / 2.9 52A343A (GCA→GCG) PP_4954membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (2/2);  total (27/25)
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.18e-01

CAGGAACTGGCCCGCTCGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCACCCAGCTCACAC  >  NC_002947/5645578‑5645875
                                                                                                                                                     |                                                                                                                                                    
cAGGAACTGGCCCGCTCGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCa                                                                                                                                                      <  1:1250890/150‑1 (MQ=255)
 aGGAACTGGCCCGCTCGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTTATGGCGGCGTTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGACTGTTGCGCTCACTGCGGGCAc                                                                                                                                                     >  1:25809/1‑150 (MQ=255)
       tGGCCCGCTCGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTcgccg                                                                                                                                               >  2:408477/1‑150 (MQ=255)
                cGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTgg                                                                                                                                      <  2:982558/150‑1 (MQ=255)
                cGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTgg                                                                                                                                      <  2:857027/150‑1 (MQ=255)
                       ccTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACt                                                                                                                               >  2:301175/1‑150 (MQ=255)
                        cTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTg                                                                                                                              <  1:280837/150‑1 (MQ=255)
                              ctaCCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGc                                                                                                                        >  1:1119464/3‑150 (MQ=255)
                              cAACCGTTCGTGCTGGACAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCaa                                                                                                                                           <  2:361023/131‑1 (MQ=255)
                              cAACCGTTCGTGCTGGACAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCaa                                                                                                                                           >  1:361023/1‑131 (MQ=255)
                               aaCCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGcc                                                                                                                       >  1:120965/1‑150 (MQ=255)
                                        tGCTGGCCAGCGGTTATCCAGAAGCCCTCTAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTg                                                                                                              <  1:1072257/150‑1 (MQ=255)
                                         gCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTgc                                                                                                             <  2:970361/150‑1 (MQ=255)
                                           tGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTgccg                                                                                                           <  2:120965/150‑1 (MQ=255)
                                               cAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCGCTCGCCGCCAACCTGGCAacac                                                                                                                                <  1:1162378/125‑1 (MQ=255)
                                               cAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCGCTCGCCGCCAACCTGGCAacac                                                                                                                                >  2:1162378/1‑125 (MQ=255)
                                                      tATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCg                                                                                                <  1:1411754/150‑1 (MQ=255)
                                                        tCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGa                                                                                              <  1:999696/150‑1 (MQ=255)
                                                          cAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACa                                                                                            <  2:25809/150‑1 (MQ=255)
                                                           aGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGc                                                                                                                        >  2:860490/1‑121 (MQ=255)
                                                           aGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGc                                                                                                                        <  1:860490/121‑1 (MQ=255)
                                                              aGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCt                                                                                        >  2:884145/1‑150 (MQ=255)
                                                               gCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTg                                                                                       <  2:1361028/150‑1 (MQ=255)
                                                                   tCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTg                                                                                                                                       >  2:90608/1‑98 (MQ=255)
                                                                   tCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTg                                                                                                                                       <  1:90608/98‑1 (MQ=255)
                                                                       aGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGc                                                                               <  1:408477/150‑1 (MQ=255)
                                                                           cTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGc                                                                           >  1:102435/1‑150 (MQ=255)
                                                                           cTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGc                                                                           >  1:1198838/1‑150 (MQ=255)
                                                                           cTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCCGGCGGTGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGc                                                                           >  1:1106472/1‑150 (MQ=255)
                                                                            tCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCg                                                                          >  1:537393/1‑150 (MQ=255)
                                                                                  ccTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGAcagca                                                                    <  1:65230/150‑1 (MQ=255)
                                                                                   cTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCACACAGCAACCCGGACAGCCTGGCCGACGCCAGCGAcagcag                                                                   >  1:672509/1‑150 (MQ=255)
                                                                                            gCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGa                                                          >  1:828907/1‑150 (MQ=255)
                                                                                            gCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGa                                                          >  1:1119585/1‑150 (MQ=255)
                                                                                                cTGCGCATCCCGACCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCCACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGc                                                      >  2:718124/1‑150 (MQ=255)
                                                                                                cTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGc                                                      >  1:566203/1‑150 (MQ=255)
                                                                                                 tGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCt                                                     >  2:129292/1‑150 (MQ=255)
                                                                                                    gcATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCGCTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTcgc                                                  >  2:174645/1‑150 (MQ=255)
                                                                                                      aTCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTcgccg                                                >  2:735733/1‑150 (MQ=255)
                                                                                                          agagCGACACGGGCTGGCGCGGCTTGTTGCGCTCACTCCGGGCCCCCGCCACCACCCTGGCAACACGGGAGCGCCTGCTCAGTGCCGCCAGCGACGCGGCCAGCCTAGCCGCCGACGCGGACAGCCGCATCCCCGTCTGCTCGGCGcgct                                            >  2:352535/1‑150 (MQ=255)
                                                                                                             gCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGcgctcgc                                         <  2:102435/150‑1 (MQ=255)
                                                                                                                 cccGGCCTGGCGCTGCCTGTTGCGCTCAGGGGGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGaa                                     <  2:537393/150‑1 (MQ=255)
                                                                                                                 cccGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGaa                                     <  1:884145/150‑1 (MQ=255)
                                                                                                                  ccGGCCTGGCGCTGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAg                                    <  1:1389980/150‑1 (MQ=255)
                                                                                                                   cGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAgg                                   >  1:147820/1‑150 (MQ=255)
                                                                                                                     gCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGAc                                 >  1:291335/1‑150 (MQ=255)
                                                                                                                     gCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGAc                                 <  2:291335/150‑1 (MQ=255)
                                                                                                                     gCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGAc                                 <  1:301175/150‑1 (MQ=255)
                                                                                                                             cgGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCTGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGAc                         <  1:388987/150‑1 (MQ=255)
                                                                                                                                       cgcTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGAGGACACAATCATCACCAGACTGGTGAACTGCTCGCCGCGCTCGCTGAAGGACGTATGGCCGCGCCGgcgc               >  1:525399/1‑150 (MQ=255)
                                                                                                                                          tCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCAcc            >  2:772867/1‑150 (MQ=255)
                                                                                                                                              tGCGGGCGCTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGAAGCCAGCGACAGCAGCCTTCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCACCCAGc        <  1:174645/150‑1 (MQ=255)
                                                                                                                                                  ggaaCTCGCAGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCACCCAGCTcac    >  2:878084/4‑150 (MQ=255)
                                                                                                                                                    cACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCACCCAGCTCacac  <  1:772867/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGAACTGGCCCGCTCGATCCGCCTGCGCCAACCGTTCGTGCTGGCCAGCGGTTATCCAGAAGCCCTCGAAGACCTCAACGCCTCGCTGGAGCACCTGCGCATCCAGAGCAACCCGGCCTGGCGCGGCCTGTTGCGCTCACTGCGGGCACTCGCCGCCAACCTGGCAACACTGGACCGCCTGCTCAGTGCCGCCAGCAACCCGGACAGCCTGGCCGACGCCAGCGACAGCAGCCTGCTCGACCGCTCGCCGCGCTCGCTGAAGGACGTGTGGACGCGCCTGCGCACCCAGCTCACAC  >  NC_002947/5645578‑5645875

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