Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,736 T→A 12.7% intergenic (+1504/‑275) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7360TA12.7% 285.1 / 29.8 110intergenic (+1504/‑275)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/60);  new base A (10/4);  total (46/64)
Fisher's exact test for biased strand distribution p-value = 2.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 131 total aligned reads displayed.

CCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTA  >  minE/2572667‑2572806
                                                                     |                                                                      
ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                        <  1:779791/70‑1 (MQ=255)
ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                        <  1:3787628/70‑1 (MQ=255)
 cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                      <  1:390349/71‑1 (MQ=255)
  aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                                                     <  1:914487/71‑1 (MQ=255)
  aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                                                     <  1:131062/71‑1 (MQ=255)
   gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                                    >  1:3445080/1‑71 (MQ=255)
     aaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                                        <  1:1237055/65‑1 (MQ=255)
        cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                                                      >  1:1617017/1‑64 (MQ=255)
         ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                                              >  1:3595513/4‑71 (MQ=255)
           gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                                            <  1:782865/71‑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:3689940/1‑71 (MQ=255)
              cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                                         >  1:3488594/1‑71 (MQ=255)
              cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                                         >  1:2362725/1‑71 (MQ=255)
               aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                                                         <  1:3629112/70‑2 (MQ=255)
               aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                                                         <  1:3263023/70‑2 (MQ=255)
                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                                      <  1:1144066/71‑1 (MQ=255)
                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                                      >  1:1169944/1‑71 (MQ=255)
                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGATCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                                   <  1:2890548/71‑1 (MQ=255)
                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                                    >  1:532667/1‑70 (MQ=255)
                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                                    >  1:1447469/1‑70 (MQ=255)
                     tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                      >  1:86929/1‑51 (MQ=255)
                     tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                                   >  1:1461899/1‑70 (MQ=255)
                       cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                                                   >  1:2670562/1‑52 (MQ=255)
                       cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                                 <  1:364774/70‑1 (MQ=255)
                        tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                                      <  1:2543605/48‑1 (MQ=255)
                        tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                                   <  1:3353523/67‑1 (MQ=255)
                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                                   >  1:2645672/1‑50 (MQ=255)
                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                                   >  1:1525899/1‑50 (MQ=255)
                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                                                   >  1:346663/1‑50 (MQ=255)
                         gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                                    >  1:2062969/1‑49 (MQ=255)
                         gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                                   <  1:504292/66‑1 (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:1950267/71‑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)
                           tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa                                            <  1:3945637/71‑1 (MQ=255)
                            aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                                    >  1:763830/1‑62 (MQ=255)
                             aaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                                      >  1:1564126/1‑43 (MQ=255)
                             aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                                        <  1:610657/57‑1 (MQ=255)
                                 ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat                                       >  1:3399694/1‑70 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                                   <  1:612605/71‑1 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                    <  1:691148/70‑1 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                    <  1:686262/70‑1 (MQ=255)
                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                    <  1:1454696/70‑1 (MQ=255)
                                       ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                                                <  1:1095040/55‑1 (MQ=255)
                                       ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                                 >  1:1536064/1‑70 (MQ=255)
                                       ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                                 >  1:1540990/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:2652447/69‑1 (MQ=255)
                                           ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                                 >  1:2355988/1‑50 (MQ=255)
                                           ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                            <  1:2932678/71‑1 (MQ=255)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                                        >  1:2983549/1‑40 (MQ=255)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                             >  1:1564481/1‑51 (MQ=255)
                                              ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                           <  1:650788/69‑1 (MQ=255)
                                                  cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                                      <  1:2540480/54‑1 (MQ=255)
                                                  cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                                      <  1:1522713/54‑1 (MQ=255)
                                                  cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGcc                     <  1:3523044/71‑1 (MQ=255)
                                                  cccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGcc                     <  1:1406826/71‑1 (MQ=255)
                                                    cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCg                                         <  1:1850087/49‑1 (MQ=255)
                                                    cAATAGGAGCCTGAGTGTTATACCGCCCCTGCTGAGGTTAATGTTACCGATAATAGTta                               <  1:3474078/59‑1 (MQ=39)
                                                       tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                                 <  1:3220334/38‑1 (MQ=255)
                                                       tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                                      >  1:687402/1‑49 (MQ=255)
                                                       tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                               <  1:2347501/56‑1 (MQ=255)
                                                       tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                <  1:3828223/71‑1 (MQ=255)
                                                        cGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATg                                               >  1:2501741/1‑39 (MQ=255)
                                                        cGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt               >  1:1954597/1‑71 (MQ=255)
                                                         ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt               >  1:1606653/1‑70 (MQ=255)
                                                         ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc              <  1:3624089/71‑1 (MQ=255)
                                                         ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc              <  1:2309267/71‑1 (MQ=255)
                                                          gAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtca             <  1:2847561/71‑1 (MQ=255)
                                                          gAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtca             <  1:2707886/71‑1 (MQ=255)
                                                             ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                                              <  1:3401310/35‑1 (MQ=255)
                                                             ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt                                              <  1:2450298/35‑1 (MQ=255)
                                                              cTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat         <  1:2675552/71‑1 (MQ=255)
                                                               tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtat                              <  1:311973/49‑1 (MQ=255)
                                                                 aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGtt                                >  1:1370847/1‑45 (MQ=255)
                                                                 aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca       >  1:167906/1‑70 (MQ=255)
                                                                  gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcat            >  1:240330/1‑64 (MQ=255)
                                                                  gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta     <  1:868433/71‑1 (MQ=255)
                                                                   tgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                                  <  1:3809870/41‑1 (MQ=255)
                                                                    gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc              <  1:1914164/60‑1 (MQ=255)
                                                                     ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                           <  1:474436/46‑1 (MQ=255)
                                                                     ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTg                  <  1:2260319/55‑1 (MQ=255)
                                                                     ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc              <  1:533125/59‑1 (MQ=255)
                                                                     ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc              <  1:3991187/59‑1 (MQ=255)
                                                                     ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtacta  >  1:3112979/1‑71 (MQ=255)
                                                                     ttATACCGACCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtacta  >  1:2729023/1‑71 (MQ=255)
                                                                     |                                                                      
CCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTA  >  minE/2572667‑2572806

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: