Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,408,944 G→A 9.7% S466N (AGC→AAC)  PP_1230 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,408,9440GA 9.7% 109.9 / 2.2 41S466N (AGC→AAC) PP_1230hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

TGCAGGCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGG  >  NC_002947/1408815‑1409085
                                                                                                                                 |                                                                                                                                             
tGCAGGCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAACAGCCTGCACTATGCCGAt                                                                                                                            >  1:177579/1‑149 (MQ=255)
cgCAGGCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAACAGCCTGCACTATGCCGAt                                                                                                                            >  1:177581/2‑149 (MQ=255)
     gCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATg                                                                                                                                 >  5:192074/1‑139 (MQ=255)
     gCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATg                                                                                                                                 <  6:192074/139‑1 (MQ=255)
                       aCCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGcc                                                                                                      >  2:290268/1‑148 (MQ=255)
                             gCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAAt                                                                                               <  1:290268/149‑1 (MQ=255)
                               ccAAGCGATTGCCATGCAGTCGAGCCGGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGa                                                                                             >  1:187380/1‑149 (MQ=255)
                                 aaGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGAcc                                                                                           >  6:355455/1‑149 (MQ=255)
                                          ccATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAAt                                                                                               >  5:115652/1‑136 (MQ=255)
                                          ccATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAAt                                                                                               <  6:115652/136‑1 (MQ=255)
                                          ccATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCAATATGCCGATGGCATGGGCGATTTCAATGGCCGAAAAATGACCGTTgaagaa                                                                                  >  7:198516/1‑149 (MQ=255)
                                                 gTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCa                                                                           <  5:355455/149‑1 (MQ=255)
                                                             gggggAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATggcgg                                                               >  8:270036/1‑149 (MQ=255)
                                                             gggggAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATgg                                                                  >  2:115156/1‑146 (MQ=255)
                                                             gggggAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATgg                                                                  <  1:115156/146‑1 (MQ=255)
                                                               gggAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCATGGCGATGATGTGGTCAGCAGCCTGCAATATGCCGATGGCATGGGCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGgcggcg                                                             <  8:198516/149‑1 (MQ=255)
                                                                ggAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGgcggcg                                                             >  5:65406/1‑148 (MQ=255)
                                                                  agaTGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCa                                                           >  6:234124/1‑148 (MQ=255)
                                                                   gaTGGTCGGCATGATGGCGCTACAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATgg                                                                  <  8:267605/140‑1 (MQ=255)
                                                                   gaTGGTCGGCATGATGGCGCTACAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATgg                                                                  >  7:267605/1‑140 (MQ=255)
                                                                     tGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCt                                                        <  7:270036/148‑1 (MQ=255)
                                                                          ggCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCAGCGCATCTGGGGGc                                                  <  2:187380/149‑1 (MQ=255)
                                                                              tgatgGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTg                                              >  7:163843/1‑149 (MQ=255)
                                                                               gatgGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGt                                             >  3:144474/1‑149 (MQ=255)
                                                                                atgGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTgg                                                      >  5:54868/1‑139 (MQ=255)
                                                                                atgGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTgg                                                      <  6:54868/139‑1 (MQ=255)
                                                                                atgGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTc                                            >  6:331118/1‑149 (MQ=255)
                                                                                   gcgcTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTcgcc                                         <  6:65406/149‑1 (MQ=255)
                                                                                     gcTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTcgccg                                        >  6:316570/1‑148 (MQ=255)
                                                                                         cAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAgg                                    <  8:163843/148‑1 (MQ=255)
                                                                                         cAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAACAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAg                                     <  2:177581/147‑1 (MQ=255)
                                                                                         cAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAACAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAg                                     <  2:177579/147‑1 (MQ=255)
                                                                                                    ggCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAg                                     >  2:175523/1‑136 (MQ=255)
                                                                                                    ggCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAg                                     <  1:175523/136‑1 (MQ=255)
                                                                                                      ccACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGgcgc                                                            <  6:271992/111‑1 (MQ=255)
                                                                                                      ccACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGgcgc                                                            >  5:271992/1‑111 (MQ=255)
                                                                                                       cACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGCCCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTtctc                      >  8:147832/1‑148 (MQ=255)
                                                                                                       cACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTtctc                      <  7:147832/148‑1 (MQ=255)
                                                                                                       cACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTcc                     <  5:316570/149‑1 (MQ=255)
                                                                                                                        atgTGGTCAGCAGCCTGCAGTATGCCGATGGCATTGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAAt                                                                                 <  8:16531/72‑1 (MQ=255)
                                                                                                                           tgGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCgg  <  5:331118/148‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                             
TGCAGGCCTTGCTGGATGGCCAGACCGACGCCCAAGCGATTGCCATGCAGTCGAGCCAGGCGGGGGAGATGGTCGGCATGATGGCGCTGCAGAGCGGCATGGCCACCGTGCAGGGCGATGATGTGGTCAGCAGCCTGCACTATGCCGATGGCATGGTCGATTTCAATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGG  >  NC_002947/1408815‑1409085

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