Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,186,703 T→C F216F (TTT→TTC yqiG → predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,186,7030TC76.7% 116.9 / 33.5 61F216F (TTT→TTCyqiGpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base C (21/25);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAC  >  W3110S.gb/3186638‑3186766
                                                                 |                                                               
cGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCcgcg                                                            <  1:400561/71‑1 (MQ=255)
    cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:2868943/1‑71 (MQ=255)
    cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:1068512/1‑71 (MQ=255)
     tGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCcgcg                                                            <  1:1863192/66‑1 (MQ=255)
     tGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCcgcg                                                            <  1:1679566/66‑1 (MQ=255)
     tGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCcgcg                                                            <  1:916947/66‑1 (MQ=255)
     tGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCcgcg                                                            <  1:1263814/66‑1 (MQ=255)
       aTTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGtattta                                                      <  1:3457350/70‑1 (MQ=255)
       aTTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGta                                                          >  1:1797393/1‑66 (MQ=255)
        ttACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGTCCAACCTGGATTTCCCGCGtattta                                                      <  1:2462304/69‑1 (MQ=255)
         tACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTatttattt                                                   >  1:262700/1‑71 (MQ=255)
          aCCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGtatttatt                                                    >  1:768005/1‑69 (MQ=255)
          aCCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGtatttatt                                                    >  1:2292120/1‑69 (MQ=255)
          aCCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGtatttatt                                                    >  1:3301384/1‑69 (MQ=255)
          aCCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGtatttatt                                                    >  1:1195366/1‑69 (MQ=255)
           ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTg                                                  <  1:3971743/70‑1 (MQ=255)
           ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGt                                                 <  1:435232/71‑1 (MQ=255)
           ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGt                                                 <  1:118076/71‑1 (MQ=255)
           ccAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGt                                                 <  1:3287714/71‑1 (MQ=255)
            cAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTg                                                  <  1:475698/69‑1 (MQ=255)
              gTACGATCAGTCATTTTACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTg                                                  <  1:3319342/67‑1 (MQ=255)
                 cGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGcc                                           >  1:444731/1‑71 (MQ=255)
                 cGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGcc                                           >  1:2104628/1‑71 (MQ=255)
                 cGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGcc                                           >  1:1219443/1‑71 (MQ=255)
                    tCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCaa                                        >  1:225768/1‑71 (MQ=255)
                    tCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCaa                                        >  1:938340/1‑71 (MQ=255)
                        tCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGccc                                          >  1:709888/1‑65 (MQ=255)
                        tCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGccc                                          >  1:1430498/1‑65 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:267031/1‑49 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:1006625/1‑49 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:313412/1‑49 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:3484828/1‑49 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:3228019/1‑49 (MQ=255)
                          aTTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                        >  1:713633/1‑49 (MQ=255)
                                      gaaaaGCCAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGt                                                 <  1:491758/43‑1 (MQ=39)
                                             cAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCATCAATGAATGCAAaactaac               >  1:1544397/1‑71 (MQ=255)
                                             cAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaac               >  1:3712872/1‑71 (MQ=255)
                                             cAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaac               <  1:1021861/71‑1 (MQ=255)
                                             cAGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaac               >  1:1755735/1‑71 (MQ=255)
                                             cAGGCGACCAACCTGGATTTACCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaac               >  1:3327853/1‑71 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:1094937/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:1430027/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:64203/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:561972/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:1489244/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:2877944/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:314864/68‑1 (MQ=255)
                                              aGGCGACCAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                 <  1:3482452/68‑1 (MQ=255)
                                                     cAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt        >  1:2258079/1‑70 (MQ=255)
                                                     cAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt        >  1:428443/1‑70 (MQ=255)
                                                     cAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt        >  1:1528742/1‑70 (MQ=255)
                                                     cAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt        >  1:344368/1‑70 (MQ=255)
                                                     cAACCTGGATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGt        >  1:642024/1‑70 (MQ=255)
                                                         cTGGATTTCCCGCGTATTTATTTGTTTCGCGCAATCCCAGCAATGAATGCAAaacta                 <  1:1440961/57‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCACTAc  <  1:1209581/69‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:1785409/69‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:3364687/69‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:3360002/69‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:1219324/69‑1 (MQ=255)
                                                            gATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:3343112/69‑1 (MQ=255)
                                                            gATTTCCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc  <  1:1577497/69‑1 (MQ=255)
                                                                 |                                                               
CGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAC  >  W3110S.gb/3186638‑3186766

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: