Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,735 G→A 7.8% intergenic (+1503/‑276) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7350GA 7.8% 408.5 / 16.7 141intergenic (+1503/‑276)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (77/53);  new base A (6/5);  total (83/58)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 159 total aligned reads displayed.

CCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAT  >  minE/2572667‑2572802
                                                                    |                                                                   
ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    <  1:1882051/70‑1 (MQ=255)
ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    <  1:4341378/70‑1 (MQ=255)
 cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                  <  1:2017498/71‑1 (MQ=255)
 cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                  >  1:28989/1‑71 (MQ=255)
   gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                                 <  1:3519381/67‑1 (MQ=255)
   gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                                 <  1:1170903/67‑1 (MQ=255)
   gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                                <  1:1731037/71‑1 (MQ=255)
    agaTCGGATTCAACGGATACTGATAAAAATTCCATCCTCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                                >  1:4143937/3‑70 (MQ=255)
       ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                                    <  1:4494340/62‑2 (MQ=255)
       ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                                            >  1:1231848/2‑71 (MQ=255)
       ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                                            >  1:1660448/2‑71 (MQ=255)
          ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc                                                         <  1:2673676/71‑1 (MQ=255)
           attCAGCGGATACGCATACAGACTCCATCCCCCCTCCAGCCCTATCGGGTCCTGAGTGATATACCGCCCct                                                        <  1:2205840/70‑1 (MQ=255)
            ttCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCctg                                                       >  1:4336450/1‑71 (MQ=255)
            ttCAACGTATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                  >  1:2577749/1‑60 (MQ=255)
            ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                  <  1:1139060/60‑1 (MQ=255)
              cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                                       >  1:3173499/1‑67 (MQ=255)
              cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                  >  1:1947126/1‑58 (MQ=255)
              aaaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    >  1:711242/2‑56 (MQ=255)
               aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGCCCctgcag                                                    >  1:1521892/1‑69 (MQ=255)
               aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                                    >  1:3513883/1‑71 (MQ=255)
                aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                                    >  1:3815513/1‑68 (MQ=255)
                aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                                    >  1:554642/1‑68 (MQ=255)
                 cGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCg                                                  <  1:3995756/71‑1 (MQ=255)
                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    <  1:1588314/53‑1 (MQ=255)
                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                                  >  1:5666993/1‑71 (MQ=255)
                  ggATACTGATAAAAATTCCATCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCCTGCAGCggg                                                >  1:2878568/1‑70 (MQ=255)
                  ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                                    <  1:5222836/52‑2 (MQ=255)
                  ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                                 >  1:5171066/1‑71 (MQ=255)
                  ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    <  1:1911403/52‑1 (MQ=255)
                  ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                    <  1:241324/52‑1 (MQ=255)
                  ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                               >  1:3432235/1‑57 (MQ=255)
                  ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                               >  1:1156018/1‑57 (MQ=255)
                   gATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGccc                                                          <  1:175360/61‑1 (MQ=255)
                   gATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                                 >  1:5212955/1‑70 (MQ=255)
                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                                >  1:2729850/1‑70 (MQ=255)
                     tACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                                    <  1:329436/65‑3 (MQ=255)
                       cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                                          >  1:4294827/1‑57 (MQ=255)
                        tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                               <  1:4968894/67‑1 (MQ=255)
                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                                 >  1:1934447/1‑64 (MQ=255)
                          aTAAAAATTCCATCCCGCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                               >  1:3379299/1‑48 (MQ=255)
                          aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                               >  1:5594015/1‑48 (MQ=255)
                          aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                               >  1:2503737/1‑48 (MQ=255)
                          aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                               >  1:1130632/1‑48 (MQ=255)
                          aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                                >  1:3729263/1‑48 (MQ=255)
                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                               >  1:4383217/1‑49 (MQ=255)
                            aaaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGatat                                                                 >  1:4746828/1‑45 (MQ=255)
                             aaaaTTCCATCACCCCTTCAGCCCAATCGGATCCTGAGTGAtata                                                                >  1:5468939/1‑45 (MQ=255)
                             aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCATGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAcc                                      >  1:4425874/1‑71 (MQ=255)
                               aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                                  >  1:1699164/1‑41 (MQ=255)
                                   ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                 <  1:5065269/70‑1 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                <  1:967637/70‑1 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                <  1:2640986/70‑1 (MQ=255)
                                     acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                                   <  1:309814/50‑1 (MQ=255)
                                       ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                                   <  1:399991/48‑1 (MQ=255)
                                       ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                                  >  1:2150441/1‑49 (MQ=255)
                                        tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                                      <  1:3548074/44‑1 (MQ=255)
                                        tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                                      <  1:158107/44‑1 (MQ=255)
                                        tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                                      <  1:5145287/44‑1 (MQ=255)
                                        tacctTCCAGCCCAATCGGAGCCTGAGAGTTATACCGCCCCAGCTGCGGTTAAAGTTACCGAAAATAGAta                           >  1:2432791/1‑71 (MQ=255)
                                         acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                           >  1:2439551/1‑70 (MQ=255)
                                          cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                           <  1:4763221/69‑1 (MQ=255)
                                          cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                           <  1:5310379/69‑1 (MQ=255)
                                            ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa                                        >  1:815782/1‑54 (MQ=255)
                                              tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                         <  1:1132602/34‑1 (MQ=37)
                                              tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                         <  1:4242852/34‑1 (MQ=37)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                                      >  1:274979/1‑38 (MQ=255)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt                       <  1:5113808/69‑1 (MQ=255)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                      <  1:5045414/70‑1 (MQ=255)
                                               cAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                                 >  1:701797/1‑42 (MQ=255)
                                               cAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                                 >  1:3910823/1‑42 (MQ=255)
                                               cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                      >  1:4094487/1‑69 (MQ=255)
                                               cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                      >  1:1811317/1‑69 (MQ=255)
                                               cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                      >  1:2144679/1‑69 (MQ=255)
                                               cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                     <  1:2124118/70‑1 (MQ=255)
                                                  ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                 <  1:1764035/71‑1 (MQ=255)
                                                  ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                 >  1:2273757/1‑71 (MQ=255)
                                                  cccAATCGGGTCCGGGGTGATATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcc                 >  1:64164/1‑71 (MQ=255)
                                                   cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc                >  1:1181916/4‑71 (MQ=255)
                                                   ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc                >  1:2761371/1‑71 (MQ=255)
                                                    cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt               <  1:28268/71‑1 (MQ=255)
                                                     caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg              <  1:2543541/70‑1 (MQ=255)
                                                     caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg              <  1:3753246/70‑1 (MQ=255)
                                                     aaTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                             <  1:4299080/40‑1 (MQ=255)
                                                      aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa             >  1:3363095/1‑71 (MQ=255)
                                                         ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                               >  1:4172927/1‑50 (MQ=255)
                                                            tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc       <  1:5675327/70‑1 (MQ=255)
                                                            tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc       <  1:2827269/70‑1 (MQ=255)
                                                             ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatca      >  1:1458927/1‑71 (MQ=255)
                                                               tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                     >  1:3606836/1‑54 (MQ=255)
                                                               tGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat     >  1:4848809/1‑70 (MQ=255)
                                                                gggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcatca   <  1:3320751/69‑1 (MQ=255)
                                                                gAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                              <  1:5317765/44‑1 (MQ=255)
                                                                 aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                              <  1:375271/43‑1 (MQ=255)
                                                                 aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca   <  1:4024081/70‑1 (MQ=255)
                                                                 aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat  <  1:9140/71‑1 (MQ=255)
                                                                  gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt           >  1:1850727/1‑61 (MQ=255)
                                                                  gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt           >  1:3075723/1‑61 (MQ=255)
                                                                  gtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcat     <  1:4888930/67‑1 (MQ=255)
                                                                    gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat     <  1:3995348/65‑1 (MQ=255)
                                                                    |                                                                   
CCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAT  >  minE/2572667‑2572802

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: