Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,771,705 A→G 8.9% D818G (GAT→GGT)  PP_1579 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,771,7050AG 8.9% 121.4 / 2.9 46D818G (GAT→GGT) PP_1579hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (2/2);  total (23/23)
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.33e-01

ATGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAGCGGTTTGATCGCATGG  >  NC_002947/1771558‑1771851
                                                                                                                                                   |                                                                                                                                                  
aTGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGAt                                                                                                                                                   <  1:22704/149‑1 (MQ=255)
aTGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGAt                                                                                                                                                   <  2:268045/149‑1 (MQ=255)
 tGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGACGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATg                                                                                                                                                  <  1:263504/149‑1 (MQ=255)
 tGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATg                                                                                                                                                  <  5:368823/149‑1 (MQ=255)
          cgacgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGa                                                                                                                                                    >  1:227754/1‑138 (MQ=255)
          cgacgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGa                                                                                                                                                    <  2:227754/138‑1 (MQ=255)
          cgacgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATg                                                                                                                                                  <  2:35564/140‑1 (MQ=255)
          cgacgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATg                                                                                                                                                  >  1:35564/1‑140 (MQ=255)
            acgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTg                                                                                                                                       >  7:36130/1‑149 (MQ=255)
             cgaGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTgg                                                                                                                                      <  3:317502/149‑1 (MQ=255)
                gCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACg                                                                                                                                   <  4:130030/149‑1 (MQ=255)
                      cGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTTGACGGTGTTc                                                                                                                             <  5:139217/149‑1 (MQ=255)
                       gCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCt                                                                                                                            >  8:120498/1‑149 (MQ=255)
                          tCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGa                                                                                                                         >  8:215438/1‑149 (MQ=255)
                               ctcTGCTCGTGGCCAGGATTTCTTAGTCGGGGTTGACGTGAAGTACGTCGTCCCGTGCACGTTGTCTGTAGAAACGGTGCGGTGGGGCTTGGGTGAATGTCTTGGGGTGTTGGTCGTTGGTATCTGGCTGGACGGTGTTCTTGATggcg                                                                                                                    <  8:272616/147‑1 (MQ=255)
                                 ttattgcGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGAc                                                                                                                  <  1:126539/143‑1 (MQ=255)
                                                aTGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATc                                                                                                   >  2:311961/1‑149 (MQ=255)
                                                   aCTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGc                                                                                                <  7:215438/149‑1 (MQ=255)
                                                              gTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGtttt                                                                                     <  8:36130/149‑1 (MQ=255)
                                                                    gTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTaa                                                                               <  1:311961/149‑1 (MQ=255)
                                                                     tGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCGTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAg                                                                              >  1:295431/1‑149 (MQ=255)
                                                                            cGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACACTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGt                                                                       >  2:71014/1‑149 (MQ=255)
                                                                             gTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTc                                                                      >  6:350046/1‑149 (MQ=255)
                                                                                  cccccGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATcgc                                                                                                                     <  8:307201/97‑1 (MQ=255)
                                                                                  cccccGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATcgc                                                                                                                     >  7:307201/1‑97 (MQ=255)
                                                                                          ccTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCg                                                         <  2:295431/149‑1 (MQ=255)
                                                                                                 ggAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGGTGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            >  1:219897/1‑107 (MQ=255)
                                                                                                 ggAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGGTGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            <  8:229103/107‑1 (MQ=255)
                                                                                                 ggAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGGTGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            >  7:229103/1‑107 (MQ=255)
                                                                                                 ggAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGGTGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            <  2:219897/107‑1 (MQ=255)
                                                                                                  gaagaaACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            <  4:49100/106‑1 (MQ=255)
                                                                                                  gaagaaACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATgc                                                                                            >  3:49100/1‑106 (MQ=255)
                                                                                                       aaCCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGcacca                                            >  5:262622/1‑149 (MQ=255)
                                                                                                        aCCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGc                                                                                   >  5:88664/1‑109 (MQ=255)
                                                                                                        aCCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGc                                                                                   <  6:88664/109‑1 (MQ=255)
                                                                                                              ccgtgccgCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCt                                     <  7:120498/149‑1 (MQ=255)
                                                                                                                    cgCTTGGGAGAACGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGCTCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGgt                               >  6:159877/1‑149 (MQ=255)
                                                                                                                           gagaAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGGGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACAGCCCTGGAGgtgctgtgc                        >  8:284340/1‑149 (MQ=255)
                                                                                                                           gagaAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTaa                                                                               >  1:189783/1‑94 (MQ=255)
                                                                                                                           gagaAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTaa                                                                               <  2:189783/94‑1 (MQ=255)
                                                                                                                           gagaAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGgtgctgtg                         >  2:19726/1‑148 (MQ=255)
                                                                                                                              aaaGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTAt                     >  6:271774/1‑149 (MQ=255)
                                                                                                                                 gCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAg                  <  6:262622/149‑1 (MQ=255)
                                                                                                                                  ccTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTc                                                    <  7:310441/114‑1 (MQ=255)
                                                                                                                                  ccTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTc                                                    >  8:310441/1‑114 (MQ=255)
                                                                                                                                           gAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAGCGGTTTGATc        >  8:222237/1‑149 (MQ=255)
                                                                                                                                               gTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAGCGGTTTGATCGCAt    >  5:259773/1‑149 (MQ=255)
                                                                                                                                                 cGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAGCGGTTTGATCGCATgg  >  3:184744/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGAGCTGACCGACGAGCTGACCGCCTCGTTCACTGCTCCGGGCCAGGATGACTTCGACGGGGTGGACGTGAAGTACGTCGACCCCCGCACCTGGTCGGAAGAAACCGTGCCGTGCCGCTTGGGAGAAAGCCTTGGCCTGAAGGTCGATGCAATCACCCTGGACGGTGTTCTTGATCGCGACCGCGCCTGGCGGATCGGCATGCGCGTTTTGCGTAAGCATCAGTCGCAGCGCTGGTCGTACTCGGGCACCACCGACCTGGAGGTGCTGTGCTATGAGCGGTTTGATCGCATGG  >  NC_002947/1771558‑1771851

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