Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,316:1 +C 9.5% intergenic (+84/‑1695) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3161.C 9.5% 342.2 / 38.4 95intergenic (+84/‑1695)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (49/37);  new base C (4/5);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 109 total aligned reads displayed.

TATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTT  >  minE/2571249‑2571385
                                                                    |                                                                      
tATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCtt                                                                     >  1:2544806/1‑71 (MQ=255)
   ggAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:819197/1‑71 (MQ=255)
    gAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:5241585/1‑70 (MQ=255)
      aCTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAatat                                                               <  1:1456148/71‑1 (MQ=255)
        tttGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCttt                                                                    <  1:1847942/64‑1 (MQ=255)
          tGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:494599/1‑64 (MQ=255)
          tGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCAt                                                           >  1:4035791/1‑71 (MQ=255)
            cTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCAtt                                                          <  1:3430721/70‑1 (MQ=255)
              ccTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATg                                                       >  1:2843252/1‑71 (MQ=255)
               cTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATAtt                                                              >  1:4395022/1‑63 (MQ=255)
               cTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATAtt                                                              >  1:1725519/1‑63 (MQ=255)
               cTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                      >  1:1917431/1‑71 (MQ=255)
                 tGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                      >  1:3885188/1‑69 (MQ=255)
                 tGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                      >  1:4428264/1‑69 (MQ=255)
                  gTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAata                                                                >  1:1343258/1‑58 (MQ=255)
                         cTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                            <  1:4524869/71‑1 (MQ=255)
                         cTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                            >  1:901165/1‑71 (MQ=255)
                          tAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCa                                                            <  1:4300056/54‑1 (MQ=255)
                          tAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCa                                                            <  1:5559931/54‑1 (MQ=255)
                           aaTTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         >  1:2212674/1‑71 (MQ=255)
                            attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATg                                                       >  1:486374/1‑57 (MQ=255)
                            attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTtaat                                                  <  1:1100764/62‑1 (MQ=255)
                            attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:1810774/70‑1 (MQ=255)
                            attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:4295718/70‑1 (MQ=255)
                            attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         >  1:3597404/1‑70 (MQ=255)
                             ttattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                              <  1:1623633/65‑1 (MQ=255)
                             ttattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                              <  1:4140373/65‑1 (MQ=255)
                              tattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                      >  1:1101002/1‑71 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:920494/1‑43 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:4532609/1‑43 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:5444221/1‑43 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTaa                                                                  >  1:702660/1‑43 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCa                                                            >  1:1684697/1‑49 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCAt                                                           >  1:1810523/1‑50 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                      >  1:2657581/1‑70 (MQ=255)
                               attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                      >  1:2151432/1‑70 (MQ=255)
                                tttttGCATATTACTTCATTAGCTTCGATGATTTTTCttt                                                                    <  1:3642417/40‑1 (MQ=255)
                                tttttGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                      >  1:3723745/1‑69 (MQ=255)
                                     gCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                               >  1:764500/1‑71 (MQ=255)
                                      cATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                               <  1:4267138/70‑1 (MQ=255)
                                        tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                            >  1:2722097/1‑71 (MQ=255)
                                        tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                            >  1:2352113/1‑71 (MQ=255)
                                         atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:3364244/57‑1 (MQ=255)
                                         atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                       >  1:4377557/1‑59 (MQ=255)
                                         atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                       >  1:4371958/1‑59 (MQ=255)
                                         atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                           >  1:4021442/1‑71 (MQ=255)
                                         atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                           >  1:5669681/1‑71 (MQ=255)
                                          ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTtaataa                                                >  1:4959213/1‑50 (MQ=255)
                                          ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCtctc                          <  1:2285486/71‑1 (MQ=255)
                                             cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTTTGTTAATAAGTTc                                            <  1:687733/51‑1 (MQ=255)
                                             cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                            <  1:4607018/51‑1 (MQ=255)
                                             cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTgca                                       >  1:5408570/1‑53 (MQ=255)
                                             cTTCATTAGCTTCGATGATTTTTATTTAATATTCATTATGTtaataa                                                >  1:905468/1‑47 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                              <  1:5675227/48‑1 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:1216093/52‑1 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:284437/52‑1 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:5198668/52‑1 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         <  1:692641/52‑1 (MQ=255)
                                              ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATg                                     <  1:5365132/56‑1 (MQ=255)
                                                cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                    <  1:1926925/55‑1 (MQ=255)
                                                cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGa                    >  1:2016698/1‑71 (MQ=255)
                                                 aTTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAt                   <  1:327776/71‑1 (MQ=255)
                                                  ttAGCTTCGATGATTTTTCTTTAATATTCATTATGtt                                                     >  1:5005218/1‑37 (MQ=255)
                                                  ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                      >  1:1240828/1‑51 (MQ=255)
                                                  ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                  >  1:598891/1‑71 (MQ=255)
                                                  ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                  >  1:4836689/1‑71 (MQ=255)
                                                   tAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGTAtt                  <  1:225473/70‑1 (MQ=255)
                                                    aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAttt                 <  1:2451652/70‑1 (MQ=255)
                                                       ttCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                    >  1:4025417/1‑48 (MQ=255)
                                                         tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGAtttctt              <  1:4252406/69‑1 (MQ=255)
                                                         tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCt                                          <  1:4702250/40‑1 (MQ=12)
                                                         tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCTgg                                        <  1:4008848/42‑1 (MQ=16)
                                                         cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                            <  1:4416152/39‑1 (MQ=255)
                                                         cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                          >  1:4491163/1‑40 (MQ=255)
                                                         cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                          >  1:3069491/1‑40 (MQ=255)
                                                          gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc          <  1:1022779/71‑1 (MQ=255)
                                                           atgatTTTTCCTTTAATATCCAGCATGTTaagaag                                               >  1:1241167/1‑35 (MQ=9)
                                                           atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTc                                            <  1:3842416/38‑1 (MQ=12)
                                                           atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTc                                            <  1:2605981/38‑1 (MQ=12)
                                                           atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATc           >  1:3425292/1‑71 (MQ=255)
                                                           atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCTgg                                        >  1:1356806/1‑41 (MQ=16)
                                                           atgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCAtt                       <  1:1042987/57‑1 (MQ=255)
                                                             gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         >  1:4875433/1‑37 (MQ=255)
                                                             gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc          <  1:3671395/68‑1 (MQ=255)
                                                             gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc       >  1:3264253/1‑71 (MQ=255)
                                                             gatTTTACTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc       >  1:734041/1‑71 (MQ=255)
                                                             gatATTTCTTTAATATTCATTATGTTAATAAGTTCTg                                         >  1:5530527/1‑37 (MQ=255)
                                                              tttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTc       >  1:5093691/2‑71 (MQ=255)
                                                              atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                       >  1:3172798/1‑38 (MQ=255)
                                                              atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc       <  1:2984344/70‑1 (MQ=255)
                                                               tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                               <  1:1745162/45‑1 (MQ=255)
                                                               tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                               <  1:4008930/45‑1 (MQ=255)
                                                               tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                               <  1:720328/45‑1 (MQ=255)
                                                               tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt     <  1:119029/71‑1 (MQ=255)
                                                               ttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt      >  1:1490904/1‑71 (MQ=255)
                                                                tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt      >  1:1847902/1‑70 (MQ=255)
                                                                 tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg    <  1:4212067/70‑1 (MQ=255)
                                                                 tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg    <  1:3487164/70‑1 (MQ=255)
                                                                 tttCTTTAATATTCATTATGTTAATAAGTTCTCGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt   >  1:5527992/1‑71 (MQ=255)
                                                                  ttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa             >  1:2233766/1‑60 (MQ=255)
                                                                    |                                                                      
TATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTT  >  minE/2571249‑2571385

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: