Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 919,675 G→C 50.0% intergenic (‑143/+48) PP_5458 ← / ← PP_0800 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947919,6750GC50.0% ‑6.2 / 11.8 10intergenic (‑143/+48)PP_5458/PP_0800hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base C (2/3);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAGGT  >  NC_002947/919537‑919824
                                                                                                                                          |                                                                                                                                                      
aaCTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGCCGAAAGCCCCg                                                                                                                                             <  1:175054/150‑1 (MQ=255)
      cttcttGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGCCGAAAGcccc                                                                                                                                              >  1:501509/1‑143 (MQ=255)
      cttcttGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGCCGAAAGcccc                                                                                                                                              <  2:501509/143‑1 (MQ=255)
                     gTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAggg                                                                                                                        >  1:359719/1‑150 (MQ=255)
                                gACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTgg                                                                                                             <  2:495706/150‑1 (MQ=255)
                                                               tCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAAc                                                                              >  1:240417/1‑150 (MQ=255)
                                                                                             ttGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTTGCGCACAACCCat                                               <  2:359719/150‑2 (MQ=255)
                                                                                                   aaaaGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCtc                                         >  2:198929/1‑149 (MQ=255)
                                                                                                                                 gggCTTTCGCCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGc                                            >  1:41903/1‑117 (MQ=255)
                                                                                                                                 gggCTTTCGCCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGc                                            <  2:41903/117‑1 (MQ=255)
                                                                                                                                          gCGAAAGCCGCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAGGt  >  2:430296/1‑150 (MQ=255)
                                                                                                                                          gCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAGGt  >  1:261554/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                      
AACTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAGGT  >  NC_002947/919537‑919824

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: