Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,696,825 T→C 19.6% 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,8250TC19.6% 97.1 / 20.0 46I321V (ATT→GTT) uidAbeta‑D‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base C (7/2);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCAT  >  W3110S.gb/1696755‑1696893
                                                                      |                                                                    
tGTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAAt                                                                      >  1:3885078/1‑71 (MQ=255)
  ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                     >  1:3053191/1‑70 (MQ=255)
  ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                     >  1:145213/1‑70 (MQ=255)
  ttCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATc                                                                     >  1:1628705/1‑70 (MQ=255)
    cATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAg                                                                  >  1:3638629/1‑71 (MQ=255)
    cATCTGCCCAGTCGAGCATCTCTTCACCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCa                                                                   <  1:3238624/70‑1 (MQ=255)
          cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAtt                                                            <  1:3362637/71‑1 (MQ=255)
          cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAtt                                                            <  1:416967/71‑1 (MQ=255)
          cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTc                                                                <  1:2517469/67‑1 (MQ=255)
          cccAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCAtt                                                            <  1:3276895/71‑1 (MQ=255)
              gTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCAtt                                                            <  1:1561536/67‑1 (MQ=255)
               tCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTc                                                                <  1:3319502/62‑1 (MQ=255)
                 gAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCg                                                      <  1:1113700/70‑1 (MQ=255)
                 gAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCg                                                      <  1:160000/70‑1 (MQ=255)
                    cATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTgg                                                   >  1:794387/1‑70 (MQ=255)
                    cATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGt                                                  >  1:3358735/1‑71 (MQ=255)
                    cATCTCTTCAGCGTAAGGGAATGCGAGGTACGGTAGGAGGTGGCCCCAATcc                                                                    <  1:19386/52‑1 (MQ=255)
                      tctcTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTc                                                                >  1:581348/1‑55 (MQ=255)
                      tctcTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTc                                                                <  1:2535130/55‑1 (MQ=255)
                        tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGt                                               >  1:2851240/1‑70 (MQ=255)
                        tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCATTAATGCGTGGTCGTg                                              >  1:442574/1‑71 (MQ=255)
                        tctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCATTAATGCGTGGTCGTg                                              >  1:1329599/1‑71 (MQ=255)
                         ctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCa                                                                   >  1:114081/1‑49 (MQ=255)
                         ctTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCa                                                                   >  1:1232477/1‑49 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGt                                                  >  1:421353/1‑63 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                           <  1:1864853/70‑1 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                           >  1:3363430/1‑70 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                           >  1:1239597/1‑70 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                           >  1:1359577/1‑70 (MQ=255)
                            cAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCAc                                           <  1:2198081/70‑1 (MQ=255)
                                   aGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCa                                   >  1:97360/1‑71 (MQ=255)
                                      gTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCATTaa                                                          >  1:1003259/1‑45 (MQ=255)
                                      gTAATGCGAGGTACGGTAGGAGTTGGCCCCAACCCAGTCCATTaa                                                          >  1:2042690/1‑45 (MQ=255)
                                           gCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTc                                                 <  1:989113/49‑1 (MQ=255)
                                            cGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATc                            >  1:941122/1‑69 (MQ=255)
                                            cGAGGTACGGTAGGAGTTGGCCACAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCg                           >  1:612725/1‑70 (MQ=255)
                                             gAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGa                          <  1:1326554/70‑1 (MQ=255)
                                                gTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGc                                                       <  1:22824/38‑1 (MQ=255)
                                                gTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATc                       >  1:315103/1‑70 (MQ=255)
                                                   cGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCttt                   >  1:3331752/1‑71 (MQ=255)
                                                     gTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGc                 <  1:1654447/71‑1 (MQ=255)
                                                      tAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGa                          <  1:2803205/61‑1 (MQ=255)
                                                       aGGAGTTGGCCCCAACCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCa               >  1:362813/1‑71 (MQ=255)
                                                           gTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGc                                             >  1:1636703/1‑37 (MQ=255)
                                                              ggCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAg         >  1:3092621/1‑70 (MQ=255)
                                                               gCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGt        <  1:1934973/70‑1 (MQ=255)
                                                                    aaTCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCAt  <  1:2382623/71‑1 (MQ=255)
                                                                      |                                                                    
TGTTCATCTGCCCAGTCGAGCATCTCTTCAGCGTAAGGGTAATGCGAGGTACGGTAGGAGTTGGCCCCAATCCAGTCCATTAATGCGTGGTCGTGCACCATCAGCACGTTATCGAATCCTTTGCCACGCAAGTCCGCAT  >  W3110S.gb/1696755‑1696893

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: