Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,796,665 C→T 27.3% 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,6650CT27.3% 53.9 / 18.7 33A68T (GCG→ACG) btuCvitamin B12 transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (7/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACGAG  >  W3110S.gb/1796599‑1796730
                                                                  |                                                                 
gtCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCacc                                                                >  1:4361151/1‑70 (MQ=255)
     ggTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCaccaac                                                             >  1:1371057/1‑68 (MQ=255)
      gTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGc                                                         <  1:1892625/71‑1 (MQ=255)
        tCTGCCATAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCaa                                                       >  1:630989/1‑71 (MQ=255)
          tGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCaa                                                       <  1:2731044/69‑1 (MQ=255)
            ccAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATa                                                     >  1:1143835/1‑69 (MQ=255)
                agGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATAcagc                                                 >  1:3829454/1‑69 (MQ=255)
                  gATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAcagccag                                              <  1:3542069/70‑1 (MQ=255)
                  gATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATAcagcc                                                >  1:3303258/1‑68 (MQ=255)
                  gATTTTAAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCa                                                                  >  1:1157004/1‑50 (MQ=255)
                     tttCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAg                                                          <  1:4510309/55‑1 (MQ=255)
                        cAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGt                                       >  1:480804/1‑71 (MQ=255)
                            caacGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCa                                   <  1:413018/71‑1 (MQ=255)
                             aacGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAAt                                                      <  1:4263194/51‑1 (MQ=255)
                                   tGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAAt                                                      <  1:4635761/45‑1 (MQ=255)
                                     cATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGc                                    <  1:4230998/61‑1 (MQ=255)
                                     cATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGc                                    <  1:4532879/61‑1 (MQ=255)
                                     cATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAAttt                          >  1:4037133/1‑71 (MQ=255)
                                      aTTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTg                         <  1:234762/71‑1 (MQ=255)
                                       ttACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCaa                                                       <  1:184977/40‑1 (MQ=255)
                                         aCAGCGCCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTAc                                         >  1:3027331/1‑52 (MQ=255)
                                         aCAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAAt                            >  1:3722317/1‑65 (MQ=255)
                                         aCAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGc                        >  1:2767110/1‑69 (MQ=255)
                                         aCAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGc                        >  1:3339057/1‑69 (MQ=255)
                                          cAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGc                        >  1:2488466/1‑68 (MQ=255)
                                          cAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGc                        >  1:2493906/1‑68 (MQ=255)
                                              gcCGGATATAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACg                  <  1:4242566/70‑1 (MQ=255)
                                                   atatAGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTa                                          >  1:1341064/1‑41 (MQ=255)
                                                   atatAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGaa                <  1:1939828/67‑1 (MQ=255)
                                                    tataGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGtt           >  1:3977366/1‑71 (MQ=255)
                                                     ataGCCAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTc          >  1:1979565/1‑71 (MQ=255)
                                                         ccAGCGCCGTACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGAc                   >  1:1607407/1‑58 (MQ=255)
                                                            gcgcCGCACCAACCAGCAATACAGCCAGCGTAAGTGGCAGGCGaa                             <  1:3751498/45‑1 (MQ=255)
                                                             cgccgtACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACgag  <  1:2405828/71‑1 (MQ=255)
                                                                cgcACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGt            >  1:3402438/1‑58 (MQ=255)
                                                                cgcACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGt            >  1:607599/1‑58 (MQ=255)
                                                                  |                                                                 
GTCCAGGTTCTGCCAGAGGATTTTCAAACAACGCCTGCATTACAGCGCCGGATATAGCCAGCGCCGCACCAACCAGCAATACAGCCAGCGTACGTGGCAGGCGAATTTGCCAGACGAACAGTTCGCCACGAG  >  W3110S.gb/1796599‑1796730

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