Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,678 G→A 7.6% intergenic (+1446/‑333) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,6780GA 7.6% 305.7 / 10.3 105intergenic (+1446/‑333)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (47/50);  new base A (3/5);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 106 total aligned reads displayed.

GAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCC  >  minE/2572609‑2572747
                                                                     |                                                                     
gAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGt                                                                      >  1:5455027/1‑71 (MQ=255)
  aaGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc                                                                    <  1:1623420/71‑1 (MQ=255)
   aGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCa                                                                   <  1:1116168/71‑1 (MQ=255)
         ccTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc                                                                    <  1:5359246/64‑1 (MQ=255)
         ccTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc                                                                    <  1:2057535/64‑1 (MQ=255)
            aaaGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACt                                                          >  1:3448779/1‑71 (MQ=255)
              aGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGa                                                        >  1:1399628/1‑71 (MQ=255)
                ttATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGAt                                                       >  1:1714589/1‑70 (MQ=255)
                     cTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAAc                                                                 >  1:701554/1‑55 (MQ=255)
                     cTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATa                                                            >  1:4163972/1‑60 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATaaa                                                    <  1:3060346/65‑1 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAtcc                                               <  1:4310736/70‑3 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTcc                                              <  1:2750938/71‑1 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTc                                               <  1:5405058/70‑1 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTc                                               <  1:3853401/70‑1 (MQ=255)
                        tataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTc                                               >  1:2850476/1‑70 (MQ=255)
                         ataAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAtt                                                >  1:5078795/1‑68 (MQ=255)
                          taAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCac                                            >  1:805141/1‑71 (MQ=255)
                            aaCAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATaaaaa                                                  >  1:1153860/1‑63 (MQ=255)
                            aaCAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacct                                          <  1:3079356/71‑1 (MQ=255)
                            aaCAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacct                                          <  1:3974861/71‑1 (MQ=255)
                               aGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGa                                                              <  1:4321657/48‑1 (MQ=255)
                                ggTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctacct                                      >  1:3298856/1‑71 (MQ=255)
                                 gTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTc                                                                    <  1:562847/40‑1 (MQ=255)
                                 gTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCtt                                     >  1:2416482/1‑71 (MQ=255)
                                  tGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    <  1:1357980/71‑1 (MQ=255)
                                    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCa                                             >  1:5379364/1‑60 (MQ=255)
                                    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    >  1:4964614/1‑69 (MQ=255)
                                    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    >  1:1142815/1‑69 (MQ=255)
                                    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                  <  1:1422235/71‑1 (MQ=255)
                                     cTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCaa                                                                  <  1:447350/38‑1 (MQ=255)
                                     cTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCaa                                                                  <  1:1561162/38‑1 (MQ=255)
                                      tAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAt                                                 <  1:2952006/54‑1 (MQ=255)
                                      tAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGc                                >  1:2094602/1‑71 (MQ=255)
                                      tAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGc                                >  1:49578/1‑71 (MQ=255)
                                       aaTCCCAAGGGGATCGATTGCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                               >  1:2985563/1‑71 (MQ=255)
                                       aaTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTc                                               >  1:4097145/1‑55 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGa                                                              >  1:2261199/1‑39 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGa                                                              >  1:2653905/1‑39 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGa                                                              >  1:5448251/1‑39 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTc                                               >  1:890351/1‑54 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                              <  1:5210802/71‑1 (MQ=255)
                                        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                               <  1:3723592/70‑1 (MQ=255)
                                         tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                              >  1:137772/1‑70 (MQ=255)
                                         tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                             <  1:2969901/71‑1 (MQ=255)
                                         tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                             <  1:5270177/71‑1 (MQ=255)
                                          cccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                            >  1:3999477/1‑71 (MQ=255)
                                          cccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                            >  1:4792942/1‑71 (MQ=255)
                                           ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAtt                                                <  1:5061806/50‑1 (MQ=255)
                                           ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                           >  1:1848093/1‑71 (MQ=255)
                                           ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACATACCTTCCAGCCCaa                            >  1:2296532/1‑70 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    <  1:2106058/61‑1 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    <  1:1032664/61‑1 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                    <  1:3296740/61‑1 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                          <  1:2416703/71‑1 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                          >  1:73009/1‑71 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                          >  1:5319246/1‑71 (MQ=255)
                                            cAAGGGGATCGATTCCAGAAATAGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                           >  1:55964/1‑70 (MQ=255)
                                             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                         >  1:1465919/1‑71 (MQ=255)
                                             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                         >  1:4834487/1‑71 (MQ=255)
                                               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCgaa                       <  1:724763/71‑3 (MQ=255)
                                               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                         >  1:2091529/1‑69 (MQ=255)
                                               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                         >  1:3756246/1‑69 (MQ=255)
                                                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATaaa                                                    <  1:3238734/39‑1 (MQ=255)
                                                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                             <  1:2036989/62‑1 (MQ=255)
                                                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGcc                    >  1:1086747/1‑71 (MQ=255)
                                                   aTCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                   <  1:423499/71‑1 (MQ=255)
                                                      gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGa                 <  1:1099265/70‑1 (MQ=255)
                                                       aTTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCaccta                                         <  1:3775956/45‑1 (MQ=255)
                                                        ttCCAGAAATCGGGTTCAACGGATACTGATAAAAAt                                                 <  1:1892832/36‑1 (MQ=255)
                                                        ttCCAGAAATCGGGTTCAACGGATACTGATAAAAAt                                                 <  1:2648276/36‑1 (MQ=255)
                                                         tCCAGAAATCGGGTTCAACGGATACTGATAAAAATTcc                                              <  1:4185870/38‑1 (MQ=255)
                                                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                  >  1:4658891/1‑49 (MQ=255)
                                                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                  >  1:5590047/1‑49 (MQ=255)
                                                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                  >  1:1602066/1‑49 (MQ=255)
                                                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt             <  1:4341378/70‑1 (MQ=255)
                                                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt             <  1:1882051/70‑1 (MQ=255)
                                                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta           <  1:4871295/71‑1 (MQ=255)
                                                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta           >  1:28989/1‑71 (MQ=255)
                                                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta           <  1:2017498/71‑1 (MQ=255)
                                                             gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat          <  1:1170903/67‑1 (MQ=255)
                                                             gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat          <  1:3519381/67‑1 (MQ=255)
                                                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                             <  1:2400665/51‑1 (MQ=255)
                                                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg              <  1:3999812/66‑1 (MQ=255)
                                                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata         <  1:1731037/71‑1 (MQ=255)
                                                              agaTCGGATTCAACGGATACTGATAAAAATTCCATCCTCCTTTCAGCCCAATCGGATCCTGAGTGAtata         >  1:4143937/3‑70 (MQ=255)
                                                              aaaTCGGGTTCAACGGATACTGATAAAAATTCCacctacct                                      >  1:300929/1‑41 (MQ=255)
                                                              aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGccc                   >  1:5141368/1‑59 (MQ=255)
                                                              aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                  >  1:2153154/1‑61 (MQ=255)
                                                                aTCTGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt               <  1:877913/62‑1 (MQ=255)
                                                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTcc                                   >  1:274741/1‑42 (MQ=255)
                                                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt               <  1:4897259/62‑1 (MQ=255)
                                                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCg      <  1:1734648/71‑1 (MQ=255)
                                                                 tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                             <  1:4050751/47‑1 (MQ=255)
                                                                 tCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGatc                     <  1:2424029/55‑3 (MQ=255)
                                                                 ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga             <  1:4494340/62‑2 (MQ=255)
                                                                 ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc     >  1:1660448/2‑71 (MQ=255)
                                                                 ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc     >  1:1231848/2‑71 (MQ=255)
                                                                  cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAg                                 <  1:171590/42‑1 (MQ=255)
                                                                    ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc  <  1:2673676/71‑1 (MQ=255)
                                                                     |                                                                     
GAAAGCATACCTAAAGTTATGCTGTATAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCC  >  minE/2572609‑2572747

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: