Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,999,934 A→G 48.7% V6V (GTT→GTC yecC ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,999,9340AG48.7% ‑6.6 / 54.5 37V6V (GTT→GTCyecCpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base G (6/12);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GCTTTACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCT  >  W3110S.gb/1999870‑1999997
                                                                |                                                               
gCTTTACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACtt                                                             <  1:1005508/69‑1 (MQ=255)
gCTTTACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACtt                                                             <  1:2815668/69‑1 (MQ=255)
  tttACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCa                                                           >  1:1156995/1‑69 (MQ=255)
  tttACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCa                                                           >  1:388187/1‑69 (MQ=255)
  tttACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCa                                                           >  1:2820032/1‑69 (MQ=255)
      ccTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCa                                                           >  1:495178/1‑65 (MQ=255)
      ccTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCa                                                           <  1:1420786/65‑1 (MQ=255)
         cAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCa                                                     >  1:1545168/1‑68 (MQ=255)
                  tACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTgg                                           <  1:2620597/69‑1 (MQ=255)
                  tACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTgg                                           <  1:2616150/69‑1 (MQ=255)
                  tACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTgg                                           <  1:2559196/69‑1 (MQ=255)
                  tACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCAATGGCACTCATTTTgg                                           <  1:92550/69‑1 (MQ=255)
                  tACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCAATGGCACTCATTTTgg                                           <  1:1223620/69‑1 (MQ=255)
                   aCCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTgg                                           >  1:3312910/1‑68 (MQ=255)
                   aCCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTGGt                                          <  1:1008516/69‑1 (MQ=255)
                    ccGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTGGt                                          <  1:2848911/68‑1 (MQ=255)
                                 tCTGACCGTGGAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGc                             <  1:3078938/68‑1 (MQ=255)
                                         ggaaaTTTTTTCACCAGGTTCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTa                        <  1:2537854/62‑1 (MQ=255)
                                           gAATTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGa                  <  1:1356897/69‑1 (MQ=255)
                                            aaTTTTTTCACCAGGTTCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGa                  <  1:700819/68‑1 (MQ=255)
                                                     aCCAGGTTCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCaaaa         <  1:1090787/68‑1 (MQ=255)
                                                      ccAGGTTCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAAt        <  1:1708738/68‑1 (MQ=255)
                                                       cAGGTTCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATgg      <  1:1319288/69‑1 (MQ=255)
                                                            tCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGa                  <  1:3281444/52‑1 (MQ=255)
                                                            tCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:88935/1‑68 (MQ=255)
                                                            tCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:2028927/1‑68 (MQ=255)
                                                            tCTTGACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:1670959/1‑68 (MQ=255)
                                                            tCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCaa           >  1:3228397/1‑59 (MQ=255)
                                                            tCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCaa           >  1:3583699/1‑59 (MQ=255)
                                                            tCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:1622790/1‑68 (MQ=255)
                                                            tCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:1115516/1‑68 (MQ=255)
                                                            tCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCt  >  1:3707544/1‑68 (MQ=255)
                                                             cTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGtt                    >  1:3161009/1‑49 (MQ=255)
                                                             cTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGtt                    >  1:2698746/1‑49 (MQ=255)
                                                             cTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGtt                    <  1:2692516/49‑1 (MQ=255)
                                                             cTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGtt                    >  1:1571204/1‑49 (MQ=255)
                                                             cTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGa                  >  1:1983260/1‑51 (MQ=255)
                                                                |                                                               
GCTTTACCTCAAGGTCGATACCGTGCAGCACCGTCTGACCGTGGAATTTTTTCACCAGGTTCTTAACTTCAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCT  >  W3110S.gb/1999870‑1999997

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: