Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,315,094 C→A 13.8% A94A (GCC→GCA)  PP_1146 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,315,0940CA13.8% 73.8 / 2.8 29A94A (GCC→GCA) PP_1146hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/16);  new base A (2/2);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GCACTTTCGTGACAGTGACGGCGAAGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAATGGAGTGGCCGAGATAACGGTGGCCAGCCA  >  NC_002947/1314965‑1315240
                                                                                                                                 |                                                                                                                                                  
gCACTTTCGTGACAGTGACGGCGAAGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGgc                                                                                                                                >  2:129299/1‑150 (MQ=255)
                 aCGGCGAAGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCgaga                                                                                                                                                >  2:322039/1‑117 (MQ=255)
                 aCGGCGAAGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCgaga                                                                                                                                                <  1:322039/117‑1 (MQ=255)
                       aaGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGaaa                                                                                                         >  1:127200/1‑150 (MQ=255)
                               ctgGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGc                                                                                                 <  1:336672/150‑1 (MQ=255)
                                   tCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGcc                                                                                                                                        >  2:424104/1‑107 (MQ=255)
                                   tCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGcc                                                                                                                                        <  1:424104/107‑1 (MQ=255)
                                     cTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTAAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCgggg                                                                                           <  1:319201/150‑1 (MQ=255)
                                       gAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCa                                                                                                              >  1:28921/1‑129 (MQ=255)
                                       gAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCa                                                                                                              <  2:28921/129‑1 (MQ=255)
                                       gAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCAGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCt                                                                                         >  2:468046/1‑150 (MQ=255)
                                           cGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGAc                                                                                     <  2:342435/150‑1 (MQ=255)
                                             cAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCt                                                                                   >  2:293052/1‑150 (MQ=255)
                                                    ggccaggccaGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGAtgtg                                                                            <  1:44212/150‑1 (MQ=255)
                                                    ggccaggccaGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTGGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGAtgtg                                                                            <  2:227017/150‑1 (MQ=255)
                                                                   cTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAg                                                                                                                                  >  2:47258/1‑81 (MQ=255)
                                                                   cTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAg                                                                                                                                  <  1:47258/81‑1 (MQ=255)
                                                                   cTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGAc                                                             <  1:197402/150‑1 (MQ=255)
                                                                        aaGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCAGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATaa                                                                 >  2:166659/1‑141 (MQ=255)
                                                                        aaGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCAGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATaa                                                                 <  1:166659/141‑1 (MQ=255)
                                                                                aGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCa                                                <  1:126477/150‑1 (MQ=255)
                                                                                     gTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGAc                                           <  1:274340/150‑1 (MQ=255)
                                                                                        gACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGGCTGGAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTg                                        >  2:299834/1‑150 (MQ=255)
                                                                                                    cTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAATGGa                            <  1:129299/150‑1 (MQ=255)
                                                                                                        aGGCCACCTTGTACGGGCGCACGGCAGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAATGGAGTgg                        <  1:468046/150‑1 (MQ=255)
                                                                                                          gCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTg                                                                                  <  2:145371/90‑1 (MQ=255)
                                                                                                          gCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTg                                                                                  >  1:145371/1‑90 (MQ=255)
                                                                                                            cACCCTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAACGGCGTGGcgtc                    <  2:127200/150‑4 (MQ=255)
                                                                                                                              ggCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAATGGAGTGGCCGAGATAACGGTGgccagcca  <  1:293052/150‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                  
GCACTTTCGTGACAGTGACGGCGAAGGGCTGCTGGTCCTGAGCCGCAGCGATGGCCAGGCCAGCGACCTGGAAAGCGAACAGGAAGTGGACAAGGTGGTGCTCAAGGCCACCTTGTACGGGCGCACGGCCGAGAGCGACGCCTTCAAGGCGCGTTGGCAGGTCGAGCAGGAAACGACCTGCCCGGGGCTGGACCTGGATGTGGACTTCTATAATGACGTCAGTGATGTCAGCGACCTGAACCAGAATGGAGTGGCCGAGATAACGGTGGCCAGCCA  >  NC_002947/1314965‑1315240

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