Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,081,156 C→T 20.5% intergenic (+31/+84) efp → / ← ohr elongation factor P/organic hydroperoxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,081,1560CT20.5% 38.8 / 4.5 24intergenic (+31/+84)efp/ohrelongation factor P/organic hydroperoxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base T (3/2);  total (14/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.99e-01

CAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACG  >  NC_002947/2081017‑2081295
                                                                                                                                           |                                                                                                                                           
cAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGtt                                                                                                                                   >  1:505395/1‑150 (MQ=255)
       gCGAACGGTACTGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGTGTTTTTtgtg                                                                                                                            >  2:379543/1‑150 (MQ=255)
       gCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTtgtg                                                                                                                            >  1:265444/1‑150 (MQ=255)
                  cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTc                                                                                                                 >  1:198019/1‑150 (MQ=255)
                                          ccAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGc                                                                                         >  2:188581/1‑150 (MQ=255)
                                                        ggATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATc                                                                           <  2:505395/150‑1 (MQ=255)
                                                              acatcgacaCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCaa                                                                     >  2:118367/1‑150 (MQ=255)
                                                                   gacaCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACACAAAAAACCCGGCTTGCGCCGGGttttt                                                                                                                                >  1:444399/1‑86 (MQ=37)
                                                                        tCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTgcgc                                                           <  2:265444/150‑1 (MQ=255)
                                                                                       cAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGTGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGt                                            <  1:379543/150‑1 (MQ=255)
                                                                                          aGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCa                                         >  1:175826/1‑150 (MQ=255)
                                                                                                 tCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGcac                                  >  1:267594/1‑150 (MQ=255)
                                                                                                   cAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACAt                                >  2:428777/1‑150 (MQ=255)
                                                                                                                     gCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTgg              <  1:428777/150‑1 (MQ=255)
                                                                                                                       gagaCGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACcg                     <  1:64194/141‑1 (MQ=255)
                                                                                                                       gagaCGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACcg                     >  2:64194/1‑141 (MQ=255)
                                                                                                                          gcacAAAAAACCCGGCTTGCGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGcac                                  >  1:307279/4‑125 (MQ=39)
                                                                                                                          gcacAAAAAACCCGGCTTGCGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGcac                                  <  2:307279/122‑1 (MQ=38)
                                                                                                                          gcacAAAAAACCCGGCTTGCGCCGGGTTTTTTGCGTCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGtgtg               <  1:25764/141‑1 (MQ=38)
                                                                                                                          gcacAAAAAACCCGGCTTGCGCCGGGTTTTTTGCGTCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGtgtg               >  2:25764/4‑144 (MQ=38)
                                                                                                                            gCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGCAGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa       <  1:118367/150‑1 (MQ=255)
                                                                                                                            gCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa       >  1:202267/1‑150 (MQ=255)
                                                                                                                            gCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa       <  2:267594/150‑1 (MQ=255)
                                                                                                                                 aaaCCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGAACATCCACATTACCGCGTGTGGCATTGGAGTACg  <  1:188581/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACG  >  NC_002947/2081017‑2081295

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: