Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,717 C→T 7.3% intergenic (+1485/‑294) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7170CT 7.3% 392.2 / 14.3 123intergenic (+1485/‑294)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (50/63);  new base T (6/3);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 133 total aligned reads displayed.

AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCC  >  minE/2572648‑2572787
                                                                     |                                                                      
aaTCCCAAGGTGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                                                                       >  1:1612121/1‑71 (MQ=255)
  tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      >  1:1021027/1‑70 (MQ=255)
  tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      >  1:1269454/1‑70 (MQ=255)
  tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                                     <  1:2772588/71‑1 (MQ=255)
  tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAATCAACCTACCTTCCAGCCCa                                                                     <  1:3674630/71‑1 (MQ=255)
    ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGGCCAAt                                                                   >  1:3284214/1‑71 (MQ=255)
    ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                                    >  1:1063436/1‑70 (MQ=255)
    ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   >  1:2753554/1‑71 (MQ=255)
     cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                                     <  1:933206/68‑1 (MQ=255)
      aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:2666170/1‑71 (MQ=255)
         gggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAg                                                              >  1:1632056/1‑71 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                          <  1:3901341/71‑1 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                          <  1:734034/71‑1 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                          <  1:4001960/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)
                 ttCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      <  1:874905/55‑1 (MQ=255)
                  tCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                                      <  1:605181/70‑1 (MQ=255)
                   ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:3787628/70‑1 (MQ=255)
                   ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:779791/70‑1 (MQ=255)
                    cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   <  1:390349/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   <  1:1872515/54‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                                  <  1:131062/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                                  <  1:914487/71‑1 (MQ=255)
                      gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                 >  1:3445080/1‑71 (MQ=255)
                          tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   >  1:1477177/1‑65 (MQ=255)
                           cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGccc                                                                      >  1:3972475/1‑45 (MQ=255)
                           cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                                   >  1:1617017/1‑64 (MQ=255)
                            gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   >  1:702095/1‑47 (MQ=255)
                            gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                           >  1:3634672/1‑71 (MQ=255)
                              gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:2946995/1‑47 (MQ=255)
                              gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                         <  1:782865/71‑1 (MQ=255)
                               ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:3534905/58‑1 (MQ=255)
                                 cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                        >  1:1241624/1‑67 (MQ=255)
                                 cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                      >  1:2362725/1‑71 (MQ=255)
                                 cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                      >  1:3689940/1‑71 (MQ=255)
                                  aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                                      <  1:3263023/70‑2 (MQ=255)
                                   aCGGATACTGATACATATTCCATCCACCTCTCAACCCAATCgg                                                                >  1:2531630/1‑43 (MQ=255)
                                   aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGtt                                                    <  1:2544727/55‑1 (MQ=255)
                                    cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                   >  1:1169944/1‑71 (MQ=255)
                                     ggATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                  >  1:1446459/1‑71 (MQ=255)
                                     ggATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      >  1:1016167/1‑35 (MQ=255)
                                      gATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCg                                                                 <  1:545833/39‑1 (MQ=25)
                                       aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGATCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                <  1:2890548/71‑1 (MQ=255)
                                       aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   >  1:2921150/1‑52 (MQ=255)
                                       aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                 >  1:1447469/1‑70 (MQ=255)
                                       aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                 >  1:532667/1‑70 (MQ=255)
                                         aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                               <  1:1864114/38‑1 (MQ=12)
                                          cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                                >  1:2670562/1‑52 (MQ=255)
                                          cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                              <  1:364774/70‑1 (MQ=255)
                                          cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                             <  1:3380212/71‑1 (MQ=255)
                                          cTGATAAAAATCCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                                >  1:2403125/1‑52 (MQ=255)
                                           tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                          <  1:328303/41‑1 (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:2645672/1‑50 (MQ=255)
                                            gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                >  1:1525899/1‑50 (MQ=255)
                                            gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                >  1:346663/1‑50 (MQ=255)
                                            gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                                    <  1:3902469/46‑3 (MQ=255)
                                            gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                 >  1:2062969/1‑49 (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: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)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                           <  1:2755390/70‑1 (MQ=255)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                           <  1:3625784/70‑1 (MQ=255)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                          <  1:3137623/70‑1 (MQ=255)
                                              tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                         >  1:2141774/1‑55 (MQ=255)
                                               aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                                        >  1:3628962/1‑39 (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)
                                                    ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                   <  1:351529/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                <  1:1199011/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                <  1:612605/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                 <  1:1454696/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                 <  1:691148/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                 <  1:686262/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)
                                                           tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGATAATCTTACCGATAATAGTta            >  1:1392362/1‑71 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            >  1:2413933/1‑70 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            >  1:2032962/1‑70 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            <  1:1552513/69‑1 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            <  1:1166089/69‑1 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                              >  1:66270/1‑50 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                              >  1:2355988/1‑50 (MQ=255)
                                                               ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                             >  1:3446/1‑50 (MQ=255)
                                                               ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                             >  1:8035/1‑50 (MQ=255)
                                                                 ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                     >  1:2983549/1‑40 (MQ=255)
                                                                     cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                   <  1:1522713/54‑1 (MQ=255)
                                                                     cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGcc  <  1:3523044/71‑1 (MQ=255)
                                                                     |                                                                      
AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCC  >  minE/2572648‑2572787

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: