Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,437,511 G→A 7.8% L79L (CTC→CTT ytfF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,437,5110GA 7.8% 188.5 / 4.4 64L79L (CTC→CTTytfFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base A (2/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CCGATAATCATCGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACG  >  W3110S.gb/4437443‑4437580
                                                                    |                                                                     
ccGATAATCATCGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAg                                                                     >  1:1046024/1‑71 (MQ=255)
  gATAATCATCGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGt                                                                   >  1:880129/1‑71 (MQ=255)
     aaTCATCGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGc                                                                >  1:1710666/1‑71 (MQ=255)
        catcGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCca                                                             >  1:3833947/1‑71 (MQ=255)
         atcGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcat                                                            >  1:2170662/1‑71 (MQ=255)
         atcGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcat                                                            >  1:166268/1‑71 (MQ=255)
         atcGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcat                                                            >  1:1606715/1‑71 (MQ=255)
         atcGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCca                                                             >  1:3221721/1‑70 (MQ=255)
         atcGTGGAAACAGGAGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcat                                                            >  1:2913676/1‑71 (MQ=255)
            gTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:1530110/70‑1 (MQ=255)
            gTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:703577/70‑1 (MQ=255)
            gTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:1525697/70‑1 (MQ=255)
            gTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          >  1:1344339/1‑70 (MQ=255)
                 aaCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:1549480/65‑1 (MQ=255)
                 aaCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:2198211/65‑1 (MQ=255)
                 aaCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:1321632/65‑1 (MQ=255)
                 aaCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCcatca                                                          <  1:1088543/65‑1 (MQ=255)
                  aCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATa                                                        >  1:1514531/1‑66 (MQ=255)
                  aCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATa                                                        >  1:772516/1‑66 (MQ=255)
                   cAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAg                                                   >  1:189061/1‑70 (MQ=255)
                      gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCAt                                                >  1:2334791/1‑70 (MQ=255)
                      gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                               <  1:3839367/71‑1 (MQ=255)
                      gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                               <  1:1838151/71‑1 (MQ=255)
                      gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                               <  1:777800/71‑1 (MQ=255)
                       cgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                               >  1:279377/1‑70 (MQ=255)
                              gTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGt                                       <  1:830044/71‑1 (MQ=255)
                              gTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGt                                       <  1:512999/71‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGa                                                                      <  1:3060139/36‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   <  1:615972/71‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   <  1:1666738/71‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   <  1:3762403/71‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   <  1:817051/71‑1 (MQ=255)
                                  gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   <  1:2077232/71‑1 (MQ=255)
                                    tGAATGGCACTTGCAAGTCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGc                                   >  1:2498504/1‑69 (MQ=255)
                                    tGAATGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGc                                   >  1:759624/1‑69 (MQ=255)
                                    tGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                   >  1:3495774/1‑69 (MQ=255)
                                     gAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATa                                                        >  1:3842781/1‑47 (MQ=255)
                                        tGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGCCAGTcg                             <  1:999264/71‑1 (MQ=255)
                                        tGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGCCAGTcg                             <  1:3078672/71‑1 (MQ=255)
                                        tGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGCCAGTcg                             <  1:2542064/71‑1 (MQ=255)
                                          gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:2626160/70‑1 (MQ=255)
                                          gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            >  1:1302912/1‑70 (MQ=255)
                                          gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:3742073/70‑1 (MQ=255)
                                          gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:45011/70‑1 (MQ=255)
                                                 cAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGcc                    <  1:1377099/71‑1 (MQ=255)
                                                 cAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGcc                    <  1:3749591/71‑1 (MQ=255)
                                                  aaGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCa                   <  1:3258513/71‑1 (MQ=255)
                                                      cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACt                >  1:137144/1‑70 (MQ=255)
                                                      cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg               <  1:1996665/71‑1 (MQ=255)
                                                      cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg               >  1:2202339/1‑71 (MQ=255)
                                                       aGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg               <  1:2506414/70‑1 (MQ=255)
                                                         aaGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacg            <  1:2130992/71‑1 (MQ=255)
                                                               tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:1283567/49‑1 (MQ=255)
                                                               tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:2155478/49‑1 (MQ=255)
                                                               tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                            <  1:3861104/49‑1 (MQ=255)
                                                               tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAc                         <  1:373621/52‑1 (MQ=255)
                                                               tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAc                         <  1:1922906/52‑1 (MQ=255)
                                                                aGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacg            >  1:3163839/1‑64 (MQ=255)
                                                                 gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcgc                       >  1:2354769/1‑52 (MQ=255)
                                                                 gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcgc                       >  1:1698869/1‑52 (MQ=255)
                                                                 gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcgc                       >  1:3616567/1‑52 (MQ=255)
                                                                 gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgcagacgca    >  1:2570562/1‑71 (MQ=255)
                                                                   tGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGacgcacg  >  1:1853272/1‑71 (MQ=255)
                                                                   tGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGacgcacg  >  1:3600983/1‑71 (MQ=255)
                                                                    |                                                                     
CCGATAATCATCGTGGAAACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACG  >  W3110S.gb/4437443‑4437580

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: