Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 713,886 A→T 8.3% I160F (ATC→TTC)  seqA → regulatory protein for replication initiation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb713,8860AT 8.3% 369.2 / 18.0 130I160F (ATC→TTC) seqAregulatory protein for replication initiation
Reads supporting (aligned to +/- strand):  ref base A (66/53);  new base T (3/8);  total (69/61)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 138 total aligned reads displayed.

GCCGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAA  >  W3110S.gb/713816‑713956
                                                                      |                                                                      
gCCGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGa                                                                        <  1:1970520/71‑1 (MQ=255)
gCCGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGa                                                                        <  1:15019/71‑1 (MQ=255)
gCCGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGa                                                                        <  1:2637038/71‑1 (MQ=255)
 ccGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAATACCAACACCGGCCGTAAATGCAGCATGAt                                                                       >  1:1302320/1‑71 (MQ=255)
 ccGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGAt                                                                       >  1:3068457/1‑71 (MQ=255)
  cGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATc                                                                      <  1:5039507/71‑1 (MQ=255)
   gAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATct                                                                     <  1:4472542/71‑2 (MQ=255)
   gAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCg                                                                     <  1:1970306/71‑1 (MQ=255)
   gAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCg                                                                     <  1:4639039/71‑1 (MQ=255)
   gAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCg                                                                     <  1:584503/71‑1 (MQ=255)
     aaCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   <  1:892533/71‑1 (MQ=255)
     aaCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   <  1:4251314/71‑1 (MQ=255)
      aCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCg                                                                     <  1:3212738/68‑1 (MQ=255)
      aCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   <  1:2369446/70‑1 (MQ=255)
      aCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   >  1:2495775/1‑70 (MQ=255)
      aCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   >  1:4019895/1‑70 (MQ=255)
      aCATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   >  1:2265712/1‑70 (MQ=255)
        aTGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCTGCCGTAAATGCAGCATGATCGAacac                                                                <  1:4874905/71‑1 (MQ=255)
        aTGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAGTGCAGCATGATCGAacac                                                                <  1:3020473/71‑1 (MQ=255)
        aTGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   >  1:4078959/1‑68 (MQ=255)
         tgtgCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAAcaca                                                               <  1:3581831/71‑1 (MQ=255)
           tgCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACAcatc                                                             <  1:1824388/71‑1 (MQ=255)
             ccAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACAcatcat                                                           <  1:4355968/71‑1 (MQ=255)
             ccAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACAcatcat                                                           >  1:4547735/1‑71 (MQ=255)
              cAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATg                                                          <  1:467779/71‑1 (MQ=255)
                  cACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGTTCGAACACATCATGCAGt                                                      <  1:2157721/71‑1 (MQ=255)
                  cACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGTTCGAACACATCATGCAGt                                                      <  1:3329531/71‑1 (MQ=255)
                  cACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGc                                                         >  1:247309/1‑68 (MQ=255)
                   aCGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGTTCGAACACATCATGCAGTc                                                     <  1:3771358/71‑1 (MQ=255)
                   aCGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGTTCGAACACATCATGCAGTc                                                     <  1:4089760/71‑1 (MQ=255)
                   aCGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAaca                                                                 >  1:4412572/1‑59 (MQ=255)
                     gCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGa                                                   <  1:1366972/71‑1 (MQ=255)
                          aTTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGc                                                <  1:3639443/69‑1 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGTTCGAACACATCATGCAGTCGATGCAAt                                             >  1:525486/1‑71 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGaa                                                                   <  1:378857/49‑1 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGc                                                <  1:2345800/68‑1 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAt                                             >  1:2638511/1‑71 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAt                                             >  1:2684655/1‑71 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAt                                             >  1:3005977/1‑71 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAt                                             >  1:3513616/1‑71 (MQ=255)
                           ttGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAt                                             >  1:3874158/1‑71 (MQ=255)
                             gggTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAAtt                                            <  1:1190564/70‑1 (MQ=255)
                                tGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcg                                        <  1:1912498/71‑1 (MQ=255)
                                tGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcg                                        >  1:1352286/1‑71 (MQ=255)
                                 gATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGc                                                <  1:3269735/62‑1 (MQ=255)
                                 gATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcgg                                       >  1:5050955/1‑71 (MQ=255)
                                 gATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcgg                                       >  1:819048/1‑71 (MQ=255)
                                    caccaacaccaaCACCGGCCGTAAATGCAGCATGATCGaa                                                                   >  1:2550989/1‑40 (MQ=255)
                                     accaacaccaacACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCaa                                              >  1:3545727/1‑60 (MQ=255)
                                     accaacaccaacACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcggcgg                                    <  1:1078362/70‑1 (MQ=255)
                                        aacaccaacaccGGCCGTAAATGCAGCATGATCGAACAcatc                                                             >  1:1280377/1‑42 (MQ=255)
                                        aacaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTcc                                          >  1:1253509/1‑61 (MQ=255)
                                        aacaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:4271090/1‑71 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACAcatca                                                            <  1:78247/42‑1 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACAcatca                                                            <  1:2045793/42‑1 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATg                                                          >  1:93591/1‑44 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGt                                                      >  1:937175/1‑48 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:1384024/1‑70 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:3886122/1‑70 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:693046/1‑70 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:3907538/1‑70 (MQ=255)
                                         acaccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAatt                                >  1:3166521/1‑70 (MQ=255)
                                           accaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCaa                                              >  1:1492762/1‑54 (MQ=255)
                                           accaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCaa                                              >  1:2067521/1‑54 (MQ=255)
                                            ccaacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            <  1:4889491/71‑1 (MQ=255)
                                              aacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            >  1:1191009/1‑69 (MQ=255)
                                              aacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            >  1:2891836/1‑69 (MQ=255)
                                              aacaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAAttgattga                          <  1:2351695/71‑1 (MQ=255)
                                               acaccGGCCGTAAATGCAGCATGTTCGAACACATCATGc                                                         >  1:1965858/1‑39 (MQ=255)
                                               acaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCcggcgg                                    <  1:2282280/60‑1 (MQ=255)
                                               acaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            >  1:743844/1‑68 (MQ=255)
                                               acaccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            >  1:2386721/1‑68 (MQ=255)
                                               acaacGGCCGTAAATGCAGCATGATCGAACACATCATg                                                          >  1:457894/1‑38 (MQ=255)
                                                caccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTgaga                        <  1:3032024/71‑1 (MQ=255)
                                                 accGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTgag                         >  1:1280818/1‑69 (MQ=255)
                                                  ccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGaa                       >  1:517065/1‑70 (MQ=255)
                                                  ccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAg                      >  1:3638421/1‑71 (MQ=255)
                                                  ccGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAg                      <  1:887880/71‑1 (MQ=255)
                                                       cGTAAATGCAGCATGTTCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTg                 <  1:2799897/71‑1 (MQ=255)
                                                         tAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAattgatt                            <  1:4657320/58‑1 (MQ=255)
                                                         tAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCg               >  1:2837875/1‑71 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGaa                                  >  1:1299656/1‑50 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGaa                                  >  1:3893930/1‑50 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGaa                       >  1:3489511/1‑61 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCgg              >  1:2521121/1‑70 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCgg              >  1:3967046/1‑70 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCgg              >  1:3932866/1‑70 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCgg              <  1:1894445/70‑1 (MQ=255)
                                                           aaTGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGa             <  1:3810735/71‑1 (MQ=255)
                                                              gcagcaTGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAAc           >  1:4035500/1‑70 (MQ=255)
                                                               cagcaTGATCGAACACATCATGCAGTCGATGCAATTCCcggcgg                                    <  1:5094773/44‑1 (MQ=255)
                                                               cagcaTGATCGAACACATCATGCAGTCGATGCAATTCCcggcgg                                    <  1:2686144/44‑1 (MQ=255)
                                                               cagcaTGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTa         >  1:3843453/1‑71 (MQ=255)
                                                                agcaTGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTAt        <  1:4057591/71‑1 (MQ=255)
                                                                  caTGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCTGATCTATc       >  1:3337692/1‑70 (MQ=255)
                                                                     gATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTaa    >  1:794229/1‑70 (MQ=255)
                                                                      aTCGAACACATCATGCAGTCGATGCAATTCCcggcgg                                    <  1:2143112/37‑1 (MQ=255)
                                                                      aTCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTaaaa  >  1:3061664/1‑71 (MQ=255)
                                                                      aTCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTaaaa  >  1:3746914/1‑71 (MQ=255)
                                                                      aTCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAcg  <  1:730402‑M1/71‑3 (MQ=255)
                                                                      |                                                                      
GCCGAAACATGTGCCAGGCACGCCGTATTGGGTGATCACCAACACCAACACCGGCCGTAAATGCAGCATGATCGAACACATCATGCAGTCGATGCAATTCCCGGCGGAATTGATTGAGAAGGTTTGCGGAACTATCTAAAA  >  W3110S.gb/713816‑713956

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: