Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,204,580 A→G 9.5% D444G (GAT→GGT)  PP_2810 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,204,5800AG 9.5% 113.5 / 2.7 42D444G (GAT→GGT) PP_2810hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (2/2);  total (22/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 = 9.95e-01

GCAAGGGCAGATCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACGAAAGCGT  >  NC_002947/3204434‑3204721
                                                                                                                                                  |                                                                                                                                             
gcaagGGCAGATCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGata                                                                                                                                             >  2:9010/1‑149 (MQ=255)
    gggCAGATCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCg                                                                                                                                         >  1:198313/1‑149 (MQ=255)
      gCAGATCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCt                                                                                                                                       <  3:98863/149‑1 (MQ=255)
          aTCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGAt                                                                                                                                   >  3:89758/1‑149 (MQ=255)
             tGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCt                                                                                                                                       >  6:63436/1‑142 (MQ=255)
             tGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCt                                                                                                                                       <  5:63436/142‑1 (MQ=255)
                        tGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTg                                                                                                                     <  2:110695/149‑1 (MQ=255)
                         gAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGaa                                                                                                                            <  7:160696/141‑1 (MQ=255)
                         gAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTAGATGGATAACCGaa                                                                                                                            >  8:160696/1‑141 (MQ=255)
                           aGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCa                                                                                                                  <  6:272985/149‑1 (MQ=255)
                              gTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAAcc                                                                                                               >  1:33005/1‑149 (MQ=255)
                                 ctcGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTc                                                                                                            >  2:206249/1‑149 (MQ=255)
                                    gATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCATACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGc                                                                                                          <  6:4552/148‑1 (MQ=255)
                                    gATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGGTATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGcc                                                                                                         >  7:157375/1‑149 (MQ=255)
                                    gATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGGTATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGcc                                                                                                         >  7:204395/1‑149 (MQ=255)
                                            tttACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAgcgc                                                                                                 <  3:210722/149‑1 (MQ=255)
                                             ttACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGaa                                                                                                                            <  4:188369/121‑1 (MQ=255)
                                             ttACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGaa                                                                                                                            >  3:188369/1‑121 (MQ=255)
                                                cGAGACGGGCTCGTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCTGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGAc                                                                                              <  1:193165/148‑1 (MQ=255)
                                                  agaCCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGActtc                                                                                           <  2:33005/149‑1 (MQ=255)
                                                                          cccAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGcc                                                                   <  6:193760/149‑1 (MQ=255)
                                                                           ccAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGATGGATAACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGAGTTCTTCACCTCGGATAACCTGCGAGCCa                                                                  >  2:232879/1‑149 (MQ=255)
                                                                              ggTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCACCCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCa                                                               >  8:128732/1‑149 (MQ=255)
                                                                              ggTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCa                                                               <  2:198313/149‑1 (MQ=255)
                                                                                    gTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCagcgagc                                                         >  2:210221/1‑149 (MQ=255)
                                                                                           ttGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTaaaa                                                  >  4:101820/1‑149 (MQ=255)
                                                                                                   tGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGGTATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTc                                          <  8:157375/149‑1 (MQ=255)
                                                                                                   tGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGGTATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTc                                          <  8:204395/149‑1 (MQ=255)
                                                                                                        gTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATc                                     <  1:164165/149‑1 (MQ=255)
                                                                                                        gTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACACTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCAGCAGCGAGCGCTACAAACGCCGTCAGAt                                      >  8:175598/1‑148 (MQ=255)
                                                                                                               ctTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGc                               >  5:122412/1‑148 (MQ=255)
                                                                                                                ttCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGActtc                                                                                           <  6:234520/87‑1 (MQ=255)
                                                                                                                ttCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGActtc                                                                                           >  5:234520/1‑87 (MQ=255)
                                                                                                                     cATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAt                         >  2:12553/1‑148 (MQ=255)
                                                                                                                           cTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGaaaaa                   >  3:77674/1‑148 (MQ=255)
                                                                                                                            tCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGaaaaaaa                 >  3:187306/1‑149 (MQ=255)
                                                                                                                             cGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAATCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGaaaaaaa                 >  6:158807/1‑148 (MQ=255)
                                                                                                                                  ggtgcaggcaAACACGATATTCGCTGGATCACCGAAGACATTGGCAACACTAGCCTGAAGCGCTACTTCTTCACCTCGTATAACCTGAGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAAtt            <  7:175598/145‑1 (MQ=255)
                                                                                                                                   gggCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTAc          >  7:214749/1‑149 (MQ=255)
                                                                                                                                    ggcaggcaAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACg         <  4:89758/149‑1 (MQ=255)
                                                                                                                                      caggcaAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTTCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACGaa       <  1:210221/149‑1 (MQ=255)
                                                                                                                                           aaaCACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACGAAAGCGt  >  8:68862/1‑149 (MQ=255)
                                                                                                                                           aaaCACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACGAAAGCGt  <  4:77674/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCAAGGGCAGATCTGGAAGGTGCGTGAAGGGTTCTCGATTCCCGTTTACGAGACCGGCGCTTGCGACGTGCAGGCCCAGGTCCAGTACAACTTGGCCGATGGTCGTTACCTCTTCGACATGACCTCGATCGGGGCAGGCAAACACGATATTCGCTGGATCACCGAAGACAGTGGCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCATGAAAAAAATAATTACGAAAGCGT  >  NC_002947/3204434‑3204721

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