Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,437,526 A→G 7.6% T74T (ACT→ACC ytfF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,437,5260AG 7.6% 175.6 / 4.9 66T74T (ACT→ACCytfFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base G (2/3);  total (29/37)
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.61e-01

ACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGC  >  W3110S.gb/4437461‑4437594
                                                                 |                                                                    
aCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATa                                                                      >  1:772516/1‑66 (MQ=255)
aCAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATa                                                                      >  1:1514531/1‑66 (MQ=255)
 cAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAg                                                                 >  1:189061/1‑70 (MQ=255)
    gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCAt                                                              >  1:2334791/1‑70 (MQ=255)
    gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                                             <  1:777800/71‑1 (MQ=255)
    gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                                             <  1:1838151/71‑1 (MQ=255)
    gcgcgcCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATc                                                             <  1:3839367/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)
                gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                                 <  1:1666738/71‑1 (MQ=255)
                gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                                 <  1:3762403/71‑1 (MQ=255)
                gTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGc                                                 <  1:615972/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:3078672/71‑1 (MQ=255)
                      tGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGCCAGTcg                                           <  1:2542064/71‑1 (MQ=255)
                      tGGCACTTGCAAGGCAGAAGTAATAGATAAGGTTGCCCATCATGGTGAGCATCAAGGCGGTAAGCCAGTcg                                           <  1:999264/71‑1 (MQ=255)
                        gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:3742073/70‑1 (MQ=255)
                        gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          >  1:1302912/1‑70 (MQ=255)
                        gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:2626160/70‑1 (MQ=255)
                        gCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:45011/70‑1 (MQ=255)
                               cAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGcc                                  <  1:3749591/71‑1 (MQ=255)
                               cAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGcc                                  <  1:1377099/71‑1 (MQ=255)
                                aaGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCa                                 <  1:3258513/71‑1 (MQ=255)
                                    cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACt                              >  1:137144/1‑70 (MQ=255)
                                    cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg                             >  1:2202339/1‑71 (MQ=255)
                                    cAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg                             <  1:1996665/71‑1 (MQ=255)
                                     aGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTg                             <  1:2506414/70‑1 (MQ=255)
                                       aaGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacg                          <  1:2130992/71‑1 (MQ=255)
                                             tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:2155478/49‑1 (MQ=255)
                                             tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:3861104/49‑1 (MQ=255)
                                             tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTcgc                                          <  1:1283567/49‑1 (MQ=255)
                                             tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAc                                       <  1:1922906/52‑1 (MQ=255)
                                             tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAc                                       <  1:373621/52‑1 (MQ=255)
                                              aGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacg                          >  1:3163839/1‑64 (MQ=255)
                                               gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcgc                                     >  1:1698869/1‑52 (MQ=255)
                                               gATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcgc                                     >  1:2354769/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)
                                                   aGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGa                                        <  1:2923452/45‑1 (MQ=255)
                                                    ggTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAAc                               >  1:3346677/1‑53 (MQ=255)
                                                     gTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAc                                       <  1:1698433/44‑1 (MQ=255)
                                                       tGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTcc           >  1:1012880/1‑70 (MQ=255)
                                                        gCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgac                           >  1:1624650/1‑53 (MQ=255)
                                                        gCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgac                           >  1:3004263/1‑53 (MQ=255)
                                                        gCCCATCATAGTGAGCATCAAGGCGGTAAGCAAGTCGCGACGCGCCAACTgac                           >  1:1886980/1‑53 (MQ=255)
                                                         cccATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgcagacg                    <  1:2491155/59‑1 (MQ=255)
                                                            atcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgca                        <  1:292757/52‑1 (MQ=255)
                                                            atcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgc                         >  1:1141704/1‑51 (MQ=255)
                                                            atcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgc                         >  1:1966571/1‑51 (MQ=255)
                                                             tcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgcag                       <  1:1814970/52‑1 (MQ=255)
                                                             tcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGacgctcgc               <  1:533736/60‑1 (MQ=255)
                                                             tcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGacgcacgc               <  1:3175043/60‑1 (MQ=255)
                                                             tcatAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGacgcacgc               <  1:9312/60‑1 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcg                                      >  1:3213933/1‑35 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGAcg                                      >  1:3713299/1‑35 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAAc                               >  1:2480740/1‑42 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgcagacgca                  >  1:3562397/1‑55 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCccagccag    >  1:2899660/1‑69 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAgg   >  1:3195315/1‑70 (MQ=255)
                                                               atAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGc  <  1:1773835/71‑1 (MQ=255)
                                                                tAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTgacgca                        <  1:1142475/48‑1 (MQ=255)
                                                                 |                                                                    
ACAGGCGCGCCAGTACGTTGAATGGCACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGC  >  W3110S.gb/4437461‑4437594

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: