Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,015,332 G→A 9.5% A303A (GCC→GCT)  PP_0875 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,015,3320GA 9.5% 115.6 / 2.7 42A303A (GCC→GCT) PP_0875hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (2/2);  total (21/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.06e-01

AGTGATGCGTACCTGTTTGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTGGCAATTGGTGTTTGCCCAGTGG  >  NC_002947/1015185‑1015478
                                                                                                                                                   |                                                                                                                                                  
aGTGATGCGTACCTGTTTGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCgg                                                                                                                                                   <  1:129831/149‑1 (MQ=255)
aGTGATGCGTACCTGTTTGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCg                                                                                                                                                    <  3:277601/148‑1 (MQ=255)
           ccTGTTTGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGAGCAGACGGTTGTCGGCAGTGCGCTg                                                                                                                                         <  4:175152/148‑1 (MQ=255)
               tttGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGcc                                                                                                                                      <  8:111331/147‑1 (MQ=255)
                        cgtgAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGAc                                                                                                                             <  4:247668/144‑1 (MQ=255)
                                 ttCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGCTAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGgctgc                                                                                                                  <  2:210297/149‑1 (MQ=255)
                                   cGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCtt                                                                                                                >  5:65710/1‑149 (MQ=255)
                                         aCAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACt                                                                                                           >  5:155579/1‑148 (MQ=255)
                                             gTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAgg                                                                                                      >  2:161421/1‑149 (MQ=255)
                                                     gcAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCCTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGccc                                                                                              >  8:158419/1‑149 (MQ=255)
                                                      cAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCt                                                                                             <  6:65710/149‑1 (MQ=255)
                                                         ttCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGc                                                                                          >  4:276145/1‑149 (MQ=255)
                                                                 cGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGGAGGCTTGGCCCTTGCTCAGGcc                                                                                   >  2:53945/1‑148 (MQ=255)
                                                                  gtgtTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAg                                                                                                                                  >  8:241149/1‑100 (MQ=255)
                                                                  gtgtTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAg                                                                                                                                  <  7:241149/100‑1 (MQ=255)
                                                                   tgtTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAAt                                                                                <  6:86154/149‑1 (MQ=255)
                                                                    gtTCGTCGACGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATATCCACATCGAGCACCGAAGACCAGACGGTTGTCGCCAGTGCGCTGGCCGTAGGTGACCTGGCGCCTGATTTCACAGAAGGCTTGGCCCTTGCTCAGGCCAAtt                                                                               >  2:218924/1‑149 (MQ=255)
                                                                       cgtcgAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTg                                                                             <  1:158997/148‑1 (MQ=255)
                                                                         tcgAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGCCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTgg                                                                                                 <  5:80610/126‑1 (MQ=255)
                                                                         tcgAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGCCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTgg                                                                                                 >  6:80610/1‑126 (MQ=255)
                                                                          cgAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCAGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTcc                                                                         >  8:184169/1‑149 (MQ=255)
                                                                          cgAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCAGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTcc                                                                         >  8:179379/1‑149 (MQ=255)
                                                                                           cgcAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGGAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAg                                                        <  1:53945/149‑1 (MQ=255)
                                                                                                     cACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGa                                                                                                         <  2:302597/90‑1 (MQ=255)
                                                                                                     cACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGa                                                                                                         >  1:302597/1‑90 (MQ=255)
                                                                                                     cACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTg                                              >  5:218868/1‑149 (MQ=255)
                                                                                                     cACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCAGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTg                                              <  7:184169/149‑1 (MQ=255)
                                                                                                     cACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCAGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTg                                              <  7:179379/149‑1 (MQ=255)
                                                                                                        cGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGc                                            >  6:227764/1‑148 (MQ=255)
                                                                                                         gAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGcca                                          >  2:117813/1‑149 (MQ=255)
                                                                                                          aaGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGccac                                         <  5:227764/149‑1 (MQ=255)
                                                                                                            gCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCg                                                 >  2:196883/1‑139 (MQ=255)
                                                                                                            gCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCg                                                 <  1:196883/139‑1 (MQ=255)
                                                                                                            gCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGccacca                                       <  1:161421/149‑1 (MQ=255)
                                                                                                            gCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGccac                                         >  3:211426/1‑147 (MQ=255)
                                                                                                               aTAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACcag                                      >  2:136452/1‑147 (MQ=255)
                                                                                                                    cACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGctc                               >  6:102379/1‑149 (MQ=255)
                                                                                                                          gAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACt                         <  6:218868/149‑1 (MQ=255)
                                                                                                                            gCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTgg                       <  1:218924/149‑1 (MQ=255)
                                                                                                                            gCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTgg                       <  4:211426/149‑1 (MQ=255)
                                                                                                                                       cAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAgg                                  <  1:198439/127‑1 (MQ=255)
                                                                                                                                       cAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAgg                                  >  2:198439/1‑127 (MQ=255)
                                                                                                                                            ggTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTGGCAATTGGTGTTTGcc        >  6:212413/1‑148 (MQ=255)
                                                                                                                                                 tCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTGGCAATTGGTGTTTGCCCAGTgg  <  5:212413/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTGATGCGTACCTGTTTGAGCTGCGGGAACAGTTCGCCGAACAGGTAGCGGCGCAGTTCGACTTCGTGTTCGTCGAGGTTGAAGTTCTCGCGCAGACGCCCACCGAAGCGATAGCCACCGCGAGCACCGAAGACCAGACGGTTGTCGGCAGTGCGCTGGCCGTAGGTGACCTGGCGGCTGCTTTCACTGAAGGCTTGGCCCTTGCTCAGGCCAATTTGCTCCCAGGTGGCCTCTGGCAGCGGCTCGGTGGCCACCAGCAGGCTCTGTACTGGCAATTGGTGTTTGCCCAGTGG  >  NC_002947/1015185‑1015478

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