Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,143,685 G→A 21.5% intergenic (‑303/‑45) metY ← / → argG tRNA‑Met/argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,143,6850GA21.5% 208.4 / 41.4 143intergenic (‑303/‑45)metY/argGtRNA‑Met/argininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (35/77);  new base A (16/15);  total (51/92)
Fisher's exact test for biased strand distribution p-value = 5.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 143 total aligned reads displayed.

TGTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATG  >  minE/2143622‑2143732
                                                               |                                               
tGTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGa                                               <  1:1041597/66‑1 (MQ=255)
tGTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGa                                               <  1:3081016/66‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGa                                               <  1:1649173/65‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGa                                               <  1:2691666/65‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAAt                                            <  1:481905/68‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAAt                                            <  1:156341/68‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAAt                                            <  1:2576889/68‑1 (MQ=255)
 gTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAAt                                            <  1:1215410/68‑1 (MQ=255)
  tctcCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATcc                                          >  1:2343353/1‑69 (MQ=255)
  tctcCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATcc                                          >  1:1072136/1‑69 (MQ=255)
  tctcCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATcc                                          >  1:1940306/1‑69 (MQ=255)
  tctcCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATcc                                          >  1:397498/1‑69 (MQ=255)
  tctcCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATcc                                          >  1:2710865/1‑69 (MQ=255)
      cGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAgg                                       >  1:2665950/1‑68 (MQ=255)
      cGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAg                                        >  1:2678721/1‑67 (MQ=255)
      cGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCAAgg                                       >  1:3392935/1‑68 (MQ=255)
       atcaaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCt                                     <  1:1041355/66‑1 (MQ=255)
       atcaaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCt                                     <  1:548843/66‑1 (MQ=255)
       atcaaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCt                                     <  1:2055812/66‑1 (MQ=255)
       atcaaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCt                                     <  1:3284269/66‑1 (MQ=255)
        aTAAGAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    <  1:43162/69‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAgg                                       <  1:189336/66‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAgg                                       <  1:330860/66‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAgg                                       <  1:2964960/66‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAgg                                       <  1:2250453/66‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    <  1:1252001/69‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    >  1:238453/1‑69 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    <  1:1183568/69‑1 (MQ=255)
        aTAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    <  1:941510/69‑1 (MQ=255)
          aaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  >  1:2162432/1‑69 (MQ=255)
          aaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  >  1:477850/1‑69 (MQ=255)
          aaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  <  1:477213/69‑1 (MQ=255)
          aaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  >  1:3431054/1‑69 (MQ=255)
          aaaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  >  1:861865/1‑69 (MQ=255)
           aaaGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcat                                  <  1:2093247/68‑1 (MQ=255)
             aGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcatca                                <  1:1521488/68‑1 (MQ=255)
             aGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtcatca                                <  1:2408293/68‑1 (MQ=255)
                  atTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGtt                          >  1:3023897/1‑69 (MQ=255)
                            aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCa                 <  1:467313/68‑1 (MQ=255)
                            aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCa                 <  1:59748/68‑1 (MQ=255)
                            aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCa                 <  1:3138077/68‑1 (MQ=255)
                            aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCa                 <  1:1748133/68‑1 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:415701/69‑1 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              >  1:2655465/1‑69 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              >  1:3400762/1‑69 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:2502856/69‑1 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:488127/69‑1 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:2480852/69‑1 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              >  1:1336480/1‑69 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              >  1:3338513/1‑69 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              >  1:2126545/1‑69 (MQ=255)
                              cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:64997/69‑1 (MQ=255)
                                cATTTATTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGgtgt           <  1:1872270/69‑1 (MQ=255)
                                 atttatttTGCATAAATATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:229067/66‑1 (MQ=255)
                                   ttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGgtgt           <  1:696822/67‑1 (MQ=255)
                                   ttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGgtgt           <  1:1227934/67‑1 (MQ=255)
                                    tatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGgtgt           <  1:2104773/66‑1 (MQ=255)
                                     atttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAt        >  1:2257504/1‑68 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1989189/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1669844/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1645211/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1602484/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1551278/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1423214/68‑1 (MQ=255)
                                       tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:349670/68‑1 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtca                                   >  1:2961848/1‑38 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtca                                   >  1:1460436/1‑38 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtca                                   >  1:1036377/1‑38 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCtc                                    >  1:2309992/1‑37 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:989516/67‑1 (MQ=255)
                                        ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1189217/67‑1 (MQ=255)
                                          gCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:2332357/68‑1 (MQ=255)
                                           cATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:498418/67‑1 (MQ=255)
                                           cATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATg  <  1:85844/68‑1 (MQ=255)
                                           cATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATg  <  1:598753/68‑1 (MQ=255)
                                           cATAAAAATTCAGTGAGAGCGGAAATCCAGGCGCATCATCAGTTAATTAAGCAGGGTGTTATTTTATg  <  1:70875/68‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:360124/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:978267/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:972012/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:964860/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:933761/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:869216/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:1667350/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:809168/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:1904523/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:2302105/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   <  1:545180/64‑1 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:2487628/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:2499011/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:2538201/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:2668143/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:3169472/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:3198246/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTAt   >  1:3434647/1‑64 (MQ=255)
                                              aaaaaTTCAGTGAAAGCAGAAATCCAGGCTCATCATCAGTTAATTAAGAAGGGTGTTATTTTAt   >  1:3059173/1‑64 (MQ=255)
                                              aaaaaTTCAGAGAAAGCAGAAATCCAGGCTCATCATCAGttaattaa                    >  1:560597/1‑47 (MQ=37)
                                               aaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGtt                          <  1:3145104/40‑1 (MQ=255)
                                               aaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGtt                          <  1:402657/40‑1 (MQ=255)
                                                         gagaGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAggg              <  1:1079327/42‑1 (MQ=255)
                                                           gagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTAttt      <  1:1697/48‑1 (MQ=255)
                                                               |                                               
TGTCTCCGATAAAAGATGATTAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATG  >  minE/2143622‑2143732

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: