Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,298,568 T→C 32.0% E49E (GAA→GAG ccmE ← periplasmic heme chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,298,5680TC32.0% 25.1 / 19.5 25E49E (GAA→GAGccmEperiplasmic heme chaperone
Reads supporting (aligned to +/- strand):  ref base T (11/6);  new base C (6/2);  total (17/8)
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.63e-01

CACTACCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAGCGCAGCGCA  >  W3110S.gb/2298504‑2298631
                                                                |                                                               
cacTACCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGc                                                           >  1:1225205/1‑71 (MQ=255)
  cTACCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCtt                                                         <  1:1467238/71‑1 (MQ=255)
  cTACCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCtt                                                         <  1:1470520/71‑1 (MQ=255)
    aCCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGc                                                       >  1:2600231/1‑71 (MQ=255)
        ggCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCg                                                     >  1:2145315/1‑69 (MQ=255)
           aTCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTa                                                   >  1:448326/1‑68 (MQ=255)
              aCCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGaa                                              >  1:2600469/1‑70 (MQ=255)
              aCCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAAt                                             >  1:641844/1‑71 (MQ=255)
              aCCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAAt                                             >  1:2597168/1‑71 (MQ=255)
                 aTCCCGCCAACGCGCAGCCGCTGACCGACTTCCGGCATTTGCTGAGTTTCa                                                              >  1:233336/1‑51 (MQ=255)
                    ccGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTccc                                        >  1:2453525/1‑70 (MQ=255)
                       ccAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGaa                                              <  1:2652158/61‑1 (MQ=255)
                         aaCGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCg                                  >  1:2370507/1‑71 (MQ=255)
                                           gACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGatat                >  1:640159/1‑71 (MQ=255)
                                                ccGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCg             <  1:201283/69‑1 (MQ=255)
                                                ccGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCg             <  1:2600412/69‑1 (MQ=255)
                                                 cGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGag           <  1:2734008/70‑1 (MQ=255)
                                                        tGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    >  1:2470509/1‑70 (MQ=255)
                                                        tGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    >  1:2111564/1‑70 (MQ=255)
                                                        tGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    >  1:2617830/1‑70 (MQ=255)
                                                        tGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    >  1:2707315/1‑70 (MQ=255)
                                                        tGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    <  1:323543/70‑1 (MQ=255)
                                                        tGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcagcg    >  1:1037747/1‑70 (MQ=255)
                                                          cTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAgcgcagcgca  <  1:2226251/70‑1 (MQ=255)
                                                          cTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAgcgcagcgca  >  1:1390722/1‑70 (MQ=255)
                                                                |                                                               
CACTACCCGGCATCACCATCCCGCCAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAGCGCAGCGCA  >  W3110S.gb/2298504‑2298631

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