Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,477,381 A→G 9.2% E48G (GAA→GGA)  yfaH → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,477,3810AG 9.2% 210.3 / 12.1 76E48G (GAA→GGA) yfaHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base G (4/3);  total (41/35)
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.08e-01

TTGTGTTCCGTTCCTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAATTTGATATGA  >  Exported/1477311‑1477450
                                                                      |                                                                     
ttGTGTTCCGTTCCTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGa                                                                       <  1:3433701/71‑1 (MQ=255)
  gtgtTCCGTTCCTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGaa                                                                      <  1:1381983/70‑1 (MQ=255)
     ttccgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAgg                                                                    >  1:3679546/1‑69 (MQ=255)
      tccgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAg                                                                 >  1:389414/1‑71 (MQ=255)
      tccgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAg                                                                 >  1:337949/1‑71 (MQ=255)
        cgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTc                                                               >  1:729827/1‑71 (MQ=255)
        cgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTc                                                               >  1:790560/1‑71 (MQ=255)
        cgttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTc                                                               >  1:2566417/1‑71 (MQ=255)
         gttccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGGAGGCAGt                                                                >  1:1973578/1‑69 (MQ=255)
           tccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCAc                                                            >  1:1825708/1‑71 (MQ=255)
           tccTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCAc                                                            >  1:3257519/1‑71 (MQ=255)
             ctctCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTg                                                          >  1:968041/1‑71 (MQ=255)
             ctctCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTg                                                          >  1:2055919/1‑71 (MQ=255)
                tcGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGt                                                         <  1:1391161/69‑1 (MQ=255)
                 cGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      >  1:3023670/1‑71 (MQ=255)
                  gAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      >  1:2222871/1‑70 (MQ=255)
                  gAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtt                                                        <  1:1697952/68‑1 (MQ=255)
                      gaaCTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGt                                                         <  1:2469717/62‑1 (MQ=255)
                      cAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      <  1:482106/66‑1 (MQ=255)
                       aaCTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGaa                                                                      <  1:1308659/49‑1 (MQ=255)
                       aaCTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTtactac                                                <  1:1426976/71‑1 (MQ=255)
                        aCTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACa                                               <  1:619827/71‑1 (MQ=255)
                        aCTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACa                                               <  1:3321110/71‑1 (MQ=255)
                         cTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtt                                                        <  1:2505479/61‑1 (MQ=255)
                          tttCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACaaa                                             >  1:2553558/1‑71 (MQ=255)
                          tttCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACaaa                                             >  1:1829472/1‑71 (MQ=255)
                          tttCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACaaa                                             >  1:3840183/1‑71 (MQ=255)
                            tCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGGAGGCAGTCCACTGTTTTTACTTCaaa                                             <  1:68300/69‑1 (MQ=255)
                                tCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAAtgtg                                         <  1:71993/69‑1 (MQ=255)
                                tCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAAtgtg                                         <  1:67498/69‑1 (MQ=255)
                                 cGACATATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCAt                                      >  1:569340/1‑71 (MQ=255)
                                 cGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCAt                                      >  1:1327711/1‑71 (MQ=255)
                                 cGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCAt                                      >  1:2816028/1‑71 (MQ=255)
                                   aCAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      <  1:2078386/53‑1 (MQ=255)
                                   aCAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      <  1:1306377/53‑1 (MQ=255)
                                   aCAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      <  1:1178987/53‑1 (MQ=255)
                                    cAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGGAGGCAGTCCACTGTTTTTACTACAAATGTGCATgg                                    >  1:3512810/1‑70 (MQ=255)
                                       aTTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATgg                                    >  1:1696064/1‑67 (MQ=255)
                                       aTTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGa                                   >  1:3751674/1‑68 (MQ=255)
                                        tttCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtt                                                        <  1:1781267/46‑1 (MQ=255)
                                         ttCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtttt                                                      >  1:2280433/1‑47 (MQ=255)
                                          tCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTc                                                               >  1:718322/1‑37 (MQ=255)
                                          tCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGtt                                                        >  1:3452550/1‑44 (MQ=255)
                                          tCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAAt                              <  1:1137569/70‑1 (MQ=255)
                                          tCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAAt                              >  1:990131/1‑70 (MQ=255)
                                            gTTGCTGGCTAAGGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAg                           >  1:163391/1‑71 (MQ=255)
                                             ttGCTGGCTAATGAAAAGCCGGTAGAAGGCAGTCCACTGttttt                                                     <  1:2919204/44‑1 (MQ=255)
                                             ttGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAAc                                 <  1:2590640/64‑1 (MQ=255)
                                             ttGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAAc                                 <  1:2301706/64‑1 (MQ=255)
                                              tGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAAc                                 <  1:238515/63‑1 (MQ=255)
                                                 tGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGt                          <  1:1079315/67‑1 (MQ=255)
                                                 tGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGt                          <  1:2199667/67‑1 (MQ=255)
                                                 tGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGt                          <  1:2949005/67‑1 (MQ=255)
                                                   gCTAAAGAAAAGCCGGTAGGAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGAt                    <  1:3317802/71‑1 (MQ=255)
                                                   gCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGa                     <  1:1903526/70‑1 (MQ=255)
                                                   gCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGa                     <  1:1285470/70‑1 (MQ=255)
                                                   gCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGAt                    >  1:3512113/1‑71 (MQ=255)
                                                    cTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAAtt                             <  1:2304052/61‑1 (MQ=255)
                                                    cTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAAtt                             <  1:2580417/61‑1 (MQ=255)
                                                    cTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAAtt                             <  1:2575846/61‑1 (MQ=255)
                                                     tAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGAtt                   >  1:1364438/1‑70 (MQ=255)
                                                     tAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTg                  >  1:3235628/1‑71 (MQ=255)
                                                     tAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTATTTACTACAAATGTGCATGGAACAATTAGTGGCGATTg                  >  1:1097500/1‑71 (MQ=255)
                                                      aaaGAAAAGCCGGTAGGAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTg                  >  1:183966/1‑70 (MQ=255)
                                                      aaaGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTg                  >  1:2532432/1‑70 (MQ=255)
                                                       aaGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTgg                 >  1:3430131/1‑70 (MQ=255)
                                                       aaGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTgg                 >  1:2453641/1‑70 (MQ=255)
                                                       aaGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTgg                 >  1:738760/1‑70 (MQ=255)
                                                       aaGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTgg                 >  1:1237709/1‑70 (MQ=255)
                                                          aaaaGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGAt                    <  1:207406/64‑1 (MQ=255)
                                                             aGCCGGTAGGAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGaaaa             <  1:2453825/68‑1 (MQ=255)
                                                              gCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAAttt         <  1:3785558/71‑1 (MQ=255)
                                                              gCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAAttt         <  1:2972804/71‑1 (MQ=255)
                                                              gCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAAtt          >  1:1291707/1‑70 (MQ=255)
                                                               ccGTTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAAttt         >  1:1542802/1‑70 (MQ=255)
                                                               ccGGTAGGAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAAt           >  1:374875/1‑68 (MQ=255)
                                                                cGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAATTTGa       <  1:1260147/71‑1 (MQ=255)
                                                                    aGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGAtt                   >  1:317238/1‑55 (MQ=255)
                                                                     gAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAATTTGATATga  <  1:2107833/71‑1 (MQ=255)
                                                                      |                                                                     
TTGTGTTCCGTTCCTCTCGAACCAACTTTCTATCGACAGATTTCGTTGCTGGCTAAAGAAAAGCCGGTAGAAGGCAGTCCACTGTTTTTACTACAAATGTGCATGGAACAATTAGTGGCGATTGGAAAAATTTGATATGA  >  Exported/1477311‑1477450

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: