Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,018,853 C→T 28.1% T346I (ACT→ATT)  ssnA → predicted chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,018,8530CT28.1% 95.9 / 39.0 57T346I (ACT→ATT) ssnApredicted chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base T (6/10);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGC  >  W3110S.gb/3018788‑3018922
                                                                 |                                                                     
tGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTAACg                                                                  >  1:3845165/1‑71 (MQ=255)
tGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTAACg                                                                  >  1:228362/1‑71 (MQ=255)
   cttcttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGt                                                                >  1:1097655/1‑70 (MQ=255)
      cttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTaa                                                              <  1:4144940/69‑1 (MQ=255)
       tttAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACg                                                            >  1:3618630/1‑70 (MQ=255)
       tttAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACg                                                            >  1:2610348/1‑70 (MQ=255)
          aaaCATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACt                                                        >  1:4768238/1‑71 (MQ=255)
           aaCATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGaa                                                          >  1:4838493/1‑68 (MQ=255)
           aaCATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACt                                                        <  1:2755191/70‑1 (MQ=255)
             cATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACtgat                                                     <  1:2905836/71‑1 (MQ=255)
                       cGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCaa                                            <  1:1434145/70‑1 (MQ=255)
                         ggtggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACtt                                         >  1:2076253/1‑71 (MQ=255)
                         ggtggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACtt                                         >  1:105339/1‑71 (MQ=255)
                          gtggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACttt                                        <  1:4469661/71‑1 (MQ=255)
                           tggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTg                                       <  1:334725/71‑1 (MQ=255)
                             gtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAAc                                                                   <  1:2190220/41‑1 (MQ=255)
                              tCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTaa                                                                    <  1:1561752/39‑1 (MQ=255)
                              tCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTaa                                                                    <  1:1463439/39‑1 (MQ=255)
                               ccGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACtga                                                      >  1:3960875/1‑52 (MQ=255)
                               ccGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACtga                                                      >  1:304239/1‑52 (MQ=255)
                                cGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACg                                                            <  1:2329325/45‑1 (MQ=255)
                                cGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCATCTTTGgcgc                                   >  1:1782115/1‑70 (MQ=255)
                                cGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGgcgc                                   >  1:4646618/1‑70 (MQ=255)
                                cGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGgcgc                                   >  1:3103423/1‑70 (MQ=255)
                                 gCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGa                                 <  1:4617280/71‑1 (MQ=255)
                                  cTGTGGCCTGACAGTTTTGCCAAAGCCCTGATTAACGGTAACGa                                                           <  1:880134/44‑1 (MQ=255)
                                    gtgGCCTGACAGTTTTGCCAAAGCCCTGATTAACGGTAACGAAc                                                         >  1:4624166/1‑44 (MQ=255)
                                    gtgGCCTGACAGTTTTGCCAAAGCCCTGATTAACGGTAACGAAc                                                         >  1:1732153/1‑44 (MQ=255)
                                     tgGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAt                              <  1:4724382/70‑1 (MQ=255)
                                       gCCTGACAGTTTTGCCAAAGCCCTGACTAACTGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTg                           <  1:3386002/71‑1 (MQ=255)
                                       gCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAttt                            <  1:4977999/70‑1 (MQ=255)
                                         cTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACg                                                            <  1:829891/36‑1 (MQ=255)
                                           gACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTg                                       <  1:1870559/55‑1 (MQ=255)
                                           gACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTg                                       <  1:1348996/55‑1 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCtt                      >  1:3599186/1‑71 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCt                       <  1:134653/70‑1 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCtt                      <  1:564749/71‑1 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCtt                      <  1:509010/71‑1 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCtt                      <  1:3353505/71‑1 (MQ=255)
                                            aCAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCtt                      <  1:1488088/71‑1 (MQ=255)
                                              aGTTTTGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCt                    >  1:4548158/1‑71 (MQ=255)
                                              aGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCt                    >  1:4972160/1‑71 (MQ=255)
                                              aGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCt                    >  1:485071/1‑71 (MQ=255)
                                                ttttGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGaaa                               <  1:1316575/58‑1 (MQ=255)
                                                ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTc                     <  1:4721406/68‑1 (MQ=255)
                                                  ttGCCAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGaaa                               <  1:455739/56‑1 (MQ=255)
                                                  ttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAttt                            >  1:4628143/1‑59 (MQ=255)
                                                  ttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAttt                            >  1:2995925/1‑59 (MQ=255)
                                                   tGCCAAAGCCCTGACTAACGCTAACGAACTGATGAGCCGCAACTTTGGCGCGaaa                               <  1:4554571/55‑1 (MQ=255)
                                                     ccAAAGCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGc              <  1:1983219/70‑1 (MQ=255)
                                                        aaGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGt          <  1:694190/71‑1 (MQ=255)
                                                         aGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAttt                            >  1:4284055/1‑52 (MQ=255)
                                                         aGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAAttt                            >  1:706077/1‑52 (MQ=255)
                                                          gCCCTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTa        <  1:4177414/71‑1 (MQ=255)
                                                            ccTGATTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACa      <  1:3318998/71‑1 (MQ=255)
                                                                aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGTAAGCCGGTTATAAAg   <  1:3184876/70‑1 (MQ=255)
                                                                aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGc  >  1:4721312/1‑71 (MQ=255)
                                                                 |                                                                     
TGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGC  >  W3110S.gb/3018788‑3018922

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: