Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,796,665 C→T 41.2% A68T (GCG→ACG)  btuC ← vitamin B12 transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,796,6650CT41.2% 5.3 / 16.1 17A68T (GCG→ACG) btuCvitamin B12 transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base T (3/4);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GTCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACGAG  >  W3110S.gb/1796599‑1796730
                                                                  |                                                                 
gtCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCacc                                                                >  1:348349/1‑70 (MQ=255)
      gTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGc                                                         <  1:1809904/71‑1 (MQ=255)
        tCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCaa                                                       >  1:3097120/1‑71 (MQ=255)
         cTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCaa                                                       <  1:287013/70‑1 (MQ=255)
            ccAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAc                                                    <  1:2978916/70‑1 (MQ=255)
            ccAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAc                                                    <  1:3083401/70‑1 (MQ=255)
                   aTTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAcagccagc                                             <  1:374482/70‑1 (MQ=255)
                      ttCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGt                                           >  1:1738450/1‑69 (MQ=255)
                      ttCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTAc                                         >  1:48659/1‑71 (MQ=255)
                      ttCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAca                                                   >  1:2253329/1‑61 (MQ=255)
                         aaaCAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTg                                      <  1:1496419/71‑1 (MQ=255)
                         aaaCAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTg                                      >  1:3671327/1‑71 (MQ=255)
                                  cTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGaa                             <  1:1243587/71‑1 (MQ=255)
                                                    tataGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGtt           >  1:791140/1‑71 (MQ=255)
                                                    tataGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGtt           >  1:323240/1‑71 (MQ=255)
                                                        gCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGcc       <  1:2523577/71‑1 (MQ=255)
                                                          cAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCAc     <  1:127564/71‑1 (MQ=255)
                                                             cgccgcACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACgag  <  1:2441468/71‑1 (MQ=255)
                                                             cgccgcACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACgag  <  1:1484651/71‑1 (MQ=255)
                                                                  |                                                                 
GTCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACGAG  >  W3110S.gb/1796599‑1796730

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