Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,194,356 C→G 10.0% P226R (CCC→CGC)  PP_4574 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,194,3560CG10.0% 114.5 / 5.0 40P226R (CCC→CGC) PP_4574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/21);  new base G (2/2);  total (17/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 = 8.70e-01

GTCGAGCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAT  >  NC_002947/5194211‑5194504
                                                                                                                                                  |                                                                                                                                                    
gTCGAGCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGcccc                                                                                                                                                     >  6:27742/1‑147 (MQ=255)
gTCGAGCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGcccc                                                                                                                                                     >  2:54691/1‑147 (MQ=255)
     gCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAgcggc                                                                                                                                              <  1:359265/149‑1 (MQ=255)
         ggCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGgcgc                                                                                                                                            <  5:27742/147‑1 (MQ=255)
         ggCACCCGCAACATCAAGCTGGCCCTGCGCAGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCaa                                                                                                                                          >  3:78083/1‑149 (MQ=255)
          gCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAac                                                                                                                                         <  1:54691/149‑1 (MQ=255)
               cGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGt                                                                                                                                    <  4:118398/149‑1 (MQ=255)
                     aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCa                                                                                                                                           >  1:69282/1‑137 (MQ=255)
                     aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCa                                                                                                                                           <  2:69282/136‑1 (MQ=255)
                     aTCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAgct                                                                                                                              >  7:341643/1‑149 (MQ=255)
                        aaGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgga                                                                                                                                     <  5:188150/139‑7 (MQ=255)
                        aaGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgga                                                                                                                                     >  6:188150/1‑133 (MQ=255)
                              gCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCtgctgttgc                                                                                                                      <  2:10506/148‑1 (MQ=255)
                                 cTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTgctgttgctgtt                                                                                                                  <  3:283202/149‑1 (MQ=255)
                                    cgcCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCg                                                                                                                <  1:222773/148‑1 (MQ=255)
                                       cGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACAt                                                                                                            >  3:328114/1‑149 (MQ=255)
                                         gCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTg                                                                                                          <  1:144498/149‑1 (MQ=255)
                                              gCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACa                                                                                                             <  8:47985/141‑1 (MQ=255)
                                              gCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACa                                                                                                             >  7:47985/1‑141 (MQ=255)
                                              gCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGt                                                                                                     >  5:178852/1‑149 (MQ=255)
                                               cAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    <  8:109355/149‑1 (MQ=255)
                                                   ttCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgg                                                                                                                                      <  7:38062/111‑6 (MQ=255)
                                                   ttCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCGCTCGGGGCGCatctgg                                                                                                                                      >  8:38062/1‑106 (MQ=255)
                                                                ggCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGt                                                                                    <  4:78083/149‑1 (MQ=255)
                                                                       ccgAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTc                                                                             <  8:187203/149‑1 (MQ=255)
                                                                                  gACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCAc                                                                                      >  8:217172/1‑129 (MQ=255)
                                                                                  gACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCAc                                                                                      <  7:217172/129‑1 (MQ=255)
                                                                                    cATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    <  1:190115/113‑1 (MQ=255)
                                                                                    cATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGtt                                                                                                    >  2:190115/1‑113 (MQ=255)
                                                                                        gATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            >  4:350653/1‑149 (MQ=255)
                                                                                            gCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGc                                                         <  6:178852/148‑1 (MQ=255)
                                                                                                 aTCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCaaa                                                                                <  4:156359/120‑1 (MQ=255)
                                                                                                    gACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGAc                                                >  5:17759/1‑149 (MQ=255)
                                                                                                     aCCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGAcc                                               >  5:146285/1‑149 (MQ=255)
                                                                                                              cgcgTGATTCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTTTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAg                                      <  5:149597/149‑1 (MQ=255)
                                                                                                                        cGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAg                             >  1:28020/1‑148 (MQ=255)
                                                                                                                              gttgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGc                                                     <  2:241502/118‑1 (MQ=255)
                                                                                                                              gttgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGc                                                     >  1:241502/1‑118 (MQ=255)
                                                                                                                                tgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            >  6:72551/1‑109 (MQ=255)
                                                                                                                                tgAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTtctc                                                            <  5:72551/109‑1 (MQ=255)
                                                                                                                                 gAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTactact                   >  3:46747/1‑149 (MQ=255)
                                                                                                                                   aCATCCAGATGCGCCCCGAGCGGCGCAACTCCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGTCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACtt                  <  8:341643/148‑1 (MQ=255)
                                                                                                                                                  ccGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAt  >  1:50636/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCGAGCTTGGCACCCGCAACATCAAGCTGGCCCTGCGCCGGCTACGCAAGTTCGCCCGCGAAGGCGCCGCCGAAGAGCTCGACATCGATGGCACCATCGACCACACCGCGCGTGATGCCGGGTTGTTGAACATCCAGATGCGCCCCGAGCGGCGCAACACCGTGAAGCTGCTGTTGCTGTTCGACATTGGCGGTTCGATGGACGCCCACGTCAAAGTCTGCGAAGAGCTGTTCTCGGCCTGCAAGACCGAGTTCAAGCACCTGGAGTATTACTACTTCCACAATTGCGTTTAT  >  NC_002947/5194211‑5194504

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