Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,418 T→A 16.7% intergenic (+186/‑1593) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,4180TA16.7% 339.4 / 54.6 138intergenic (+186/‑1593)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (62/53);  new base A (11/12);  total (73/65)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 139 total aligned reads displayed.

GAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATC  >  minE/2571350‑2571488
                                                                    |                                                                      
gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                                      >  1:463623/1‑71 (MQ=255)
gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                                      >  1:1537116/1‑71 (MQ=255)
 aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                                      <  1:1820505/70‑1 (MQ=255)
     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                                                 <  1:4960781/71‑1 (MQ=255)
     tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:4984534/70‑1 (MQ=255)
       ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                                               <  1:151787/71‑1 (MQ=255)
       atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                                                      <  1:1855778/63‑1 (MQ=255)
       atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                                                      <  1:4993226/63‑1 (MQ=255)
         ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                                             >  1:1943869/1‑71 (MQ=255)
          tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             >  1:4314960/1‑70 (MQ=255)
          tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                                                 <  1:1135986/66‑4 (MQ=255)
           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              <  1:1937277/68‑1 (MQ=255)
           cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                                           >  1:2811266/1‑71 (MQ=255)
            aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:3581123/68‑1 (MQ=255)
            aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                          <  1:3135749/71‑1 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                                              <  1:4527043/67‑1 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                           <  1:2542554/70‑3 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                                           <  1:4734283/70‑3 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:1715371/1‑71 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:3556466/1‑71 (MQ=255)
            aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                          >  1:2753796/1‑71 (MQ=255)
             ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                                                    <  1:1107040/60‑1 (MQ=255)
             ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:5173566/62‑1 (MQ=255)
               ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:4112445/1‑62 (MQ=255)
               ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTTTGCATTGTTGAGAACGCACTAATGACATGa                                                        <  1:5040939/70‑1 (MQ=255)
               ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                                        <  1:947627/70‑1 (MQ=255)
                gATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgt                                                      >  1:383292/1‑71 (MQ=255)
                gATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGaat                                                      >  1:3122013/1‑69 (MQ=255)
                  ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                                                         <  1:231130/66‑1 (MQ=255)
                  ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                                                         <  1:2994731/66‑1 (MQ=255)
                  tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                                    >  1:1970997/1‑55 (MQ=255)
                  tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                                    >  1:5134465/1‑55 (MQ=255)
                  tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                                    >  1:4599096/1‑55 (MQ=255)
                   tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGc                                                                     <  1:4813418/53‑1 (MQ=255)
                   tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                                              >  1:2172702/1‑60 (MQ=255)
                    ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt                                                  >  1:622813/1‑71 (MQ=255)
                     tttATCGGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                                   >  1:731348/1‑53 (MQ=255)
                       tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                                                <  1:2406726/70‑1 (MQ=255)
                         tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:5384567/55‑1 (MQ=255)
                         tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:3553576/55‑1 (MQ=255)
                         tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:1538276/55‑1 (MQ=255)
                           cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           >  1:3330444/1‑71 (MQ=255)
                           cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           >  1:5491644/1‑71 (MQ=255)
                             tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           <  1:3229031/69‑1 (MQ=255)
                             tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                           <  1:5372862/69‑1 (MQ=255)
                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                       >  1:1955615/1‑71 (MQ=255)
                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                       >  1:1260997/1‑71 (MQ=255)
                               ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg                                      >  1:1625714/1‑71 (MQ=255)
                                 gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                                       <  1:1230706/69‑1 (MQ=255)
                                  ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:4301107/46‑1 (MQ=255)
                                  ttCTGCCCCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             >  1:4292406/1‑46 (MQ=255)
                                   tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                                   >  1:2466197/1‑39 (MQ=255)
                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                             <  1:1971366/44‑1 (MQ=255)
                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                  >  1:4362226/1‑71 (MQ=255)
                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                                  >  1:5671987/1‑71 (MQ=255)
                                    cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct                                   >  1:3810355/1‑70 (MQ=255)
                                     tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:3935302/38‑1 (MQ=255)
                                     tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                                  <  1:64923/38‑1 (MQ=255)
                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:2298830/1‑38 (MQ=255)
                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:4033043/1‑38 (MQ=255)
                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                                >  1:1664118/1‑38 (MQ=255)
                                       ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct                                 <  1:5208430/69‑1 (MQ=255)
                                        cGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCg                                         >  1:3029933/1‑60 (MQ=255)
                                          cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            <  1:4786608/71‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                                     <  1:3708620/43‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                                     <  1:2142377/43‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTa                                <  1:5339379/64‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:258067/71‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:5505732/71‑1 (MQ=255)
                                             gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:3290580/71‑1 (MQ=255)
                                              aTGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                                     <  1:4218532/42‑1 (MQ=255)
                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                               >  1:5210211/1‑63 (MQ=255)
                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         >  1:2546351/1‑69 (MQ=255)
                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:711638/1‑66 (MQ=255)
                                               tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:4208551/1‑66 (MQ=255)
                                               tGAGACGATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                               >  1:1459/1‑63 (MQ=255)
                                                 agaCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                      <  1:3733606/70‑1 (MQ=255)
                                                   aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa                       >  1:4354855/1‑67 (MQ=255)
                                                   aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAagag                    >  1:242859/1‑70 (MQ=255)
                                                     tataCATCGTCGAGATCGCAATGATGTTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 >  1:3801726/1‑71 (MQ=255)
                                                       taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCTACGGCTCTAATCAAAAAAGAGCAgg                <  1:3670887/70‑1 (MQ=255)
                                                       taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAAt                              <  1:5045480/56‑1 (MQ=255)
                                                       taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg                 >  1:2680435/1‑69 (MQ=255)
                                                         cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:3092649/1‑56 (MQ=255)
                                                         cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                            >  1:4450503/1‑56 (MQ=255)
                                                         cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGt             >  1:4219500/1‑71 (MQ=255)
                                                          aTCGTCGAGATCGCAATGATGATATGAGTTGTtcct                                               >  1:4080725/1‑36 (MQ=255)
                                                           tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:1351916/1‑71 (MQ=255)
                                                           tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:224432/1‑71 (MQ=255)
                                                           tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:4684665/1‑71 (MQ=255)
                                                           tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:4018276/1‑71 (MQ=255)
                                                           tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc           >  1:3781592/1‑71 (MQ=255)
                                                                gagaTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAg              >  1:5419034/1‑63 (MQ=255)
                                                                  gaTCGCAATGATGATATGAGTTGTTCCTCCTGCGAAGGCTCTaa                               >  1:1148926/1‑44 (MQ=255)
                                                                   aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:2593939/49‑1 (MQ=255)
                                                                   aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                      <  1:4729998/52‑1 (MQ=255)
                                                                   aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag                      <  1:124382/52‑1 (MQ=255)
                                                                   aTCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAgg                <  1:84878/58‑1 (MQ=255)
                                                                   aTCGCAATCATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                         <  1:4850145/49‑1 (MQ=255)
                                                                    tCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATc  <  1:4509072/71‑1 (MQ=255)
                                                                    |                                                                      
GAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGCTTTTCCATC  >  minE/2571350‑2571488

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: