Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,115,678 A→T 8.3% R697R (CGT→CGA yghJ ← predicted inner membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,115,6780AT 8.3% 218.7 / 11.3 84R697R (CGT→CGAyghJpredicted inner membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (35/42);  new base T (2/5);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCA  >  W3110S.gb/3115612‑3115745
                                                                  |                                                                   
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:1056748/1‑69 (MQ=255)
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:3389919/1‑69 (MQ=255)
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:2876618/1‑69 (MQ=255)
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:1547495/1‑69 (MQ=255)
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:1459796/1‑69 (MQ=255)
tttGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACg                                                                   >  1:3164159/1‑69 (MQ=255)
    ttCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCgg                                                             <  1:901259/71‑1 (MQ=255)
    ttCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCgg                                                             <  1:2421646/71‑1 (MQ=255)
          tAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTg                                                        <  1:3045163/70‑1 (MQ=255)
          tAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGc                                                       <  1:1697788/71‑1 (MQ=255)
             gCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATa                                                    <  1:2989864/71‑1 (MQ=255)
                aTCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAgg                                                 <  1:1959426/71‑1 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTacac                                                                    >  1:1931337/1‑50 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCAc                                                                    >  1:3196135/1‑50 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGa                                                           >  1:928767/1‑59 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCa                                                      >  1:2974352/1‑64 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                               >  1:1287734/1‑71 (MQ=255)
                  cAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTCGTATCTGCACGCAGCGGGATTGCAt                                                     >  1:912968/1‑65 (MQ=255)
                   aGATCGGTCACATCCTGCTTAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATc                                              >  1:2220741/1‑71 (MQ=255)
                   aGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGAt                                               <  1:1998627/70‑1 (MQ=255)
                     aTCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGt                                            <  1:3591250/71‑1 (MQ=255)
                       cGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTa                                          >  1:1906553/1‑71 (MQ=255)
                       cGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTa                                          >  1:2465289/1‑71 (MQ=255)
                       cGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTa                                          >  1:2129219/1‑71 (MQ=255)
                        ggTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTa                                          <  1:1602838/70‑1 (MQ=255)
                        ggTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTAc                                         <  1:1874370/71‑1 (MQ=255)
                         gTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTAc                                         <  1:2799517/70‑1 (MQ=255)
                            acaTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                     >  1:1297999/1‑71 (MQ=255)
                            acaTCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                     >  1:910355/1‑71 (MQ=255)
                                 ctgctgAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGa                                 <  1:2145942/70‑1 (MQ=255)
                                    ctgAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTca                             <  1:1892937/71‑1 (MQ=255)
                                     tgAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTca                             <  1:2014781/70‑1 (MQ=255)
                                     tgAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTca                             <  1:1444824/70‑1 (MQ=255)
                                       aGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGt                                            >  1:545824/1‑53 (MQ=255)
                                       aGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGt                                            >  1:2678896/1‑53 (MQ=255)
                                       aGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTcac                            <  1:2898208/69‑1 (MQ=255)
                                       aGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTcac                            <  1:957597/69‑1 (MQ=255)
                                       aGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTcac                            <  1:2696601/69‑1 (MQ=255)
                                         tCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTcaca                           <  1:2799909/68‑1 (MQ=255)
                                         tCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTcaca                           <  1:1194582/68‑1 (MQ=255)
                                           agattCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTc                      <  1:2834054/68‑1 (MQ=255)
                                           aGCTTCGGTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGATCGt                                            >  1:385380/1‑49 (MQ=255)
                                           aGCTTCGGTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGATCGt                                            >  1:2410814/1‑49 (MQ=255)
                                           aGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATc                                              >  1:2599402/1‑47 (MQ=255)
                                           aGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATc                                              >  1:1866639/1‑47 (MQ=255)
                                           aGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGa                                 <  1:3429764/60‑1 (MQ=255)
                                            gCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTTAGTCACATATTc                      <  1:3716286/70‑1 (MQ=255)
                                            gCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTc                      <  1:325465/70‑1 (MQ=255)
                                            gCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTc                      <  1:1944635/70‑1 (MQ=255)
                                            gCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTc                      <  1:3075239/70‑1 (MQ=255)
                                              ttCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTg                                  <  1:3317700/56‑1 (MQ=255)
                                              ttCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCaa                    >  1:1003935/1‑70 (MQ=255)
                                               tCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAg                  >  1:1140424/1‑71 (MQ=255)
                                               tCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAg                  >  1:3198305/1‑71 (MQ=255)
                                               tCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAg                  >  1:2814184/1‑71 (MQ=255)
                                                cGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTAc                                         >  1:2106468/1‑47 (MQ=255)
                                                  gTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGAt                                               <  1:1885676/39‑1 (MQ=255)
                                                  gTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGAt                                               <  1:3619194/39‑1 (MQ=255)
                                                  gTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGAt                                               <  1:1298514/39‑1 (MQ=255)
                                                  gTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGAt                                               <  1:3813349/39‑1 (MQ=255)
                                                  gTTTGCTGGTATCTGCTCGCAGCGGGATTGCATAAGGAt                                               <  1:3390837/39‑1 (MQ=255)
                                                  gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTca                             <  1:2019202/57‑1 (MQ=255)
                                                   tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                     >  1:2585701/1‑48 (MQ=255)
                                                   tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                     >  1:1319697/1‑48 (MQ=255)
                                                   tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCAc                                     >  1:2884695/1‑48 (MQ=255)
                                                   tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGt              <  1:537013/71‑1 (MQ=255)
                                                   tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGt              <  1:2077304/71‑1 (MQ=255)
                                                   cttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGt              <  1:3859442/70‑1 (MQ=255)
                                                      gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTa            <  1:946509/70‑1 (MQ=255)
                                                      gCTGGTATCTGCACGCAGCGGGATTGCATAAGGACCGTTACCCACCTGAGTCACATATTCAAAGCCGTTa            <  1:1607940/70‑1 (MQ=255)
                                                        tGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTca                             <  1:2920947/51‑1 (MQ=255)
                                                          gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAgg         <  1:989590/69‑1 (MQ=255)
                                                          gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGa        <  1:874501/70‑1 (MQ=255)
                                                          gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGa        <  1:1155445/70‑1 (MQ=255)
                                                           tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCGAAGCCGTTAAGGAt       >  1:2801051/1‑70 (MQ=255)
                                                            aTCTGCACGCAGCGGGATTGCATAAGGCTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGa        >  1:2438877/1‑68 (MQ=255)
                                                            aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGAt       >  1:2100905/1‑69 (MQ=255)
                                                             tCTGCACGCAGCGGGATTGCATAAGGATCGTTAcc                                        <  1:3254793/35‑1 (MQ=255)
                                                             tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcag    >  1:1887576/1‑71 (MQ=255)
                                                               tGCACGCAGCGGGATTGCATAAGGATCGTTACCCAcc                                    >  1:1177336/1‑37 (MQ=255)
                                                                 cACGCAGCGGGATTGCATAAGTATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGTATcag    <  1:3357055/67‑1 (MQ=255)
                                                                 cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATAtt                       >  1:541780/1‑48 (MQ=255)
                                                                 cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcagca  <  1:3715046/69‑1 (MQ=255)
                                                                 cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcag    <  1:837036/67‑1 (MQ=255)
                                                                 cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATcag    <  1:1022453/67‑1 (MQ=255)
                                                                  |                                                                   
TTTGTTCAGATAGGCGATCAGATCGGTCACATCCTGCTGAGTCAGCTTCGGTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCA  >  W3110S.gb/3115612‑3115745

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: