Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,639,281 A→G 8.7% intergenic (+51/+59) ypfN → / ← ypfH hypothetical protein/predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,639,2810AG 8.7% 251.0 / 11.3 92intergenic (+51/+59)ypfN/ypfHhypothetical protein/predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base G (2/6);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

AGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTC  >  Exported/1639211‑1639351
                                                                      |                                                                       
aGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCaa                                                                         <  1:1008489/71‑1 (MQ=255)
aGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCaa                                                                         <  1:1590836/71‑1 (MQ=255)
   aTCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATgg                                                                      >  1:4104017/1‑71 (MQ=255)
      aaCCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTg                                                                    <  1:1121924/70‑1 (MQ=255)
        ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                 <  1:3704294/71‑1 (MQ=255)
        ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                 >  1:4502322/1‑71 (MQ=255)
        ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                 <  1:3442558/71‑1 (MQ=255)
         cGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                 <  1:773098/70‑1 (MQ=255)
          gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGa                                                               >  1:511664/1‑71 (MQ=255)
             gaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGc                                                           <  1:2470358/71‑1 (MQ=255)
                gTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                         >  1:841144/1‑70 (MQ=255)
                  aaGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGAt                                                      >  1:4334209/1‑71 (MQ=255)
                    gCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGc                                                    <  1:2063037/71‑1 (MQ=255)
                      tGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGa                                                  >  1:217416/1‑71 (MQ=255)
                        aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGa                                                  <  1:1472496/69‑1 (MQ=255)
                        aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACg                                                <  1:170100/71‑1 (MQ=255)
                        aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACg                                                <  1:176796/71‑1 (MQ=255)
                         gTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGc                                               <  1:754492/71‑1 (MQ=255)
                          tGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCt                                              <  1:3028274/71‑1 (MQ=255)
                               attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          >  1:1694301‑M1/4‑46 (MQ=255)
                               attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          <  1:4143815‑M1/68‑26 (MQ=255)
                               attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          <  1:2626015‑M1/68‑26 (MQ=255)
                               attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          >  1:2786583‑M1/4‑46 (MQ=255)
                               actcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          <  1:4576253‑M1/68‑26 (MQ=255)
                               actcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgcta                                                 >  1:800526‑M1/4‑46 (MQ=255)
                               actcgtaGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                          <  1:3539727‑M1/68‑26 (MQ=255)
                               actcgtaGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                          <  1:3142512‑M1/68‑26 (MQ=255)
                                ctcgtAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                          >  1:766880‑M1/3‑45 (MQ=255)
                                ctcgtAGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatc                                         <  1:3266569‑M1/69‑27 (MQ=255)
                                 tcgTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatg                                                             >  1:3737166‑M1/2‑44 (MQ=255)
                                 tcgTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaa                                                                >  1:417648‑M1/2‑44 (MQ=255)
                                  cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:3560763‑M1/1‑43 (MQ=255)
                                  cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                 <  1:2398275‑M1/45‑3 (MQ=255)
                                   gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcagg                                      >  1:3620561‑M1/1‑42 (MQ=255)
                                   gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcag                                       <  1:2527874‑M1/70‑29 (MQ=255)
                                   gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcag                                       <  1:4664754‑M1/70‑29 (MQ=255)
                                   gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                        <  1:2292632‑M1/69‑28 (MQ=255)
                                   gTAGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatca                                        <  1:2061615‑M1/69‑28 (MQ=255)
                                     aGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                   >  1:3526884/1‑71 (MQ=255)
                                     aGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                          <  1:2228200‑M1/65‑26 (MQ=255)
                                      ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                  <  1:3121612/71‑1 (MQ=255)
                                      ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                  <  1:1179540/71‑1 (MQ=255)
                                       gCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatcaggccta                                  <  1:2665404‑M1/71‑34 (MQ=255)
                                         cGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatcaggccta                                  >  1:3745869‑M1/1‑36 (MQ=255)
                                          ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggccta                                  >  1:1051323‑M1/1‑35 (MQ=255)
                                          ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                        >  1:2540654‑M1/1‑35 (MQ=255)
                                          ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                        >  1:1768925‑M1/1‑35 (MQ=255)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                              <  1:415688‑M1/71‑38 (MQ=255)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                              <  1:3744905‑M1/71‑38 (MQ=255)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                 <  1:860661‑M1/36‑3 (MQ=255)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                 <  1:2947988‑M1/36‑3 (MQ=255)
                                           gATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                 <  1:2174536/36‑1 (MQ=255)
                                            aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                          <  1:1256639‑M1/42‑10 (MQ=255)
                                            aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                          <  1:3769719‑M1/42‑10 (MQ=255)
                                            aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                          <  1:1726267‑M1/42‑10 (MQ=255)
                                            aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                            <  1:4210460/71‑1 (MQ=255)
                                            aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                            <  1:4354912/71‑1 (MQ=255)
                                            aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                            <  1:1902689/71‑1 (MQ=255)
                                             tAAGGCGTTTACGCCGCATCCGGCAACAGTGCCGACTGCCTGATGCGACGCTCGCGCGTCTTATCAGGccg                           >  1:718344/1‑70 (MQ=1)
                                              aaGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctga                                                        >  1:60934‑M1/1‑31 (MQ=255)
                                                 gCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTa                       <  1:4425751/71‑1 (MQ=255)
                                                  cGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTaa                      >  1:5151097/1‑71 (MQ=255)
                                                   gTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                              <  1:651551‑M1/63‑38 (MQ=255)
                                                   gTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                         <  1:4731550/35‑1 (MQ=255)
                                                    ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                  <  1:2366829/57‑1 (MQ=255)
                                                    ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                             <  1:3638156/62‑1 (MQ=255)
                                                    ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                    >  1:1261933/1‑71 (MQ=255)
                                                    ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                    >  1:3411478/1‑71 (MQ=255)
                                                     tAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                             <  1:841510/61‑1 (MQ=255)
                                                      aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCa                                >  1:2613281/1‑57 (MQ=255)
                                                      aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                  >  1:931845/1‑71 (MQ=255)
                                                      aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                  >  1:3506308/1‑71 (MQ=255)
                                                      aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                  >  1:1180973/1‑71 (MQ=255)
                                                      aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                  >  1:4503588/1‑71 (MQ=255)
                                                          ccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTa                       >  1:3891078/1‑62 (MQ=255)
                                                           cgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCttt             >  1:5414319/1‑71 (MQ=255)
                                                            gcATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctaca                                >  1:2911165‑M1/1‑17 (MQ=255)
                                                            gcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                   <  1:5360181/48‑1 (MQ=255)
                                                             cATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGt                        >  1:2983973/1‑58 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:3830605‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:4036882‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:2782819‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:4226855‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:3345929‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                           >  1:3260030‑M1/1‑13 (MQ=255)
                                                                ccGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTa           >  1:5385728/1‑68 (MQ=255)
                                                                 cGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                          >  1:1302286‑M1/1‑12 (MQ=255)
                                                                 cGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcat        <  1:4598506/70‑1 (MQ=255)
                                                                 cGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcat        >  1:2368769/1‑70 (MQ=255)
                                                                  ggCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcatca      <  1:3014156/71‑1 (MQ=255)
                                                                    cAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctaca                                <  1:62743‑M1/44‑36 (MQ=255)
                                                                    cAATGGTGAcaaatgcctgatgcgctacgcttatca                                        <  1:3056331‑M1/36‑28 (MQ=255)
                                                                    cAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                   >  1:2206842/1‑40 (MQ=255)
                                                                    cAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTGTACATCAtt    >  1:3445305/1‑71 (MQ=255)
                                                                      aTGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTc  <  1:1408944/71‑1 (MQ=255)
                                                                      |                                                                       
AGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTC  >  Exported/1639211‑1639351

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