Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,632,128 G→T 8.7% E327* (GAG→TAG) ‡ PP_5590 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,632,1280GT 8.7% 134.6 / 3.3 47E327* (GAG→TAG) ‡PP_5590hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base T (2/2);  total (23/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 = 4.97e-01

CCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCT  >  NC_002947/3631981‑3632274
                                                                                                                                                   |                                                                                                                                                  
ccGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                  <  1:769887/150‑1 (MQ=255)
      ccGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCgctg                                                                                                                                            >  1:234073/1‑150 (MQ=255)
      ccGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCAACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCgctg                                                                                                                                            >  2:214756/1‑150 (MQ=255)
           cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                         <  1:529360/148‑1 (MQ=255)
           cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                         >  2:529360/1‑148 (MQ=255)
           cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTAAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGc                                                                                                                                       >  1:249667/1‑150 (MQ=255)
                       cATGAGCTGCGATGAGCTGTACGTGTTCACCGGTCGCGTCGACCGCTTCTAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCt                                                                                                                           <  2:567325/150‑1 (MQ=255)
                                      gCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCAt                                                                                                            >  2:482858/1‑150 (MQ=255)
                                       cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGc                                                                                                                                       <  1:693834/122‑1 (MQ=255)
                                       cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGc                                                                                                                                       >  2:693834/1‑122 (MQ=255)
                                       cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATc                                                                                                           >  1:453036/1‑150 (MQ=255)
                                               gtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGa                                                                                                   <  2:232865/150‑1 (MQ=255)
                                                             tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGACGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACGgtgtgt                                                                                     >  2:82770/1‑150 (MQ=255)
                                                                   gCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtgggtgg                                                                               >  2:521878/1‑150 (MQ=255)
                                                                    cTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGggtggt                                                                              >  2:202992/1‑150 (MQ=255)
                                                                     ttCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGgtggtg                                                                             >  2:404376/1‑150 (MQ=255)
                                                                       cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGa                                                                           <  1:521878/150‑1 (MQ=255)
                                                                       cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGa                                                                           <  1:202992/150‑1 (MQ=255)
                                                                        gAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCGCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGGGGAg                                                                          >  1:127309/1‑150 (MQ=255)
                                                                         aGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAgg                                                                         <  1:214756/150‑1 (MQ=255)
                                                                          gTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGc                                                                        >  1:129057/1‑150 (MQ=255)
                                                                             gAAGCTGGTGGGCAAGAAGGGGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGAc                                                                     >  1:746597/1‑150 (MQ=255)
                                                                             gAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                  <  2:113125/73‑1 (MQ=255)
                                                                             gAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                  >  1:113125/1‑73 (MQ=255)
                                                                                   ggtggGCAAGAAGGAGATGTACCTCCTCTAAAACCCGAACAAGGTGGTCACCGCCAGCAACGCCGCGGCCCTGCGGGCTACCCATCACCTGAACCCCGACGCTATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                               <  2:2160/150‑1 (MQ=255)
                                                                                        gCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                          <  1:482858/150‑1 (MQ=255)
                                                                                        gCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                          <  2:462962/150‑1 (MQ=255)
                                                                                           agaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGcc                                                       >  2:708918/1‑150 (MQ=255)
                                                                                               ggAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAgcg                                                   >  1:640487/1‑150 (MQ=255)
                                                                                                   aTGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGCAGCTGCACCGTGTGTCGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAt                                               <  2:249667/150‑1 (MQ=255)
                                                                                                    tGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCCTGTGACGACGTGCCCCCCATCACCTGAACCCCGTCGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtg                                              <  2:663555/150‑1 (MQ=255)
                                                                                                     gTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgt                                             >  1:186026/1‑150 (MQ=255)
                                                                                                     gTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgt                                             <  2:453036/150‑1 (MQ=255)
                                                                                                       aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGa                                           <  2:129057/150‑1 (MQ=255)
                                                                                                       aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGCCCAGCGCCATGTGa                                           >  1:442335/1‑150 (MQ=255)
                                                                                                        cATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGAt                                          <  2:234073/150‑1 (MQ=255)
                                                                                                           ccccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                       >  2:14223/1‑150 (MQ=255)
                                                                                                           ccccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCCGCCCCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                       <  1:321060/150‑1 (MQ=255)
                                                                                                             ccTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGcg                                                                                                                                      <  2:409472/53‑2 (MQ=17)
                                                                                                             ccTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGcg                                                                                                                                      >  1:409472/1‑52 (MQ=17)
                                                                                                             ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcgc                                     >  2:828437/1‑150 (MQ=255)
                                                                                                                aCAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTa                                  <  2:186026/150‑1 (MQ=255)
                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTa                                  >  1:171868/1‑136 (MQ=255)
                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTa                                  <  2:171868/136‑1 (MQ=255)
                                                                                                                                    aCCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa              <  1:708918/150‑1 (MQ=255)
                                                                                                                                       gCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGac           <  1:14223/150‑1 (MQ=255)
                                                                                                                                                gCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCt  <  1:485265/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCT  >  NC_002947/3631981‑3632274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑