Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,218,735 G→T 8.5% W288C (TGG→TGT metA → homoserine transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,218,7350GT 8.5% 176.9 / 2.2 58W288C (TGG→TGTmetAhomoserine transsuccinylase
Reads supporting (aligned to +/- strand):  ref base G (33/20);  new base T (3/2);  total (36/22)
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.97e-01

GCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCT  >  W3110S.gb/4218666‑4218804
                                                                     |                                                                     
gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGc                                                                      >  1:261384/1‑71 (MQ=255)
gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGc                                                                      <  1:241698/71‑1 (MQ=255)
gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGc                                                                      >  1:4217823/1‑71 (MQ=255)
   caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCa                                                                   <  1:2572589/71‑1 (MQ=255)
   caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCa                                                                   <  1:2206360/71‑1 (MQ=255)
     aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAc                                                                 <  1:2283861/71‑1 (MQ=255)
        aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTAt                                                              >  1:1063701/1‑71 (MQ=255)
        aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                              >  1:3379495/1‑71 (MQ=255)
        aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                              <  1:1081655/71‑1 (MQ=255)
        aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                              <  1:3651373/71‑1 (MQ=255)
           cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTAc                                                           >  1:2393546/1‑71 (MQ=255)
             cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCaa                                                                  <  1:3440631/62‑1 (MQ=255)
              aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                        <  1:3234193/71‑1 (MQ=255)
                    cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCaa                                                                  >  1:4161440/1‑55 (MQ=255)
                       cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTc                                                                    <  1:2609131/50‑1 (MQ=255)
                         cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                        >  1:3485904/1‑60 (MQ=255)
                         cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                        >  1:119013/1‑60 (MQ=255)
                          gagcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACg                                            <  1:1420597/71‑1 (MQ=255)
                             cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCa                                         <  1:2618298/71‑1 (MQ=255)
                                    gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                  <  1:3940113/71‑1 (MQ=255)
                                    gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                  >  1:4067610/1‑71 (MQ=255)
                                         gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACgg                             >  1:97013/1‑71 (MQ=255)
                                         gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACgg                             >  1:2450181/1‑71 (MQ=255)
                                          tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACg                                     <  1:1282160/62‑1 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGTCTCAACTAt                                                              >  1:980830/1‑35 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                              >  1:3489381/1‑35 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                        >  1:4424816/1‑41 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:2413148/1‑70 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:2125228/1‑70 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:3614661/1‑70 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:3187272/1‑70 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:1829157/1‑70 (MQ=255)
                                            aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           >  1:1790897/1‑70 (MQ=255)
                                               gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                       <  1:2026043/71‑1 (MQ=255)
                                                tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATc                                               >  1:2017424/1‑46 (MQ=255)
                                                tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                       >  1:3917280/1‑70 (MQ=255)
                                                tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                      <  1:3908763/71‑1 (MQ=255)
                                                   tttACTGTTTACCAACTGGCTCAACTATTACGTCt                                                       <  1:1589024/35‑1 (MQ=255)
                                                   tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCa                                                   <  1:1697136/39‑1 (MQ=255)
                                                       cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc               >  1:3108999/1‑71 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACCAGATc                                               >  1:2781614/1‑36 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAc               <  1:3435516/68‑1 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGc             <  1:1622104/70‑1 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCt            >  1:1309539/1‑71 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAAAGCt            >  1:407455/1‑71 (MQ=255)
                                                          tttACCAACTGGCTCAACTATTACGTCTACAAGATCACGCCATACGATCTACGGCACATGAATCCAAc               <  1:1580752/68‑1 (MQ=255)
                                                           ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg           >  1:3765738/1‑71 (MQ=255)
                                                           ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg           >  1:2433254/1‑71 (MQ=255)
                                                            tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                      >  1:57395/1‑59 (MQ=255)
                                                            tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                      >  1:377629/1‑59 (MQ=255)
                                                            tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg           <  1:2866237/70‑1 (MQ=255)
                                                            tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGATg           >  1:3932016/1‑70 (MQ=255)
                                                              ccAACTGTCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGa         <  1:2910611/70‑1 (MQ=255)
                                                              ccAACTGTCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt        >  1:1710811/1‑71 (MQ=255)
                                                                aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt        >  1:1591318/1‑69 (MQ=255)
                                                                 aCTGGCTCAACTATTACGTCTACCAGATCACGCCATa                                       >  1:311682/1‑37 (MQ=255)
                                                                 aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTa                                >  1:1366768/1‑44 (MQ=255)
                                                                 aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa      >  1:287286/1‑70 (MQ=255)
                                                                     gCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                           <  1:1080957/45‑1 (MQ=255)
                                                                     gCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct  <  1:3008053/70‑1 (MQ=255)
                                                                     |                                                                     
GCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCT  >  W3110S.gb/4218666‑4218804

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: