Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,415 G→A 9.7% G194S (GGC→AGC)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,4150GA 9.7% 111.7 / 2.4 41G194S (GGC→AGC) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/22);  new base A (2/2);  total (17/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 = 1.00e+00

CCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAAT  >  NC_002947/4301271‑4301560
                                                                                                                                                |                                                                                                                                                 
ccACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCAcggcg                                                                                                                                                <  5:119740/148‑1 (MQ=255)
  aCTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGCCAGTACCTGGGCTTCACGCCTGTGCGAGTTCCACGGCgggg                                                                                                                                             >  8:319862/1‑149 (MQ=255)
   cTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCAGGCCTGTTCGAGTTCCACGGCGGGGt                                                                                                                                            <  2:359015/149‑1 (MQ=255)
      ggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCGACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTc                                                                                                                                           <  7:169047/147‑1 (MQ=255)
      ggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTc                                                                                                                                           <  3:163468/147‑1 (MQ=255)
      ggTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTa                                                                                                                                         >  5:822/1‑149 (MQ=255)
                     aTCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTa                                                                                                                          >  2:338868/1‑149 (MQ=255)
                     aTCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGGCGACTa                                                                                                                          >  8:325190/1‑149 (MQ=255)
                         gCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGTGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACAGCGGGGTCTACAACCAGGTCGACTACAcc                                                                                                                      >  4:62298/1‑149 (MQ=255)
                              cacaCCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACg                                                                                                                  <  3:242614/148‑1 (MQ=255)
                               acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGAc                                                                                                                   <  6:329585/146‑1 (MQ=255)
                               acacCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGAc                                                                                                                   >  5:329585/1‑146 (MQ=255)
                                caccacGTCTGAAGGCGTGGTAACCGTGACGGCCAAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCaa                                                                                                               <  7:2513/149‑1 (MQ=255)
                                  ccacGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGTGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACAGCGGGGTCTACAACCAGGTCGACTACACCGACGCAAcc                                                                                                             <  3:62298/149‑1 (MQ=255)
                                   cacGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAAccc                                                                                                            <  1:152567/149‑1 (MQ=255)
                                   cacGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAAccc                                                                                                            <  1:158048/149‑1 (MQ=255)
                                             ggCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGa                                                                                                  <  1:338868/149‑1 (MQ=255)
                                                gTGGTAACCGTGCCGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTTCACCGACGCAACCCCCGCCGAGGACgg                                                                                               >  1:64937/1‑149 (MQ=255)
                                                   gTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTc                                                                                            <  7:325190/149‑1 (MQ=255)
                                                    tAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCg                                                                                           <  3:164393/149‑1 (MQ=255)
                                                      aCCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGAt                                                                                         <  2:214891/149‑1 (MQ=255)
                                                        cGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACAGCGGGGTCTACAACCAGGTCGACTaca                                                                                                                        <  4:117001/116‑1 (MQ=255)
                                                        cGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACAGCGGGGTCTACAACCAGGTCGACTaca                                                                                                                        >  3:117001/1‑116 (MQ=255)
                                                         gTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTg                                                                                      >  3:86666/1‑149 (MQ=255)
                                                          tGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGc                                                                                     <  5:129931/149‑1 (MQ=255)
                                                             cGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGcatc                                                                                  <  7:304668/149‑1 (MQ=255)
                                                                 cGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGcatcat                                                                                >  6:42862/1‑147 (MQ=255)
                                                                          cGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGCTCTGCATCATCAGCGAGctct                                                                     >  8:176355/1‑149 (MQ=255)
                                                                             ggTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGc                                                                  <  4:51707/149‑1 (MQ=255)
                                                                                 tCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTca                                                               <  5:42862/148‑1 (MQ=255)
                                                                                   cGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTcatc                                                             >  5:197501/1‑148 (MQ=255)
                                                                                                    tGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGAt                                            <  6:822/148‑1 (MQ=255)
                                                                                                      cAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGcc                                         >  5:259752/1‑149 (MQ=255)
                                                                                                                  aacTGTGCTTCCCGCCTGTTCGTTTTCCACGGCTGGGGCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCTAGCTCTCGCTCTTCCTTCGCTGAACGTCGATGGCCAGCCCTGCTa                              <  4:86666/146‑1 (MQ=255)
                                                                                                                    cTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGAt                                            <  4:170326/132‑1 (MQ=255)
                                                                                                                    cTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGAt                                            >  3:170326/1‑132 (MQ=255)
                                                                                                                    cTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATc                            >  7:178208/1‑148 (MQ=255)
                                                                                                                               gCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAg                >  8:322893/1‑149 (MQ=255)
                                                                                                                                  tGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAGCGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGc             >  2:100487/1‑149 (MQ=255)
                                                                                                                                       gAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCg         <  8:178208/148‑1 (MQ=255)
                                                                                                                                         gTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCAt      <  1:143910/149‑1 (MQ=255)
                                                                                                                                             cACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGACCGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAAt  <  1:100487/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCACTGGGTGCGACGTTCAACATCAGCATTCACACCACGTCTGAAGGCGTGGTAACCGTGACGGCCGAATGCAACGGGGTATCCGCTTCGTCGGGCGATTTGCAGCTCGACAGTACCTGGGCTTCACGCCTGTTCGAGTTCCACGGCGGGGTCTACAACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAAT  >  NC_002947/4301271‑4301560

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