Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,594,973 A→G 10.2% pseudogene (807/1401 nt) ydeU ← ECK1502:JW1502:b1509; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,594,9730AG10.2% 150.6 / 9.5 59pseudogene (807/1401 nt)ydeUECK1502:JW1502:b1509; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/17);  new base G (4/2);  total (40/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGAT  >  W3110S.gb/1594909‑1595040
                                                                |                                                                   
tATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAg                                                               >  1:714935/1‑71 (MQ=255)
tATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAg                                                               >  1:694386/1‑71 (MQ=255)
    aaaCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATa                                                                   >  1:2289851/1‑63 (MQ=255)
    aaaCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATAc                                                                  >  1:2979198/1‑64 (MQ=255)
    aaaCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATAc                                                                  >  1:2983935/1‑64 (MQ=255)
     aaCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTcgc                                                          <  1:1905400/71‑1 (MQ=255)
      aCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTcgcg                                                         <  1:894632/71‑1 (MQ=255)
        cAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTa                                                       <  1:239841/71‑1 (MQ=255)
        cAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTa                                                       <  1:3553586/71‑1 (MQ=255)
           ttATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTa                                                       <  1:1919783/68‑1 (MQ=255)
             ataCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAGGTCGCGTACAGcc                                                  >  1:1626679/1‑71 (MQ=255)
             ataCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                  >  1:1957080/1‑71 (MQ=255)
             ataCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                  >  1:2341776/1‑71 (MQ=255)
             ataCAGAGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGca                                                  >  1:9095/1‑70 (MQ=255)
              taCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                  >  1:2014111/1‑70 (MQ=255)
               aCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                  >  1:1742898/1‑69 (MQ=255)
                   cgcCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATa                                                                   >  1:715466/1‑48 (MQ=255)
                         gCTGTCAACATAAGCGCCGTCTTATTCGCATCGTTCTGa                                                                     <  1:371447/39‑1 (MQ=38)
                           tGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATa                                                                   >  1:772609/1‑40 (MQ=255)
                           tGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATa                                                                   >  1:1115320/1‑40 (MQ=255)
                           tGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATa                                                                   >  1:2888725/1‑40 (MQ=255)
                             tCAACGTAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  >  1:4029839/1‑71 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:1183239/71‑1 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:1417692/71‑1 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGc                                   >  1:2388932/1‑70 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:2653178/71‑1 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:819788/71‑1 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:2877547/71‑1 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  >  1:2340690/1‑71 (MQ=255)
                             tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:773579/71‑1 (MQ=255)
                              cAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  <  1:1410047/70‑1 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCa                                                                >  1:950701/1‑39 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCa                                                                >  1:3120767/1‑39 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCa                                                                >  1:2681959/1‑39 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGt                                                        >  1:2261868/1‑47 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCt                                           >  1:1638405/1‑60 (MQ=255)
                               aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                  >  1:1637782/1‑69 (MQ=255)
                                    aaGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATgcgc                           >  1:2539925/1‑71 (MQ=255)
                                     aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTg                                          >  1:2815375/1‑55 (MQ=255)
                                     aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                          >  1:234052/1‑71 (MQ=255)
                                     aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                          >  1:54839/1‑71 (MQ=255)
                                     aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                          >  1:2394252/1‑71 (MQ=255)
                                       cgcCGGTCTTATTCGCATCTTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCa                        <  1:3332752/71‑1 (MQ=255)
                                        gcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCa                        <  1:1577995/70‑1 (MQ=255)
                                        gcCGGTCTTATTCGAATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCa                        <  1:538898/70‑1 (MQ=255)
                                          cGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGc                                   <  1:3254642/57‑1 (MQ=255)
                                           ggTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAAc                                      >  1:331544/1‑53 (MQ=255)
                                                 aTTGGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTta               >  1:3698374/1‑70 (MQ=255)
                                                 aTTCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAGCGCTGTAAcc                                     >  1:1552408/1‑48 (MQ=255)
                                                 aTTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGAtt                 >  1:2114354/1‑68 (MQ=255)
                                                 aTTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTta               >  1:4052495/1‑70 (MQ=255)
                                                    cGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCt                                           >  1:309453/1‑39 (MQ=255)
                                                    cGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCt                                           >  1:1697324/1‑39 (MQ=255)
                                                    cGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCt                                           >  1:3557511/1‑39 (MQ=255)
                                                          gTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac     >  1:522123/1‑71 (MQ=255)
                                                          gTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac     <  1:727833/71‑1 (MQ=255)
                                                          gTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac     <  1:3737855/71‑1 (MQ=255)
                                                          gTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac     <  1:3699518/71‑1 (MQ=255)
                                                             cTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                          >  1:2276251/1‑47 (MQ=255)
                                                             cTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCAATCCGATTTATAACCCACACGAt  >  1:3565085/1‑71 (MQ=255)
                                                                |                                                                   
TATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGAT  >  W3110S.gb/1594909‑1595040

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: