Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,065,341 G→T 38.7% I163I (ATC→ATA ygfI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,065,3410GT38.7% 95.8 / 53.7 62I163I (ATC→ATAygfIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (15/9);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTTGCTCTTGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAAC  >  W3110S.gb/3065272‑3065407
                                                                     |                                                                  
ctcttgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                   <  1:3111293/71‑1 (MQ=255)
    tgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                               >  1:1390218/1‑71 (MQ=255)
    tgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                               >  1:4314267/1‑71 (MQ=255)
     gctcttgTCGCTCTTCATAACGGCGATAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGt                                                               >  1:790438/1‑70 (MQ=255)
     gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTtata                                                                 <  1:1462132/68‑1 (MQ=255)
     gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGt                                                               <  1:877697/70‑1 (MQ=255)
     gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGt                                                               >  1:2007472/1‑70 (MQ=255)
     gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGt                                                               >  1:4454670/1‑70 (MQ=255)
         ttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAg                                                            <  1:2226425/69‑1 (MQ=255)
           gTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAg                                                            >  1:3043329/1‑67 (MQ=255)
           gTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGttt                                                        >  1:2645369/1‑71 (MQ=255)
              gCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTc                                                      >  1:4358246/1‑70 (MQ=255)
              gCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTc                                                      <  1:3413053/70‑1 (MQ=255)
                  ttCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAg                                                   <  1:615341/69‑1 (MQ=255)
                   tCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAg                                                                <  1:617716/55‑1 (MQ=255)
                   tCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGtt                                                 >  1:4637804/1‑70 (MQ=255)
                   tCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTc                                                >  1:3990896/1‑71 (MQ=255)
                    cATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAg                                                                <  1:1675802/54‑1 (MQ=255)
                    cATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAg                                                                <  1:1669902/54‑1 (MQ=255)
                       aaCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTg                                             >  1:4204939/1‑70 (MQ=255)
                       aaCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTg                                             >  1:3542750/1‑70 (MQ=255)
                        aCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGaa                                           >  1:251804/1‑71 (MQ=255)
                        aCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGaa                                           <  1:1110611/71‑1 (MQ=255)
                         cGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAAc                                          <  1:4528322/71‑1 (MQ=255)
                            cgagaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCg                                       <  1:4414084/71‑1 (MQ=255)
                            cgagaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCg                                       <  1:810570/71‑1 (MQ=255)
                               gaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                     >  1:3566019/1‑70 (MQ=255)
                                aaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                   <  1:3586965/71‑1 (MQ=255)
                                aaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                   <  1:486921/71‑1 (MQ=255)
                                aaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                   <  1:3209333/71‑1 (MQ=255)
                                aaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                   <  1:2981661/71‑1 (MQ=255)
                                 aacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAg                                  <  1:2323667/71‑1 (MQ=255)
                                  acgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                     >  1:4165088/1‑67 (MQ=255)
                                   cgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGa                                 >  1:1522040/1‑70 (MQ=255)
                                    gagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTcc                                                     >  1:5002840/1‑49 (MQ=255)
                                    gagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTcc                                                     >  1:387481/1‑49 (MQ=255)
                                    gagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTcc                                                     >  1:2773998/1‑49 (MQ=255)
                                    gagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               >  1:2635893/1‑71 (MQ=255)
                                     agaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               >  1:4668131/1‑70 (MQ=255)
                                      gaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACa                             >  1:1776160/1‑71 (MQ=255)
                                          gCCTAATCTCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCt                         >  1:4631960/1‑71 (MQ=255)
                                           ccTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtg                        >  1:2313324/1‑71 (MQ=255)
                                           ccTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtg                        <  1:4497840/71‑1 (MQ=255)
                                            cTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtgt                       >  1:2681399/1‑71 (MQ=255)
                                              aaTCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGc                                    >  1:4541392/1‑56 (MQ=255)
                                                tCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtgtg                      <  1:1070543/68‑1 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATg                           >  1:390043/1‑61 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATg                           >  1:2831926/1‑61 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTggg                    >  1:1280783/1‑68 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTggg                    >  1:675746/1‑68 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTggg                    >  1:3747744/1‑68 (MQ=255)
                                                  gCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACa                             >  1:4026059/1‑59 (MQ=255)
                                                     gCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                   <  1:931000/50‑1 (MQ=255)
                                                        gtggtAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               <  1:2763387/51‑1 (MQ=255)
                                                        gtggtAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               <  1:5065338/51‑1 (MQ=255)
                                                           gtAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGa                                 >  1:5033988/1‑46 (MQ=255)
                                                           gtAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATaaa         <  1:3089536/70‑1 (MQ=255)
                                                               ttCACTTATAGTCAGGTTTTCCAGTTCCTTAACTCGTGCAGATACATGCTGTGGGGAGATAATaaa         <  1:4516027/66‑1 (MQ=255)
                                                                 cACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGa                                 >  1:2484185/1‑40 (MQ=255)
                                                                 cACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAAc  >  1:5007124/1‑71 (MQ=255)
                                                                 cACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAAc  >  1:2354077/1‑71 (MQ=255)
                                                                  aCTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATaaaa        <  1:1733570/64‑1 (MQ=255)
                                                                  aCTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATaaaa        <  1:2849075/64‑1 (MQ=255)
                                                                     gATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               <  1:4944029/38‑1 (MQ=255)
                                                                     gATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               <  1:275600/38‑1 (MQ=255)
                                                                     gATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                               <  1:3783261/38‑1 (MQ=255)
                                                                     |                                                                  
CTCTTGCTCTTGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAAC  >  W3110S.gb/3065272‑3065407

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: