Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,656,712 C→A 15.4% Q48H (CAG→CAT)  rspA ← predicted dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,656,7120CA15.4% 173.2 / 11.4 58Q48H (CAG→CAT) rspApredicted dehydratase
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base A (5/4);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCG  >  W3110S.gb/1656644‑1656778
                                                                    |                                                                  
taGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATcc                                                                    <  1:4254781/69‑1 (MQ=255)
taGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTg                                                                  >  1:193828/1‑71 (MQ=255)
   aaaaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCa                                                                >  1:2214278/1‑70 (MQ=255)
    aaaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCaaa                                                              <  1:4091908/71‑1 (MQ=255)
    aaaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCaaa                                                              <  1:4263118/71‑1 (MQ=255)
     aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCaaa                                                              >  1:3302962/1‑70 (MQ=255)
     aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAAt                                                             <  1:821310/71‑1 (MQ=255)
           ccAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAgg                                                       >  1:2854930/1‑71 (MQ=255)
               atatCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATaa                                                           <  1:2161663/63‑1 (MQ=255)
                tatCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACg                                                  <  1:1178709/71‑1 (MQ=255)
                 atCTTCGATACGGTGCGCATCGCGCCCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACg                                                  >  1:2543089/1‑70 (MQ=255)
                 atCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAaga                                                         >  1:2248278/1‑63 (MQ=255)
                 atCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGAGATCCTGCAAATAAGAGGCCACg                                                  >  1:2157170/1‑70 (MQ=255)
                      cGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAg                                             <  1:4147381/70‑1 (MQ=255)
                       gATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCt                                           <  1:3273400/71‑1 (MQ=255)
                        aTACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATaa                                                           <  1:3260042/54‑1 (MQ=255)
                          aCGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCAc                                        >  1:216874/1‑71 (MQ=255)
                           cGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATaa                                                           <  1:3930785/51‑1 (MQ=255)
                             gTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTc                                     >  1:3138301/1‑71 (MQ=255)
                             gTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTc                                     >  1:728905/1‑71 (MQ=255)
                              tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTcc                                    <  1:2671344/71‑1 (MQ=255)
                                  cATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTg                                <  1:4054283/71‑1 (MQ=255)
                                  cATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTg                                <  1:4059503/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:2384276/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:4002499/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:2035163/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:4282857/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:1493690/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:3750162/71‑1 (MQ=255)
                                   aTCGCGACCAATAAGCTGCGGACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGa                               <  1:1490101/71‑1 (MQ=255)
                                    tCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATaa                                                           <  1:3838934/42‑1 (MQ=255)
                                    tCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAg                              <  1:1840776/71‑1 (MQ=255)
                                    tCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAg                              >  1:3197500/1‑71 (MQ=255)
                                      gcgACCAATAAGCTGCGGACAAAGGTGATCCTGCaaa                                                              <  1:821014/37‑1 (MQ=255)
                                      gcgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAggg                            >  1:4210872/1‑71 (MQ=255)
                                      gcgACCAATAAGCTGCGGACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAggg                            <  1:1539817/71‑1 (MQ=255)
                                         aCCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTgg                         >  1:3861894/1‑71 (MQ=255)
                                         aCCAATAAGCTGCGGACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTgg                         >  1:3675354/1‑71 (MQ=255)
                                           cAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGCCCATTGAGGGTGGCa                       >  1:2247381/1‑71 (MQ=255)
                                             aTAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGTAAAGCTc                                          <  1:3744461/50‑1 (MQ=255)
                                             aTAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTc                                          <  1:1630165/50‑1 (MQ=255)
                                              tAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCa                                                    >  1:1372383/1‑39 (MQ=255)
                                              tAAGCTGCGGACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAgg                             <  1:4355284/62‑1 (MQ=255)
                                               aaGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATccc                   >  1:4134207/1‑71 (MQ=255)
                                                aGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATcccc                  <  1:1745721/71‑1 (MQ=255)
                                                     cGGACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGc             >  1:2342938/1‑71 (MQ=255)
                                                       gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCaa                >  1:3560617/1‑66 (MQ=255)
                                                       gACAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGccc           >  1:4260232/1‑71 (MQ=255)
                                                        aCAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTc                                          >  1:532697/1‑39 (MQ=255)
                                                         cAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCg          >  1:1346467/1‑70 (MQ=255)
                                                         cAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCg          >  1:1690528/1‑70 (MQ=255)
                                                         cAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGt         <  1:271286/71‑1 (MQ=255)
                                                         cAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGt         <  1:2551384/71‑1 (MQ=255)
                                                         cAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGc                        >  1:2449932/1‑56 (MQ=255)
                                                         cAAAGGTGATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGt         <  1:333388/71‑1 (MQ=255)
                                                          aaaGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGt         >  1:3202752/1‑70 (MQ=255)
                                                           aaGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTaa       >  1:721181/1‑71 (MQ=255)
                                                              gTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTcc                                    <  1:251406/39‑1 (MQ=255)
                                                                gATCATGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCg  >  1:4606100/1‑71 (MQ=255)
                                                                    |                                                                  
TAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCG  >  W3110S.gb/1656644‑1656778

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