Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,696,825 T→C 10.1% I321V (ATT→GTT)  uidA ← beta‑D‑glucuronidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,696,8250TC10.1% 186.0 / 11.7 69I321V (ATT→GTT) uidAbeta‑D‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base T (39/23);  new base C (2/5);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCATCT  >  W3110S.gb/1696756‑1696895
                                                                     |                                                                      
gTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       >  1:1239529/1‑71 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATcc                                                                      <  1:3912498/71‑1 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATcc                                                                      <  1:2788034/71‑1 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       >  1:45707/1‑70 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       <  1:825288/70‑1 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       <  1:275421/70‑1 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       >  1:311986/1‑70 (MQ=255)
 ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       >  1:3932937/1‑70 (MQ=255)
  tCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                       >  1:2203389/1‑69 (MQ=255)
      cTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTc                                                                  >  1:826147/1‑70 (MQ=255)
         cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTc                                                                  <  1:1234424/67‑1 (MQ=255)
         cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTc                                                                  <  1:2521005/67‑1 (MQ=255)
         cccAGTCGAGCAGCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCAtt                                                              <  1:335289/71‑1 (MQ=255)
                gAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTaa                                                            >  1:4530710/1‑66 (MQ=255)
                   cATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGt                                                    >  1:354998/1‑71 (MQ=255)
                    aTCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATcc                                                                      >  1:4612348/1‑52 (MQ=255)
                     tctcTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTa                                                             <  1:897849/60‑1 (MQ=255)
                       tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAt                                                               >  1:2383480/1‑56 (MQ=255)
                       tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCg                                                        >  1:1632206/1‑63 (MQ=255)
                       tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCACTCCa                                                                >  1:369542/1‑55 (MQ=255)
                         ttCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGt                                                 <  1:234917/68‑1 (MQ=255)
                         ttCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTg                                                <  1:799090/69‑1 (MQ=255)
                         ttCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTg                                                <  1:2698285/69‑1 (MQ=255)
                         ttCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCACAATCCAGTCCATTAATGCGTGGTCGTg                                                <  1:2377080/69‑1 (MQ=255)
                          tCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCg                                                        <  1:453290/60‑1 (MQ=255)
                          tCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCATTAATg                                                          >  1:883988/1‑58 (MQ=255)
                           cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCa                                                                     >  1:460028/1‑46 (MQ=255)
                           cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGc                                               <  1:1658792/68‑1 (MQ=255)
                                tAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAtt                                                              <  1:1449583/48‑1 (MQ=255)
                                 aaGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                             >  1:3324583/1‑64 (MQ=255)
                                  aGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCa                                                                >  1:3857243/1‑44 (MQ=255)
                                  aGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGc                                      <  1:1469124/70‑1 (MQ=255)
                                  aGGGTAATGCGAGGTACGGTAGGAGTTGGCCCAATCCAGTCCATTaa                                                            >  1:3864166/1‑47 (MQ=255)
                                    ggTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAtt                                                              <  1:3398088/44‑1 (MQ=255)
                                      tAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGtt                                 <  1:3672408/71‑1 (MQ=255)
                                        aTGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTa                                <  1:4293749/70‑1 (MQ=255)
                                        aTGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTAt                               >  1:1023355/1‑71 (MQ=255)
                                        aTGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTAt                               >  1:1547091/1‑71 (MQ=255)
                                         tGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATc                              >  1:1580553/1‑71 (MQ=255)
                                         tGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATc                              >  1:2501916/1‑71 (MQ=255)
                                          gCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGc                                                         >  1:1226418/1‑43 (MQ=255)
                                          gCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTAt                               >  1:1168549/1‑69 (MQ=255)
                                          gCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATc                              >  1:4279550/1‑70 (MQ=255)
                                           cGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCAt                                          >  1:3588067/1‑57 (MQ=255)
                                           cGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCg                             >  1:3605500/1‑70 (MQ=255)
                                               gTACGGTAGGAGTTGGCCCCAATCCAGTCCATTaa                                                            >  1:4471149/1‑35 (MQ=255)
                                                tACGGTAGGAGTTGGCCCCAACCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGt                                  <  1:2033325/60‑1 (MQ=255)
                                                 aCGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCa                                        >  1:1846238/1‑53 (MQ=255)
                                                 aCGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCAc                                    >  1:252742/1‑57 (MQ=255)
                                                 aCGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCtt                      >  1:2229687/1‑71 (MQ=255)
                                                  cGGTAGGAGTTGGCCCCAACCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCttt                     <  1:2928594/71‑1 (MQ=255)
                                                    gTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAcc                                            >  1:4210225/1‑46 (MQ=255)
                                                     tAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATc                                         <  1:1065054/48‑1 (MQ=255)
                                                       ggAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCttt                     >  1:541161/1‑66 (MQ=255)
                                                          gTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATc                         >  1:1727878/1‑59 (MQ=255)
                                                           ttGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCt                       <  1:134285/60‑1 (MQ=255)
                                                           ttGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGc              <  1:1553963/69‑1 (MQ=255)
                                                           ttGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGc              <  1:720041/69‑1 (MQ=255)
                                                           ttGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGc              <  1:418945/69‑1 (MQ=255)
                                                             ggCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCttt                     >  1:2142212/1‑60 (MQ=255)
                                                             ggCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAg           >  1:1793277/1‑70 (MQ=255)
                                                             ggCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAg           <  1:1527175/70‑1 (MQ=255)
                                                             ggCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAg           >  1:3579155/1‑70 (MQ=255)
                                                                  cAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTg                    >  1:4369937/1‑56 (MQ=255)
                                                                  cAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCa     <  1:2290322/71‑1 (MQ=255)
                                                                  cAACCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGc      >  1:277437/1‑70 (MQ=255)
                                                                   aaTCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGaa                           >  1:1726497/1‑48 (MQ=255)
                                                                   aaTCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGaa                           >  1:1659522/1‑48 (MQ=255)
                                                                   aaTCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGa                            >  1:1884858/1‑47 (MQ=255)
                                                                     tCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCAc                <  1:3343021/57‑1 (MQ=255)
                                                                     tCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCATCt  >  1:2703954/1‑71 (MQ=255)
                                                                     |                                                                      
GTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCATCT  >  W3110S.gb/1696756‑1696895

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: