Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,018,630 G→A 40.0% D272N (GAC→AAC)  ssnA → predicted chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,018,6300GA40.0% 51.1 / 78.6 70D272N (GAC→AAC) ssnApredicted chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base A (19/9);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCAT  >  W3110S.gb/3018560‑3018698
                                                                      |                                                                    
aCGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATg                                                                      <  1:501808/71‑1 (MQ=255)
aCGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATg                                                                      <  1:1889429/71‑1 (MQ=255)
  gACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAAc                                                                    <  1:474546/71‑1 (MQ=255)
       gcgcAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAc                                                               <  1:1426451/71‑1 (MQ=255)
         gcAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCAccc                                                             >  1:2234726/1‑71 (MQ=255)
         gcAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCAccc                                                             >  1:405082/1‑71 (MQ=255)
          cAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAccc                                                             >  1:2036258/1‑70 (MQ=255)
          cAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCAccc                                                             >  1:152318/1‑70 (MQ=255)
             ttCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACAt                                                                  >  1:2439459/1‑62 (MQ=255)
              tCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACt                                                         >  1:1809761/1‑70 (MQ=255)
               cGATCTCTTCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTc                                                        >  1:2463032/1‑70 (MQ=255)
               cGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGa                                                                     >  1:127202/1‑57 (MQ=255)
               cGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGCCATCACCCTACTc                                                        >  1:2718864/1‑70 (MQ=255)
               cGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGAc                                                                    >  1:332435/1‑58 (MQ=255)
               cGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTc                                                        >  1:106588/1‑70 (MQ=255)
                gATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCa                                                       <  1:1894135/70‑1 (MQ=255)
                   ctcATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGAAATCACCCTACt                                                         >  1:2775907/1‑65 (MQ=255)
                   ctcATCGACAGCAAAACGCTGGTCGCTCATGGGCAGTACTTGTCGAAAGATGACATCACCCTACt                                                         >  1:208822/1‑65 (MQ=255)
                      aTCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGATAGATGACATCACCCTACTCAATCAgcg                                                >  1:1361580/1‑71 (MQ=255)
                      aTCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAgcg                                                >  1:1454138/1‑71 (MQ=255)
                      aTCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAgcg                                                >  1:1641296/1‑71 (MQ=255)
                      aTCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCa                                                                >  1:1419840/1‑55 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTGCTTGTCGAAAGATGACa                                                                   >  1:2363489/1‑49 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:1309952/1‑71 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGcgcg                                              >  1:139418/1‑70 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGa                                             >  1:2778625/1‑71 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGa                                             >  1:2642463/1‑71 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGa                                             >  1:1267907/1‑71 (MQ=255)
                         gACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGa                                             >  1:52402/1‑71 (MQ=255)
                              cAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGt                                        <  1:491750/71‑1 (MQ=255)
                               aaaaCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGt                                        <  1:2305361/70‑1 (MQ=255)
                                 aaCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTcc                                     <  1:2608643/71‑1 (MQ=255)
                                 aaCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTcc                                     <  1:679207/71‑1 (MQ=255)
                                 aaCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTcc                                     <  1:1139662/71‑1 (MQ=255)
                                  aCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCt                                    <  1:574741/71‑1 (MQ=255)
                                  aCGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCt                                    <  1:578424/71‑1 (MQ=255)
                                   cGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCaa                                                      >  1:2569867/1‑52 (MQ=255)
                                         tCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCAt                             <  1:1820612/71‑1 (MQ=255)
                                         tCGCTCATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGt                                        <  1:1270127/60‑1 (MQ=255)
                                          cGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:1411216/1‑54 (MQ=255)
                                              cATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTg                                   <  1:2000926/60‑1 (MQ=255)
                                              cATGGGCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGc                        >  1:2566698/1‑71 (MQ=255)
                                                 gggCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCg                     >  1:912494/1‑71 (MQ=255)
                                                 gggCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCg                     >  1:2365108/1‑71 (MQ=255)
                                                  ggCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGt                    >  1:1359684/1‑71 (MQ=255)
                                                  ggCTGTACTTGTCGAAAGATGACATCACCCAACTCAATCAGCGCGATGCGTTCCTGGT‑‑ATAACGCCCGtt                   >  1:892047/1‑70 (MQ=255)
                                                  ggCTGTACTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGt                    >  1:242893/1‑71 (MQ=255)
                                                       tACTTGTCGAAAGATGACATCACCCTACTCAATCAgcgc                                               <  1:751204/39‑1 (MQ=255)
                                                       tACTTGTCGAAAGATGACATCACCCTACTCAATCAgcgc                                               <  1:1417515/39‑1 (MQ=255)
                                                       tACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTcc                                     <  1:437403/49‑1 (MQ=255)
                                                       tACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCt                                    <  1:2510119/50‑1 (MQ=255)
                                                       tACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCaaa               <  1:735423/71‑1 (MQ=255)
                                                        aCTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGt                                        <  1:1173086/45‑1 (MQ=255)
                                                         cTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTg                                   >  1:1189885/1‑49 (MQ=255)
                                                         cTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTg                                   >  1:1697205/1‑49 (MQ=255)
                                                         cTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACa             >  1:2732303/1‑71 (MQ=255)
                                                         cTTGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACa             >  1:2249472/1‑71 (MQ=255)
                                                          ttGTCGAAAGATGACATCACCCTACTCAATCAgcgc                                               <  1:2492296/36‑1 (MQ=255)
                                                          ttGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATg                                           >  1:798768/1‑40 (MQ=255)
                                                          ttGTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGt                                        <  1:1546094/43‑1 (MQ=255)
                                                           tGTCGAAAGATAACATCACCCTACTCAATCAGCGCGa                                             >  1:2091488/1‑37 (MQ=255)
                                                            gTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:384343/1‑36 (MQ=255)
                                                            gTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:1749463/1‑36 (MQ=255)
                                                            gTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:529420/1‑36 (MQ=255)
                                                            gTCGAAAGATGACATCACCCTACTCAATCAGCGCGa                                             >  1:1996877/1‑36 (MQ=255)
                                                            gTCGAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATg           >  1:2558525/1‑70 (MQ=255)
                                                               gAAAGATAACATCACCCTACTCAATCAGGGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGaac        <  1:72822/70‑1 (MQ=255)
                                                               gAAAGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATaa                           >  1:1920110/1‑51 (MQ=255)
                                                                 aaGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGaacaac     <  1:340353/71‑1 (MQ=255)
                                                                 aaGATAACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGaacaac     <  1:680633/71‑1 (MQ=255)
                                                                     tGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCAt  <  1:192620/70‑1 (MQ=255)
                                                                     tGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCAt  <  1:2265796/70‑1 (MQ=255)
                                                                      |                                                                    
ACGACTGGCGCAATTCGATCTCATCGACAGCAAAACGCTGGTCGCTCATGGGCTGTACTTGTCGAAAGATGACATCACCCTACTCAATCAGCGCGATGCGTTCCTGGTGCATAACGCCCGTTCAAACATGAACAACCAT  >  W3110S.gb/3018560‑3018698

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: