Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 219,687 A→T 14.3% A103A (GCT→GCA rcsF ← predicted outer membrane protein, signal

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb219,6870AT14.3% 104.7 / 11.0 42A103A (GCT→GCArcsFpredicted outer membrane protein, signal
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base T (4/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

CCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAG  >  W3110S.gb/219619‑219756
                                                                    |                                                                     
ccGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAg                                                                      >  1:687570/1‑70 (MQ=255)
    tacatacaGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTg                                                                 <  1:1525268/71‑1 (MQ=255)
     acatacaGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTgg                                                                >  1:124040/1‑71 (MQ=255)
      catacaGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGc                                                               <  1:1575598/71‑1 (MQ=255)
       atacaGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGc                                                                     <  1:1790827/64‑1 (MQ=255)
          caGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCttt                                                            >  1:1797206/1‑70 (MQ=255)
                  cGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTaga                                                   <  1:1047962/71‑1 (MQ=255)
                  cGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTa                                                     >  1:1730168/1‑69 (MQ=255)
                   gATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTagag                                                  >  1:955123/1‑71 (MQ=255)
                           ccTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAt                                          <  1:231563/71‑1 (MQ=255)
                            cTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAtt                                         <  1:280759/71‑1 (MQ=255)
                             tGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTc                                                           >  1:2091035/1‑52 (MQ=255)
                             tGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAttt                                        >  1:847256/1‑71 (MQ=255)
                               gCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTg                                       >  1:420453/1‑70 (MQ=255)
                                  tACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTaga                                                   <  1:194655/55‑1 (MQ=255)
                                   aCCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGtt                                             <  1:339947/60‑1 (MQ=255)
                                   aCCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATcc                                  >  1:1957774/1‑71 (MQ=255)
                                       cTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCtt                              <  1:768973/71‑1 (MQ=255)
                                        tGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGc                                                               <  1:1203444/37‑1 (MQ=255)
                                             aCTTCGCTGCTATGCAGTAATACAGCATTGGCTTTCAt                                                         >  1:1651401/1‑38 (MQ=255)
                                              cTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCg                       <  1:180920/71‑1 (MQ=255)
                                                 cGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGt                     <  1:2010876/70‑1 (MQ=255)
                                                  gCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTg                   <  1:975730/71‑1 (MQ=255)
                                                   cAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTgg                  <  1:83197/71‑1 (MQ=255)
                                                   cAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTgg                  <  1:1012874/71‑1 (MQ=255)
                                                    aGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTgg                  >  1:1309891/1‑70 (MQ=255)
                                                     gCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAttt                                        <  1:1036433/47‑1 (MQ=255)
                                                      caaTGCTGTATTACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGaa                <  1:1627525/68‑1 (MQ=255)
                                                      caaTGCTGTATTACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGaa                <  1:1164331/68‑1 (MQ=255)
                                                       tATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCt                               <  1:1660435/54‑1 (MQ=255)
                                                         tGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc             <  1:96515/70‑1 (MQ=255)
                                                            aGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAt                                          <  1:1279883/38‑1 (MQ=255)
                                                            aGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAAt               >  1:1410412/1‑65 (MQ=255)
                                                              tattACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc             >  1:2074263/4‑65 (MQ=255)
                                                              tattACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc             >  1:1915067/4‑65 (MQ=255)
                                                              tattACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc             >  1:177740/4‑65 (MQ=255)
                                                              tAATACTGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc             >  1:1088587/1‑65 (MQ=255)
                                                              tAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCgg                      <  1:957603/56‑1 (MQ=255)
                                                              tAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTcggcg       <  1:779333/71‑1 (MQ=255)
                                                                aTACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATc                                   >  1:697683/1‑41 (MQ=255)
                                                                   cAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGt                          >  1:1464904/1‑47 (MQ=255)
                                                                   cAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGagag  >  1:915768/1‑71 (MQ=255)
                                                                    |                                                                     
CCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAG  >  W3110S.gb/219619‑219756

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: