Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,727 G→T 15.4% intergenic (+1495/‑284) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7270GT15.4% 388.3 / 33.3 123intergenic (+1495/‑284)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/66);  new base T (12/7);  total (50/73)
Fisher's exact test for biased strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 125 total aligned reads displayed.

GGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCA  >  minE/2572657‑2572795
                                                                      |                                                                    
gggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAg                                                                      >  1:1632056/1‑71 (MQ=255)
    tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                                  <  1:4001960/71‑1 (MQ=255)
    tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                                  <  1:3901341/71‑1 (MQ=255)
    tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                                  <  1:734034/71‑1 (MQ=255)
     cGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                                   >  1:3383275/1‑69 (MQ=255)
      gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                                                >  1:1264899/1‑71 (MQ=255)
       aTTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                               <  1:3702400/71‑1 (MQ=255)
         tCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                                              <  1:605181/70‑1 (MQ=255)
          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                             <  1:3787628/70‑1 (MQ=255)
          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                             <  1:779791/70‑1 (MQ=255)
          ccAGAAATCGGGGTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCCGAGCCTGAgtgt                                                             <  1:3865189/70‑1 (MQ=255)
           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           <  1:567787/71‑1 (MQ=255)
           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           <  1:390349/71‑1 (MQ=255)
            aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                                          <  1:914487/71‑1 (MQ=255)
               aaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                             <  1:1237055/65‑1 (MQ=255)
                 tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           >  1:1477177/1‑65 (MQ=255)
                  cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                                           >  1:1617017/1‑64 (MQ=255)
                   gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                                   >  1:3634672/1‑71 (MQ=255)
                   ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                                   >  1:3595513/4‑71 (MQ=255)
                      ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                   <  1:1275520/52‑1 (MQ=255)
                      ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                             <  1:3534905/58‑1 (MQ=255)
                        cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                                               <  1:3681761/70‑1 (MQ=255)
                        cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                                >  1:1241624/1‑67 (MQ=255)
                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                              >  1:2362725/1‑71 (MQ=255)
                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                              >  1:3488594/1‑71 (MQ=255)
                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                              >  1:3689940/1‑71 (MQ=255)
                         aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                                              <  1:3629112/70‑2 (MQ=255)
                         aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                                              <  1:3263023/70‑2 (MQ=255)
                          aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGtt                                                            <  1:2544727/55‑1 (MQ=255)
                           cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                           <  1:1144066/71‑1 (MQ=255)
                           cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                           >  1:1169944/1‑71 (MQ=255)
                            ggATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                          >  1:1446459/1‑71 (MQ=255)
                              aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGATCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                        <  1:2890548/71‑1 (MQ=255)
                              aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                                              <  1:3245108/49‑1 (MQ=255)
                              aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           >  1:2921150/1‑52 (MQ=255)
                               tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           >  1:86929/1‑51 (MQ=255)
                               tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                        >  1:1461899/1‑70 (MQ=255)
                                 cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                                        >  1:2670562/1‑52 (MQ=255)
                                 cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                                     <  1:3380212/71‑1 (MQ=255)
                                 cTGATAAAAATCCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                                        >  1:2403125/1‑52 (MQ=255)
                                  tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           <  1:2543605/48‑1 (MQ=255)
                                  tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                        <  1:3353523/67‑1 (MQ=255)
                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgtgat                                                            >  1:3220834/1‑44 (MQ=255)
                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                        >  1:346663/1‑50 (MQ=255)
                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                        >  1:2645672/1‑50 (MQ=255)
                                   gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                                            <  1:3902469/46‑3 (MQ=255)
                                   gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                        <  1:504292/66‑1 (MQ=255)
                                    aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                                            >  1:4064961/1‑43 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGATTGTTATACCGCCCCTGCTGCGGTTAATGtt                                  <  1:2329000/71‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                           <  1:2963637/46‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                             <  1:3080931/60‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                             <  1:156782/60‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                         <  1:3131141/64‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                         <  1:2890094/64‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                  <  1:3621980/71‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                  <  1:1950267/71‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                                   <  1:3625784/70‑1 (MQ=255)
                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                                   <  1:2755390/70‑1 (MQ=255)
                                     tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                                 >  1:2141774/1‑55 (MQ=255)
                                      aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                                                >  1:3628962/1‑39 (MQ=255)
                                      aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                         >  1:763830/1‑62 (MQ=255)
                                       aaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                           >  1:1564126/1‑43 (MQ=255)
                                       aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                             <  1:610657/57‑1 (MQ=255)
                                         aaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                             >  1:2694211/1‑71 (MQ=255)
                                           ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                           <  1:1348121/71‑1 (MQ=255)
                                           ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat                            >  1:3399694/1‑70 (MQ=255)
                                             ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                          >  1:429961/1‑70 (MQ=255)
                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                        <  1:612605/71‑1 (MQ=255)
                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                         <  1:1454696/70‑1 (MQ=255)
                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                         <  1:686262/70‑1 (MQ=255)
                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                         <  1:691148/70‑1 (MQ=255)
                                                cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                      <  1:723006/71‑1 (MQ=255)
                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                                     <  1:1095040/55‑1 (MQ=255)
                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                      >  1:1540990/1‑70 (MQ=255)
                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                      >  1:1536064/1‑70 (MQ=255)
                                                  tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGATAATCTTACCGATAATAGTta                    >  1:1392362/1‑71 (MQ=255)
                                                   acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                    >  1:2032962/1‑70 (MQ=255)
                                                   acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                    >  1:2413933/1‑70 (MQ=255)
                                                    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                    <  1:1166089/69‑1 (MQ=255)
                                                    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                    <  1:1552513/69‑1 (MQ=255)
                                                     ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                      >  1:2355988/1‑50 (MQ=255)
                                                     ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                      >  1:66270/1‑50 (MQ=255)
                                                      ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                                     >  1:8035/1‑50 (MQ=255)
                                                      ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                                     >  1:3446/1‑50 (MQ=255)
                                                        ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                             >  1:2983549/1‑40 (MQ=255)
                                                            cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                           <  1:1522713/54‑1 (MQ=255)
                                                            cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                           <  1:2540480/54‑1 (MQ=255)
                                                            cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGcc          <  1:3523044/71‑1 (MQ=255)
                                                            cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGcc          <  1:1406826/71‑1 (MQ=255)
                                                              cAATAGGAGCCTGAGTGTTATACCGCCCCTGCTGAGGTTAATGTTACCGATAATAGTta                    <  1:3474078/59‑1 (MQ=39)
                                                               caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg       <  1:1951602/70‑1 (MQ=255)
                                                                 tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                      <  1:3220334/38‑1 (MQ=255)
                                                                 tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                           >  1:687402/1‑49 (MQ=255)
                                                                 tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                    <  1:2342542/56‑1 (MQ=255)
                                                                  cGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt    >  1:1954597/1‑71 (MQ=255)
                                                                   ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt    >  1:1606653/1‑70 (MQ=255)
                                                                   ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc   <  1:2309267/71‑1 (MQ=255)
                                                                   ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc   <  1:939085/71‑1 (MQ=255)
                                                                    gAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtca  <  1:2707886/71‑1 (MQ=255)
                                                                    gAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtca  <  1:2847561/71‑1 (MQ=255)
                                                                      |                                                                    
GGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCA  >  minE/2572657‑2572795

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: