Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,081,155 G→T 11.1% intergenic (+30/+85) 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,1550GT11.1% 93.5 / 3.8 36intergenic (+30/+85)efp/ohrelongation factor P/organic hydroperoxide resistance protein
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base T (2/2);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

TGCCAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCA  >  NC_002947/2081014‑2081299
                                                                                                                                             |                                                                                                                                                
tGCCAAGCTGGCGAGCGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCgg                                                                                                                                          <  1:333216/150‑1 (MQ=255)
tGCCAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCgg                                                                                                                                          <  2:528885/150‑1 (MQ=255)
    aaGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGttt                                                                                                                                      <  1:7075/150‑1 (MQ=255)
           cGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCTTGCGCCGGGTTTTTtgcgt                                                                                                                               <  2:161239/150‑4 (MQ=255)
            gAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGcc                                                                                                                              >  2:575096/1‑150 (MQ=255)
                    ccGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGCTTTTTGTGCCTGTGtttt                                                                                                                      >  2:831584/1‑150 (MQ=255)
                     cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCTTGCGCCGGGtttttt                                                                                                                                   <  1:499899/136‑1 (MQ=255)
                     cGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCTTGCGCCGGGtttttt                                                                                                                                   >  2:499899/1‑136 (MQ=255)
                           ctcGGTTGCCGACGTCATCCAGGTCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTTAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATc                                                                                                               <  1:376132/150‑1 (MQ=255)
                                    cGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGcgcg                                                                                                      >  2:44705/1‑150 (MQ=255)
                                        ttCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCt                                                                                                  >  2:584496/1‑150 (MQ=255)
                                           aTCCAGATCGACGAGTGGATCGACATCGCCACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAGTCCTCAGCGCGTGCTtgt                                                                                               <  1:627406/150‑1 (MQ=255)
                                             ccAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGc                                                                                             <  2:405267/150‑1 (MQ=255)
                                                      cgaGTGGGTCGACATCGACCCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAgg                                                                                    <  1:474275/150‑1 (MQ=255)
                                                             atcgacatcgacACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGc                                                                                                                               >  1:463765/1‑100 (MQ=255)
                                                             atcgacatcgacACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGc                                                                                                                               <  2:463765/100‑1 (MQ=255)
                                                              tcgacatcgacaCTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATcgcg                                                                            >  1:712520/1‑150 (MQ=255)
                                                                   atcgacaCTCGCGACAACAGCTTCAAAGACCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGc                                                                       >  1:238827/1‑150 (MQ=255)
                                                                            cgcgACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGcc                                                              <  1:608431/150‑1 (MQ=255)
                                                                              ggACAATAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTTCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCtt                                                            >  2:646885/2‑150 (MQ=255)
                                                                                  aacaGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTATTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGGCTTAGAc                                                        >  1:277804/1‑150 (MQ=255)
                                                                                              gaccGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAg                                            <  2:238827/148‑1 (MQ=255)
                                                                                                    cccAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGcac                                      >  2:669821/2‑150 (MQ=255)
                                                                                                       aagaagTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTg                                                                  <  2:790429/119‑1 (MQ=255)
                                                                                                       aagaagTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTg                                                                  >  1:790429/1‑119 (MQ=255)
                                                                                                        agaagTAATGCTTCTTGCGAGACGCAAAAAACCAGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCATGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATcc                                  >  2:403722/1‑150 (MQ=255)
                                                                                                                    tcttGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTATGCTCCCACAGGGGATCGCGCAACCTTTGGCGCCTTAGGCGGTTACGTGCAGGCGCACATCCACATTACCGCgt                      <  1:575096/150‑1 (MQ=255)
                                                                                                                           aggcacAAAAAACCCGGCTTGCGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCAt               >  2:722182/6‑150 (MQ=255)
                                                                                                                             aCGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCACAATCCTCAGCGCGTGCTTGTGCTCCCACAGGTGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTGCCGCGTGTGGCATTg             <  1:584496/150‑1 (MQ=255)
                                                                                                                               gCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGa           >  2:828184/1‑150 (MQ=255)
                                                                                                                                      ccccGGCGCAAGCCGTGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCCCCCATGGGGGATCGCGACCGCCTTTGCGCCTAAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACggg    <  1:831584/149‑1 (MQ=255)
                                                                                                                                      aCCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACggg    <  1:51880/150‑1 (MQ=255)
                                                                                                                                        ccGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTg                                               >  2:156455/1‑105 (MQ=255)
                                                                                                                                        ccGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTg                                               <  1:156455/105‑1 (MQ=255)
                                                                                                                                        ccGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGca  <  1:669821/150‑1 (MQ=255)
                                                                                                                                        ccGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGca  >  2:168793/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                
TGCCAAGCTGGCGAACGGTACCGAGATCTCGGTTGCCGACTTCATCCAGATCGACGAGTGGATCGACATCGACACTCGCGACAACAGCTTCAAAGGCCGCTCCAAGAAGTAATGCTTCTTGCGAGACGCAAAAAACCCGGCGCAAGCCGGGTTTTTTGTGCCTGTGTTTTCCCAATCCTCAGCGCGTGCTTGTGCTCCCACAGGGGATCGCGCAAGCCTTTGCGCCTTAGACGGTTACGTGCAGGCGCACATCCACATTACCGCGTGTGGCATTGGAGTACGGGCA  >  NC_002947/2081014‑2081299

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: