Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,104 C→T 19.0% intergenic (‑34/+19) PP_5445 ← / ← PP_t71 hypothetical protein/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,1040CT19.0% 38.0 / 3.4 21intergenic (‑34/+19)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base C (9/8);  new base T (2/2);  total (11/10)
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.34e-01

GCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGC  >  NC_002947/466957‑467241
                                                                                                                                                   |                                                                                                                                         
gcgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgc                                                                                                                                         >  1:388039/1‑150 (MQ=255)
  gcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTAGCTGCGATTCCGATTGGGGGCTTCTGCCTAGGTCATGCGCCAGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgc                                                                                                                                       <  1:291370/150‑1 (MQ=255)
           ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                              >  2:187674/1‑150 (MQ=255)
            cGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCg                                                                                                                             >  1:36573/1‑150 (MQ=255)
                         gggCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                <  2:36573/150‑1 (MQ=255)
                          ggCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGccc                                                                                                               <  2:388039/150‑1 (MQ=255)
                           gCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCca                                                                                                              >  2:474538/1‑150 (MQ=255)
                               aTTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACg                                                                                                          <  1:474538/150‑1 (MQ=255)
                                              cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                           >  2:282263/1‑150 (MQ=255)
                                                   tAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                      <  1:282263/150‑1 (MQ=255)
                                                   tAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGAAATCGCGCGGCTTCGTCGTAACTGGTGGGCCCAGACGGATCCGAACCGTGGACCaaa                                                                                      <  2:179641/150‑1 (MQ=255)
                                                            ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                             <  1:187674/150‑1 (MQ=255)
                                                                                           tgGATGTACGCGCCCCTGCCCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCCTGGACCAAAAGATGATGAGGCCTCTGCTCTAACCAACTGAGCTATAgcc                                              >  2:586304/1‑148 (MQ=255)
                                                                                            ggATGTACGCGCTCCTGACCATTTATTTCTTGCGGGCATGAAACAGCCGCGGAATCGTGCGGCGTTGTCGTATCTGGTAGGCCCACACGGACTAGAACCGTGGGACAAAGGATTATGAGTCGTCTGCTCTAAACAACTGAACTATAGGcc                                             >  2:31636/1‑150 (MQ=255)
                                                                                            ggATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGcc                                             <  1:590348/150‑1 (MQ=255)
                                                                                                                         ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAggg                >  2:150653/1‑150 (MQ=255)
                                                                                                                                  tGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCAcc       >  1:52980/1‑150 (MQ=255)
                                                                                                                                      aaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                              >  2:221398/1‑91 (MQ=255)
                                                                                                                                      aaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                              <  1:221398/91‑1 (MQ=255)
                                                                                                                                      aaaGCCGCGCGATTGCGCGGCTTTGACGTAGCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgtg   >  1:570687/1‑150 (MQ=255)
                                                                                                                                       aaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGc  <  1:233430/150‑1 (MQ=255)
                                                                                                                                            gcgcGATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                      <  2:353207/61‑1 (MQ=255)
                                                                                                                                            gcgcGATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                      >  1:353207/1‑61 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGC  >  NC_002947/466957‑467241

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