Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,218,901 T→A 8.1% intergenic (+100/‑169) metA → / → aceB homoserine transsuccinylase/malate synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,218,9010TA 8.1% 203.3 / 9.0 74intergenic (+100/‑169)metA/aceBhomoserine transsuccinylase/malate synthase A
Reads supporting (aligned to +/- strand):  ref base T (37/31);  new base A (4/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGATTTTGCA  >  W3110S.gb/4218832‑4218971
                                                                     |                                                                      
cAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAAt                                                                        <  1:3337164/70‑1 (MQ=255)
    aCCTTACCTAAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACATCAGCAGGTGAATAAATTTTa                                                                    <  1:3150940/70‑1 (MQ=255)
     ccTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTAtt                                                                  <  1:2593063/71‑1 (MQ=255)
      cTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTAt                                                                   >  1:5301181/1‑69 (MQ=255)
      cTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAAATTTATTc                                                                 <  1:3838897/71‑1 (MQ=255)
       ttACCTCAGGCACCTTAGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCa                                                                >  1:4410659/1‑71 (MQ=255)
        tACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCa                                                                >  1:5167122/1‑70 (MQ=255)
          ccTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCatat                                                             >  1:5016912/1‑71 (MQ=255)
            tCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTAAGCAGGTGAATAAATTTTATTCATATTg                                                           >  1:7477575/1‑71 (MQ=255)
              aGGCACCTTCTGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGt                                                          >  1:6139177/1‑70 (MQ=255)
              aGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGtt                                                         <  1:7782892/71‑1 (MQ=255)
                  aCCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTAt                                                       >  1:1529616/1‑69 (MQ=255)
                  aCCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATca                                                     <  1:1027430/71‑1 (MQ=255)
                  aCCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATc                                                      >  1:7572859/1‑70 (MQ=255)
                     ttCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCatat                                                             >  1:6326743/1‑60 (MQ=255)
                      tCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac                                                   >  1:549190/1‑69 (MQ=255)
                       cGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAAt                                                                        <  1:2788859/47‑1 (MQ=255)
                       cGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAg                                                >  1:4545714/1‑71 (MQ=255)
                        gggTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTAtt                                                                  >  1:4417992/1‑52 (MQ=255)
                          gTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTa                                             <  1:5395553/71‑1 (MQ=255)
                          gTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTa                                             <  1:5866983/71‑1 (MQ=255)
                               ttttttATTTCCGAAACGTACCTCAGCAGGTGAATAAAATTTATTc                                                                 <  1:7586769/46‑1 (MQ=255)
                                tttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATc                                           >  1:1270283/1‑67 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAAt                                                                        >  1:1072583/1‑37 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCa                                                                >  1:5499340/1‑45 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGtt                                                         >  1:792770/1‑52 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGtt                                                         >  1:2676666/1‑52 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaacaa                                                 >  1:1634966/1‑60 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTa                                      >  1:2353294/1‑71 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTa                                      >  1:6272123/1‑71 (MQ=255)
                                 ttttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTa                                      >  1:2572621/1‑71 (MQ=255)
                                  tttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac                                                   >  1:3934172/1‑57 (MQ=255)
                                  tttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac                                                   >  1:491317/1‑57 (MQ=255)
                                  tttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTa                                      <  1:4198888/70‑1 (MQ=255)
                                  tttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTAt                                     >  1:6268680/1‑71 (MQ=255)
                                   ttATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAACTAtt                                    >  1:6587966/1‑71 (MQ=255)
                                      tttCCGAAATGTACCTCAGCAGGTGAATAAATTTTATTCatat                                                             <  1:4558336/43‑1 (MQ=255)
                                      tttCCGAAACGTACCTCAGCAGGTGAATAAATTTTAtt                                                                  <  1:6295503/38‑1 (MQ=255)
                                      tttCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCatat                                                             <  1:2534668/43‑1 (MQ=255)
                                      tttCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCatat                                                             <  1:4967256/43‑1 (MQ=255)
                                      tttCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTAttttt                                 <  1:1271344/71‑1 (MQ=255)
                                       ttCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaaca                                                  <  1:7333843/53‑1 (MQ=255)
                                         ccGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaat                              >  1:3359873/1‑71 (MQ=255)
                                           gAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAg                                        <  1:1534566/59‑1 (MQ=255)
                                           gAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaatta                            <  1:5735408/71‑1 (MQ=255)
                                           gAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaatta                            <  1:4899292/71‑1 (MQ=255)
                                           gAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTtta                                <  1:7583168/67‑1 (MQ=255)
                                            aaaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaattaa                           <  1:14983/71‑1 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATTAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTaaa                          <  1:1528053/71‑1 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaattaa                           >  1:1690825/1‑70 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTttaattaa                           >  1:2970357/1‑70 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTaaa                          <  1:2022345/71‑1 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTaaa                          <  1:5955320/71‑1 (MQ=255)
                                             aaCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTaaa                          <  1:6236193/71‑1 (MQ=255)
                                              aCGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTaaaa                         >  1:6222613/1‑71 (MQ=255)
                                                 tACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTAttttt                                 <  1:1154592/60‑1 (MQ=255)
                                                    cTCATCAGGTGAATAAATTTTATTCATATTGTTATcaacaa                                                 <  1:6692475/41‑1 (MQ=255)
                                                    cTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaacaa                                                 <  1:1651897/41‑1 (MQ=255)
                                                    cTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaacaa                                                 <  1:2574033/41‑1 (MQ=255)
                                                      cagcagGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTAttttt                                 <  1:5349259/55‑1 (MQ=255)
                                                       agcagGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTg                >  1:4505180/1‑71 (MQ=255)
                                                       agcagGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTg                <  1:4314419/71‑1 (MQ=255)
                                                        gcagGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGt               >  1:6207216/1‑71 (MQ=255)
                                                         cagGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGtt              <  1:2737431/71‑1 (MQ=255)
                                                            gTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCaa                                         >  1:1227308/1‑41 (MQ=255)
                                                            gTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGtttt            <  1:4302248/70‑1 (MQ=255)
                                                             tGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGttttt           >  1:7479713/1‑70 (MQ=255)
                                                             tGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGttttt           >  1:7551297/1‑70 (MQ=255)
                                                             tGAATAAAATTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGttttt           >  1:3290439/1‑70 (MQ=255)
                                                              gAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGa         >  1:1585118/1‑71 (MQ=255)
                                                              gAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGa         >  1:198173/1‑71 (MQ=255)
                                                               aaTAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAAtt                 <  1:7219195/62‑1 (MQ=255)
                                                                aTAAAATTTATTCACATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGAt        >  1:5672689/1‑70 (MQ=255)
                                                                aTAAAATTTATTCACATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGAt        >  1:6977824/1‑70 (MQ=255)
                                                                aTAAAATTTATTCACATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGAt        >  1:2326967/1‑70 (MQ=255)
                                                                 tAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGAttt      >  1:7096267/1‑71 (MQ=255)
                                                                 tAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGAttt      >  1:1761156/1‑71 (MQ=255)
                                                                     ttttATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGt               <  1:5828416/58‑1 (MQ=255)
                                                                     ttttATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGt               <  1:5223379/58‑1 (MQ=255)
                                                                     ttttATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGATTTTGCa  >  1:3923235/1‑71 (MQ=255)
                                                                     |                                                                      
CAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTATCAAGTATTTTTAATTAAAATGGAAATTGTTTTTGATTTTGCA  >  W3110S.gb/4218832‑4218971

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: