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% 210.2 / 4.3 71W288C (TGG→TGTmetAhomoserine transsuccinylase
Reads supporting (aligned to +/- strand):  ref base G (38/25);  new base T (4/2);  total (42/29)
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.05e-01

CGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATC  >  W3110S.gb/4218665‑4218803
                                                                      |                                                                    
cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTgg                                                                      >  1:2119309/1‑71 (MQ=255)
cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTgg                                                                      >  1:4086917/1‑71 (MQ=255)
 gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGc                                                                     <  1:322322/71‑1 (MQ=255)
  cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCt                                                                    >  1:1286394/1‑71 (MQ=255)
  cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCt                                                                    >  1:2896532/1‑71 (MQ=255)
   acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGCCTc                                                                   <  1:3025247/71‑1 (MQ=255)
    caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTc                                                                   >  1:2983589/1‑70 (MQ=255)
      aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAc                                                                >  1:2759493/1‑71 (MQ=255)
      aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAc                                                                >  1:1754823/1‑71 (MQ=255)
      aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAc                                                                <  1:1776311/71‑1 (MQ=255)
      aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAc                                                                >  1:3983228/1‑71 (MQ=255)
        gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTa                                                              >  1:3741704/1‑71 (MQ=255)
         aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTa                                                              <  1:3414343/70‑1 (MQ=255)
         aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTAt                                                             <  1:1354855/71‑1 (MQ=255)
         aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                             <  1:1518906/71‑1 (MQ=255)
          tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTa                                                              <  1:1199251/69‑1 (MQ=255)
            cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTAc                                                          >  1:2478453/1‑71 (MQ=255)
               aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                       >  1:4261271/1‑71 (MQ=255)
                aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAtt                                                            <  1:473140/65‑1 (MQ=255)
                aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                       >  1:2600070/1‑70 (MQ=255)
                  aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTATTACGTCTAc                                                    >  1:346261/1‑71 (MQ=255)
                  aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAc                                                    >  1:3620831/1‑71 (MQ=255)
                     cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAc                                                    >  1:499119/1‑68 (MQ=255)
                        cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCt                                                                    >  1:2458733/1‑49 (MQ=255)
                        cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACt                                                               >  1:1932804/1‑54 (MQ=255)
                         gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCa                                             >  1:4061157/1‑71 (MQ=255)
                          cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTAc                                                          >  1:2740247/1‑57 (MQ=255)
                          cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCAc                                            >  1:583005/1‑71 (MQ=255)
                          cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCAc                                            >  1:3733751/1‑71 (MQ=255)
                            agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTAt                                                             <  1:3150566/52‑1 (MQ=255)
                             gcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGcc                                         <  1:2941841/71‑1 (MQ=255)
                              cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCa                                                  >  1:4113345/1‑61 (MQ=255)
                              cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCa                                        <  1:1618362/71‑1 (MQ=255)
                               gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCa                                        >  1:2645456/1‑70 (MQ=255)
                                 gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATAc                                     <  1:2990404/71‑1 (MQ=255)
                                 gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATAc                                     <  1:1048512/71‑1 (MQ=255)
                                  cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATAc                                     <  1:3208781/70‑1 (MQ=255)
                                  cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACg                                    <  1:617708/71‑1 (MQ=255)
                                   tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTATTACGTCTACCAGATCACGCCATACGa                                   >  1:1548395/1‑71 (MQ=255)
                                     gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAt                                  >  1:1930864/1‑70 (MQ=255)
                                     gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                 >  1:51939/1‑71 (MQ=255)
                                     gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                 <  1:3409921/71‑1 (MQ=255)
                                        tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAg                                                 <  1:3514290/52‑1 (MQ=255)
                                        tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCa                                        >  1:647062/1‑61 (MQ=255)
                                            cACGGTAATTTACTGTTTACCAACTGTCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                          >  1:212637/1‑71 (MQ=255)
                                            cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAcc                                                   >  1:984707/1‑46 (MQ=255)
                                            cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                          >  1:3451347/1‑71 (MQ=255)
                                             aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTc                                                       >  1:757567/1‑41 (MQ=255)
                                             aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCa                                                  >  1:3837808/1‑46 (MQ=255)
                                               ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACAt                       <  1:2601297/71‑1 (MQ=255)
                                               ggTAATTTACTGTTTACCAACTGGCTCAACTATTACCTCTACCAGATCACGCCATACGATCTACGGCACAt                       <  1:3059187/71‑1 (MQ=255)
                                                gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGAt                                               <  1:715544/46‑1 (MQ=255)
                                                 tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                     <  1:968815/71‑1 (MQ=255)
                                                 tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                     <  1:156223/71‑1 (MQ=255)
                                                 tAATTTACTGTTTACCAACTGCCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                     <  1:3626097/71‑1 (MQ=255)
                                                  aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                    >  1:2249195/1‑71 (MQ=255)
                                                   aTTTACTGTTTACCTACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                    >  1:1928773/1‑70 (MQ=255)
                                                   aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                   >  1:1170580/1‑71 (MQ=255)
                                                     ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAg                                                 <  1:2913705/39‑1 (MQ=255)
                                                      tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                 >  1:3844216/1‑70 (MQ=255)
                                                      tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                 <  1:3984597/70‑1 (MQ=255)
                                                       aCTGTTTACCAACTGTCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCaa               >  1:2492745/1‑71 (MQ=255)
                                                           tttACCAACTGGCTCAACTATTACGTCTACCAGAt                                               >  1:244533/1‑35 (MQ=255)
                                                           tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                    >  1:3337046/1‑62 (MQ=255)
                                                           tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg             <  1:2013130/69‑1 (MQ=255)
                                                           tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg             <  1:344473/69‑1 (MQ=255)
                                                             tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg          >  1:1000152/1‑70 (MQ=255)
                                                             tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg          <  1:646761/70‑1 (MQ=255)
                                                             tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg          <  1:1719704/70‑1 (MQ=255)
                                                             tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg          >  1:3562671/1‑70 (MQ=255)
                                                               ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                  >  1:1447841/1‑60 (MQ=255)
                                                                  aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa    >  1:2177339/1‑71 (MQ=255)
                                                                  aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa    >  1:801489/1‑71 (MQ=255)
                                                                     ggCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGa                     >  1:2489201/1‑51 (MQ=255)
                                                                      gCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc  <  1:1566849/69‑1 (MQ=255)
                                                                      |                                                                    
CGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATC  >  W3110S.gb/4218665‑4218803

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: