Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,363,041 T→C 15.4% L127S (TTA→TCA)  yhcD → predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,363,0410TC15.4% 94.8 / 11.0 39L127S (TTA→TCA) yhcDpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (3/3);  total (21/18)
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.20e-01

GTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACA  >  W3110S.gb/3362978‑3363093
                                                               |                                                    
gTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGatat                                                >  1:457077/1‑70 (MQ=255)
gTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAAACAGGAAGCGTTAGata                                                 >  1:386307/1‑69 (MQ=255)
      aCTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGa                                                   >  1:1678407/1‑61 (MQ=255)
       cTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTa                                                     <  1:652063/58‑1 (MQ=255)
        tatataTTGTCCTACCCGGGTACCACAACCACACCATTACCTAACCAGGAAGCGTTAGata                                                 >  1:1016847/1‑61 (MQ=255)
              ttGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTcacca                                   >  1:517960/1‑69 (MQ=255)
                      cccGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTCAGATATTATCGTTTCACCACAGGCGAt                           <  1:955080/69‑1 (MQ=255)
                       ccGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAt                        <  1:1440399/71‑1 (MQ=255)
                        cGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTCAGAtattat                                             >  1:689168/1‑49 (MQ=255)
                            tACCACAATCACACCATTACCTGACCAGGAAGCGTTAGATATTATCGTTTcaccac                                  <  1:832208/56‑1 (MQ=255)
                            tACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGt                                          <  1:1377067/48‑1 (MQ=255)
                             aCCACAATCACACCATTACCTAACCAGGAAGCTTTAGATATTATCGTTTcac                                     <  1:1744906/52‑1 (MQ=255)
                             aCCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTcac                                     <  1:1266378/52‑1 (MQ=255)
                             aCCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTcac                                     <  1:1447275/52‑1 (MQ=255)
                              ccACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTcacc                                    >  1:1047066/1‑52 (MQ=255)
                                acaATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGtt                                         >  1:1542252/1‑45 (MQ=255)
                                acaATCACACCATTACCTAACCAGGAAGCGTCAGATATTATCGTTTCACCACAGGCGATCATTCCCATAgg               >  1:1101037/1‑71 (MQ=255)
                                  aaTCACACCATTACCTAACCAGGAAGCGTCAGATATTATCGTTTCACCACAGGCGATCATTCCCATAg                <  1:1070812/68‑1 (MQ=255)
                                    tCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTg           >  1:779022/1‑71 (MQ=255)
                                       cacCATTACCTAACCAGGAAGCGTTAGATATTATCGt                                          >  1:37640/1‑37 (MQ=255)
                                       cacCATTACCTAACCAGGAAGCGTTAGATATTATCGt                                          >  1:447286/1‑37 (MQ=255)
                                       cacCATTACCTAACCAGGAAGCGTTAGATATTATCGt                                          >  1:428049/1‑37 (MQ=255)
                                       cacCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAt        <  1:1051186/71‑1 (MQ=255)
                                        acCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtc       >  1:959914/1‑71 (MQ=255)
                                         ccATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct      <  1:200915/71‑1 (MQ=255)
                                         ccATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct      <  1:424534/71‑1 (MQ=255)
                                         ccATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct      <  1:1592990/71‑1 (MQ=255)
                                          cATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCa                   <  1:896299/57‑1 (MQ=255)
                                             tACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAtt                       >  1:384712/1‑50 (MQ=255)
                                             tACCTAACCAGGAAGCGTTAGATATTATCGTTTCACAACAGGc                              >  1:777181/1‑43 (MQ=255)
                                              aCCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCa                         <  1:506807/47‑1 (MQ=255)
                                              aCCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCa                         <  1:1636220/47‑1 (MQ=255)
                                              aCCTAACCAGGAAGCGTCAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcaca  >  1:977263/1‑70 (MQ=255)
                                                cTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAtt                       >  1:1388816/1‑47 (MQ=255)
                                                cTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcaca  <  1:467509/68‑1 (MQ=255)
                                                cTAACCAGGAAGCGTCAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcaca  <  1:1684913/68‑1 (MQ=255)
                                                  aaCCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAt        >  1:1185653/1‑60 (MQ=255)
                                                  aaCCAGGAAGCGTTAGATATTATCGTTGCACCACAGGCGATCATTCCCATAGGGTTGGAt        >  1:1694853/1‑60 (MQ=255)
                                                      aGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTca    >  1:1263268/1‑60 (MQ=255)
                                                               |                                                    
GTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACA  >  W3110S.gb/3362978‑3363093

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: