Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,706,900 T→C 9.9% I159V (ATT→GTT)  yphB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,706,9000TC 9.9% 267.7 / 17.3 91I159V (ATT→GTT) yphBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (51/31);  new base C (5/4);  total (56/35)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTT  >  Exported/1706834‑1706968
                                                                  |                                                                    
gctCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCcgc                                                                  >  1:215346/1‑71 (MQ=255)
 ctCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCcgc                                                                  <  1:2704092/70‑1 (MQ=255)
  tCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAACCCGCGt                                                                >  1:199702/1‑71 (MQ=255)
  tCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAACCCGCGt                                                                >  1:2868304/1‑71 (MQ=255)
    gCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAACCCGCGtt                                                               <  1:4121172/70‑1 (MQ=255)
        aaCTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:1275226/1‑71 (MQ=255)
        aaCTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:1678297/1‑71 (MQ=255)
        aaCTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:4486607/1‑71 (MQ=255)
             aCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGAc                                                      >  1:3032225/1‑70 (MQ=255)
                 gccagccaCTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGGGCGACAACgg                                                 >  1:554484/1‑70 (MQ=255)
                 gccagccaCTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAACCCGCGTTTGCGGCGACAACg                                                  >  1:1660497/1‑70 (MQ=255)
                    agccaCTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaa                                               >  1:2829318/1‑70 (MQ=255)
                    agccaCTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaa                                               >  1:4608772/1‑70 (MQ=255)
                        gcTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaaaa                                             >  1:2968423/2‑68 (MQ=255)
                        aCTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGa                                                >  1:4630636/1‑65 (MQ=255)
                          tGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGtt                                                               <  1:661234/48‑1 (MQ=255)
                          tGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGtt                                                               <  1:1897492/48‑1 (MQ=255)
                          tGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGtt                                                               <  1:510840/48‑1 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaaaa                                             >  1:807323/1‑65 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAgg                                        >  1:1563595/1‑70 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGa                                       >  1:3169644/1‑71 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGa                                       >  1:2412551/1‑71 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGa                                       >  1:2590193/1‑71 (MQ=255)
                           gCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAACCCGCGTTTGCGGCGACAACGGAAAATAAGGa                                       >  1:2593519/1‑71 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATc                                                                     >  1:3057353/1‑37 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATc                                                                     >  1:2054887/1‑37 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATc                                                                     >  1:2245645/1‑37 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATc                                                                     >  1:2413661/1‑37 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATc                                                                     >  1:2507033/1‑37 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGt                                                                >  1:2352690/1‑42 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAAt                                            >  1:3533829/1‑62 (MQ=255)
                               ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAAt                                            >  1:4703444/1‑62 (MQ=255)
                                cGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGc                                    >  1:1560706/1‑69 (MQ=255)
                                 gCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAg                                 >  1:3215483/1‑71 (MQ=255)
                                  cTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaaaa                                             <  1:988864/58‑1 (MQ=255)
                                  cTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGc                                <  1:2254585/71‑1 (MQ=255)
                                  cTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGc                                <  1:1457712/71‑1 (MQ=255)
                                    ccagccagTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCg                              >  1:1406212/1‑71 (MQ=255)
                                     cagccagTAACCGCTCGCCTGTGCCTGAATCcgcg                                                                 <  1:854962/35‑1 (MQ=255)
                                     cagccagTAACCGCTCGCCTGTGCCTGAATCcgcg                                                                 <  1:2709315/35‑1 (MQ=255)
                                     cagccagTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGcc                                   >  1:4274153/1‑65 (MQ=255)
                                     cagccagTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGcc                               <  1:2980717/69‑1 (MQ=255)
                                       gccagTAACCGCTCGCCTGTGCCTGAATCCGGGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTa                           >  1:3768907/1‑71 (MQ=255)
                                       gccagTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTa                           >  1:982249/1‑71 (MQ=255)
                                          agTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTAcc                         >  1:1344322/1‑70 (MQ=255)
                                          agTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTAcc                         >  1:260000/1‑70 (MQ=255)
                                           gTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaaaa                                             <  1:248166/49‑1 (MQ=255)
                                           gTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAAt                                            <  1:3984798/50‑1 (MQ=255)
                                           gTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCaa                       <  1:3780981/71‑1 (MQ=255)
                                            tAACCGCTCGCCTGTGCCTTAATCCGCGTTTGCCGCGACAACGGAAAATAAGGATGCCAGCCGGTACCaaa                      >  1:738254/1‑71 (MQ=255)
                                            tAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:507738/1‑35 (MQ=255)
                                            tAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:2412214/1‑35 (MQ=255)
                                            tAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg                                                          >  1:2033558/1‑35 (MQ=255)
                                            tAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTAc                          <  1:1093304/67‑1 (MQ=255)
                                              aCCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:956958/71‑1 (MQ=255)
                                              aCCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:3654685/71‑1 (MQ=255)
                                              aCCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:770478/71‑1 (MQ=255)
                                              aCCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:2559962/71‑1 (MQ=255)
                                               ccGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATgg                   <  1:1869404/71‑1 (MQ=255)
                                                cGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGc                  <  1:2036074/71‑1 (MQ=255)
                                                 gCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATgg                   >  1:3631662/1‑69 (MQ=255)
                                                      ccTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtc            >  1:3581153/1‑71 (MQ=255)
                                                       cTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtct           >  1:561134/1‑71 (MQ=255)
                                                       cTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtct           >  1:4757015/1‑71 (MQ=255)
                                                       cTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtct           >  1:4627965/1‑71 (MQ=255)
                                                       cTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtct           >  1:3396580/1‑71 (MQ=255)
                                                       cTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtc            >  1:1157652/1‑70 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGcc                               <  1:194277/50‑1 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGcc                               <  1:2381412/50‑1 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCgg                             >  1:3200658/1‑52 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCaaa                      >  1:3596549/1‑59 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCaaa                      >  1:2295792/1‑59 (MQ=255)
                                                        tgtgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtctc          >  1:3599168/1‑71 (MQ=255)
                                                          tgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTg        <  1:4243830/71‑1 (MQ=255)
                                                          tgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTg        <  1:1843206/71‑1 (MQ=255)
                                                          tgCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTg        <  1:2041965/71‑1 (MQ=255)
                                                          tgCCTGAACCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtctc          <  1:624312/69‑1 (MQ=255)
                                                          tgCCTGAACCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtctc          <  1:894824/69‑1 (MQ=255)
                                                          tgCCTGAACCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtctc          <  1:3619560/69‑1 (MQ=255)
                                                            ccTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGc               >  1:4553127/1‑62 (MQ=255)
                                                            ccTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGc               >  1:492261/1‑62 (MQ=255)
                                                            ccTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGcc      >  1:4661523/1‑71 (MQ=255)
                                                            ccTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGcc      <  1:3465243/71‑1 (MQ=255)
                                                            ccTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGcc      <  1:4378063/71‑1 (MQ=255)
                                                               gAACCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCt   >  1:3925332/1‑71 (MQ=255)
                                                                aaTCCGCGTTTGCGGCGACAACGGAAAATAAGGATGc                                    <  1:3739196/37‑1 (MQ=255)
                                                                aaTCCGCGTTTGCGGCGACAACGGAAAATAAGGATGc                                    <  1:530884/37‑1 (MQ=255)
                                                                aaTCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGc                  >  1:3592657/1‑55 (MQ=255)
                                                                aaTCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCATCGTCTCTGCCCCtt  <  1:189340/71‑1 (MQ=255)
                                                                aaTCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCtt  >  1:1850411/1‑71 (MQ=255)
                                                                aaTCCGAGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCtt  <  1:2339288/71‑1 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCa                                  >  1:842557/1‑37 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCa                                  >  1:933414/1‑37 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCa                                  >  1:2645901/1‑37 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCa                        <  1:2494186/47‑1 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:4354659/51‑1 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:4349426/51‑1 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATg                    <  1:3075662/51‑1 (MQ=255)
                                                                  tCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGtctc          >  1:1713028/1‑61 (MQ=255)
                                                                  |                                                                    
GCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTT  >  Exported/1706834‑1706968

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: