Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,903,275 T→C 7.0% I9V (ATT→GTT)  yoaE ← fused predicted membrane proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,903,2750TC 7.0% 161.9 / 2.9 57I9V (ATT→GTT) yoaEfused predicted membrane proteins
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAAGACCAGGTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGAAATAAT  >  W3110S.gb/1903207‑1903344
                                                                    |                                                                     
gATGAAGACCAGGTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTg                                                                     <  1:1365880/71‑1 (MQ=255)
gATGAAGACCAGGTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTg                                                                     <  1:2051521/71‑1 (MQ=255)
 aTGAAGACCAGGTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAAtt                                                                      <  1:76001/69‑1 (MQ=255)
           ggTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCAtt                                                          <  1:45386/71‑1 (MQ=255)
            gTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAg                                                                   <  1:1727440/61‑1 (MQ=255)
            gTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCAt                                                           >  1:3815944/1‑69 (MQ=255)
              tATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAgg                                                                  >  1:1050097/1‑60 (MQ=255)
              tATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAgg                                                                  >  1:4910818/1‑60 (MQ=255)
              tATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAgg                                                                  >  1:2430018/1‑60 (MQ=255)
               aTCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAg                                                                   >  1:5047042/1‑58 (MQ=255)
                tCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGaa                                                     >  1:4756946/1‑71 (MQ=255)
                 cGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAAtt                                                                      >  1:1520364/1‑53 (MQ=255)
                  gATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTATCCCCGCCCAAATTGAGGGGTCCATTAATAAtt                                                   <  1:4638522/71‑1 (MQ=255)
                  gATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTTAGaa                                                     >  1:1231024/1‑69 (MQ=255)
                  gATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAAtt                                                   >  1:2261841/1‑71 (MQ=255)
                    tACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTcc                                                 <  1:2607597/71‑1 (MQ=255)
                    tACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTcc                                                 <  1:1303486/71‑1 (MQ=255)
                     aCCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAAtt                                                                      <  1:3108601/49‑1 (MQ=255)
                     aCCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAAtt                                                                      <  1:3751334/49‑1 (MQ=255)
                      cccAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTg                                                                     <  1:207523/49‑1 (MQ=255)
                      cccAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTg                                                                     <  1:2768934/49‑1 (MQ=255)
                      cccAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTg                                                                     <  1:871738/49‑1 (MQ=255)
                      cccAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCa                                                <  1:827887/70‑1 (MQ=255)
                         aGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAACTGAGGGGTCCATTAAGAATTCCATGAc                                            <  1:1653568/71‑1 (MQ=255)
                         aGCACAATTTCGAGAACAACAAGCGTGAGTAGCCACGCCCAAACTGAGGGGTCCATTAAGAATTCCATGAc                                            <  1:1678521/71‑1 (MQ=255)
                          gCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGAc                                            >  1:3613464/1‑70 (MQ=255)
                          gCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGAc                                            >  1:3930207/1‑70 (MQ=255)
                          gCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACa                                           <  1:2736998/71‑1 (MQ=255)
                          gCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACa                                           <  1:3343436/71‑1 (MQ=255)
                           cacaATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGAc                                            >  1:3619030/1‑69 (MQ=255)
                             caATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGa                                                                    >  1:4728303/1‑43 (MQ=255)
                               aTTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGAc                                            >  1:3652324/1‑65 (MQ=255)
                               aTTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGAc                                            >  1:565686/1‑65 (MQ=255)
                                 ttCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAg                                                                   >  1:1800711/1‑40 (MQ=255)
                                    gagaACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCt                                 >  1:1473805/1‑71 (MQ=255)
                                    gagaACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCt                                 >  1:896472/1‑71 (MQ=255)
                                    gagaACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCt                                 >  1:1626326/1‑71 (MQ=255)
                                        acaacaAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTc                                                              <  1:4971827/38‑1 (MQ=255)
                                         caacaaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGaa                                                     >  1:4847116/1‑46 (MQ=255)
                                         caacaaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGaa                                                     >  1:4257558/1‑46 (MQ=255)
                                         caacaaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAgg                            <  1:1067350/71‑1 (MQ=255)
                                         caacaaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAgg                            <  1:3300775/71‑1 (MQ=255)
                                         caacaaGCGTGAGTAGCCCCGCCCAAACTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAg                             >  1:615175/1‑70 (MQ=255)
                                          aacaaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTc                                                  <  1:1990217/48‑1 (MQ=255)
                                          aacaaGCGTGAGTAGCCCCGCCCAAACTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAgg                            >  1:4694707/1‑70 (MQ=255)
                                            caaGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAAt                         <  1:2514203/71‑1 (MQ=255)
                                               gCGTGAGTAGCCCCGCCCAAATTGAGTGGTCCATTAAGAATTCCATGACa                                           <  1:2477429/50‑1 (MQ=255)
                                               gCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGc                      >  1:4273081/1‑71 (MQ=255)
                                                cGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGc                      >  1:609582/1‑70 (MQ=255)
                                                cGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGc                      >  1:2401615/1‑70 (MQ=255)
                                                 gTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGa                                                      >  1:899467/1‑37 (MQ=255)
                                                    aGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTc                 <  1:3508340/71‑1 (MQ=255)
                                                      tAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGaa                          >  1:396905/1‑60 (MQ=255)
                                                           ccGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGc            >  1:1498108/1‑69 (MQ=255)
                                                               ccAAATTGAGGGGTCCATTAAGAATTCCATGACaac                                         >  1:289228/1‑35 (MQ=255)
                                                               ccAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATa                        <  1:404771/53‑1 (MQ=255)
                                                               ccAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGaaa      >  1:283264/1‑71 (MQ=255)
                                                                cAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAgg                            <  1:1880946/48‑1 (MQ=255)
                                                                 aaaTTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGAaata    <  1:3882403/71‑1 (MQ=255)
                                                                   aTTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGAaataat  >  1:2779855/1‑71 (MQ=255)
                                                                   aTTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGAaataat  >  1:170578/1‑71 (MQ=255)
                                                                    ttGAGGGGTCGATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCt                     <  1:890327/51‑1 (MQ=255)
                                                                    ttGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCt                     <  1:161042/51‑1 (MQ=255)
                                                                    |                                                                     
GATGAAGACCAGGTTATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTCCATGACAAGCTCCTGCTTAAGGAATAGCTATTCGACGCCAGAAATAAT  >  W3110S.gb/1903207‑1903344

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: