Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 747,612:1 +G 100% intergenic (‑739/‑55) PP_0638 ← / → PP_5448 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947747,6121.G92.3% 35.7 / ‑2.6 13intergenic (‑739/‑55)PP_0638/PP_5448transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (12/0);  total (13/0)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ATCATCGGTTAGCGTGGCCACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATC  >  NC_002947/747468‑747749
                                                                                                                                                 |                                                                                                                                          
atcatcGGTTAGCGTGGCCACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTgggg                                                                                                                                        <  4:280898/149‑1 (MQ=255)
   atcGGTTAGCGTGGCCACTCACCATGTAGGAGGATCTCCGTATGAGATCGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGTCTCTCACTGGCCGACACCACTCTGGGTCATCTTCGTGTTGTGTGTTGTgggggg                                                                                                                                     <  6:326543/148‑1 (MQ=255)
        ttAGCGTGGCCACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGTGGGgcgc                                                                                                                                >  3:185918/1‑149 (MQ=255)
                 ccACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTgggggggg                                                                                                                                   <  3:250501/136‑1 (MQ=255)
                 ccACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTgggggggg                                                                                                                                   >  4:250501/1‑136 (MQ=255)
                                   gaTCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGc                                                                                                     >  1:930/1‑149 (MQ=255)
                                                          gttgtGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTAGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATgc                                                                              >  3:274163/1‑149 (MQ=255)
                                                           ttgtGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGGTGAATGGGTGGAGCGGAGGGGGCGCAACCAATGATgcg                                                                             >  3:184318/1‑149 (MQ=255)
                                                                  ttGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCaa                                                                       >  7:22519/1‑148 (MQ=255)
                                                                    gggCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGGTGAATGGCTGGGGCGGTGGGGGCGCAATCAATGTTGCGCGCCAAtg                                                                     >  5:195409/1‑147 (MQ=255)
                                                                    gggCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCg                                                                                                    >  7:37584/1‑117 (MQ=255)
                                                                                    gCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATc                                                    >  3:40646/1‑149 (MQ=255)
                                                                                      ttCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGAc                                                                >  3:231977/1‑135 (MQ=255)
                                                                                      ttCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGt                                                  >  1:253472/1‑149 (MQ=255)
                                                                                                                          gTCATCGTCGTGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTTGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGTTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCaatatc              >  1:226381/1‑144 (MQ=255)
                                                                                                                                 tcgtGTTGTGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTgcg       >  3:8043/1‑149 (MQ=255)
                                                                                                                                      ttgtGGGTGGTGGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCTATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATc  >  5:233746/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
ATCATCGGTTAGCGTGGCCACTCACCATGTAGGAGGATCTCCGTATGAGAACGGTTTTGTTGTGTCTTGGGCTTGCAAGTGCAGGCTTCTGTGTGGGCTCTCACTGGCCGACACCACTCTGGGTCATCGTCGTGTTGTGGGTGGTGGGGGGGGCGCTGTGGGTTCTGGCTGAATGGCTGGAGCGGAGGGGGCGCAATCAATGATGCGCGCCAATCTAGACGAGCTGGTGATCGTGAGCTGTTTGTTTCCAGCTATGAGGTGGTCAGCTTCTATTGCGCGATC  >  NC_002947/747468‑747749

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