Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,632,123 A→G 8.5% N325S (AAC→AGC)  PP_5590 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,632,1230AG 8.5% 128.2 / 2.9 47N325S (AAC→AGC) PP_5590hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base G (2/2);  total (24/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 = 7.26e-01

ACACCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGAC  >  NC_002947/3631978‑3632265
                                                                                                                                                 |                                                                                                                                              
acacCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGcc                                                                                                                                            >  1:2160/1‑150 (MQ=255)
   ccGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                         <  1:769887/150‑1 (MQ=255)
         ccGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCgctg                                                                                                                                   >  1:234073/1‑150 (MQ=255)
         ccGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCAACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCgctg                                                                                                                                   >  2:214756/1‑150 (MQ=255)
              cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                >  2:529360/1‑148 (MQ=255)
              cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                <  1:529360/148‑1 (MQ=255)
              cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTAAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGc                                                                                                                              >  1:249667/1‑150 (MQ=255)
                          cATGAGCTGCGATGAGCTGTACGTGTTCACCGGTCGCGTCGACCGCTTCTAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCt                                                                                                                  <  2:567325/150‑1 (MQ=255)
                                         gCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCAt                                                                                                   >  2:482858/1‑150 (MQ=255)
                                          cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGc                                                                                                                              >  2:693834/1‑122 (MQ=255)
                                          cTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGc                                                                                                                              <  1:693834/122‑1 (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)
                                                                                                                                                 |                                                                                                                                              
ACACCGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGAC  >  NC_002947/3631978‑3632265

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