Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,138,447 G→A 36.8% L182L (TTG→TTA ydjM → predicted inner membrane protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,138,4470GA36.8% 74.5 / 70.1 114L182L (TTG→TTAydjMpredicted inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base G (51/21);  new base A (25/17);  total (76/38)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 114 total aligned reads displayed.

GATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCAAAAATCACCTTT  >  minE/1138384‑1138511
                                                               |                                                                
gATGCCCCATTCATTACCCGAGAAGAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:274129/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCa                                                                >  1:1045480/1‑66 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:3128530/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:676369/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:3345055/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:1130312/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGa                                                              >  1:3258397/1‑68 (MQ=255)
gATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCAACGCAAATGATCAATACCTTGCAg                                                               >  1:2320035/1‑67 (MQ=255)
  tGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATcc                                                           <  1:2997942/69‑1 (MQ=255)
        attcattACCCGAAAACAGAGCTGTTCGTTGGTCATCGCAAATGATAATACCTTACAGATCCAGTTTCa                                                    <  1:1266655/69‑1 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:2015636/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:3342727/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:3250096/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:3009701/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1832741/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:3075470/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:708490/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1616830/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1297797/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1294925/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1517227/1‑69 (MQ=255)
         ttcattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa                                                    >  1:1293352/1‑69 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:1109259/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:2155671/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:3038502/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:2407157/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:2507251/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:2641839/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:901420/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:3193761/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:732250/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:298581/69‑1 (MQ=255)
         ttcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:3159974/69‑1 (MQ=255)
          tcattACCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:482360/68‑1 (MQ=255)
           tattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCAt                                                   >  1:1025661/2‑68 (MQ=255)
           tattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCAt                                                   >  1:2838014/2‑68 (MQ=255)
           tattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATc                                                  >  1:3420500/2‑69 (MQ=255)
           tattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATc                                                  >  1:1396948/2‑69 (MQ=255)
           cattACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCAt                                                   >  1:708999/1‑68 (MQ=255)
               aCCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGttt                                                      >  1:1925147/1‑61 (MQ=255)
               aCCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTc                                                     >  1:3396498/1‑62 (MQ=255)
               aCCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTc                                                     >  1:867911/1‑62 (MQ=255)
               aCCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTc                                                     >  1:668494/1‑62 (MQ=255)
               aCCCGAAAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:1406311/63‑1 (MQ=255)
                       aCAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGc                                      <  1:1591552/69‑1 (MQ=255)
                          gcgcTGTTCGTTGGTCATCGCAAATGATCAATACCTTACAGATCCAGTTTCa                                                    <  1:118677/52‑1 (MQ=255)
                                  cGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTg                             <  1:3151039/67‑1 (MQ=255)
                                  cGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTg                             <  1:1516760/67‑1 (MQ=255)
                                  cGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTg                             <  1:1694420/67‑1 (MQ=255)
                                  cGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTg                             <  1:285598/67‑1 (MQ=255)
                                   gTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTg                             <  1:1068350/66‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:245124/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:3300866/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:1554987/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:2348446/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:358780/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:3211805/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGaa                           <  1:354794/58‑1 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:2381308/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:3170706/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:2717632/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:1215536/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:2465869/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:1940520/1‑68 (MQ=255)
                                             gCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAg                 >  1:818958/1‑68 (MQ=255)
                                               aaaTGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCa              <  1:1026199/69‑1 (MQ=255)
                                               aaaTGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCa              <  1:1698441/69‑1 (MQ=255)
                                               aaaTGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCa              <  1:3029392/69‑1 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:2922793/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:564995/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:1610509/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:2284760/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:825370/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaaa           >  1:2868969/1‑69 (MQ=255)
                                                  tGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCaaa            >  1:260926/1‑68 (MQ=255)
                                                   gATCAATAGCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:3200707/1‑41 (MQ=38)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAg                                       >  1:2171314/1‑40 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1482266/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:824990/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1588583/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1592427/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1130592/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1755850/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1862938/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:646516/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:639233/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:937747/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:559340/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1933521/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:376152/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:1945515/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:219484/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:2270404/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:3289121/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:3285270/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:2478951/1‑41 (MQ=255)
                                                   gATCAATACCTTACAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:2706899/1‑41 (MQ=255)
                                                   gATCAATACCTTAAAGATCCAGTTTCATCGGCTTATTAAGc                                      >  1:2339070/1‑41 (MQ=38)
                                                           ccTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCAAAAATCACCttt  >  1:2474973/1‑69 (MQ=255)
                                                           ccTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCAAAAATCACCttt  >  1:2796960/1‑69 (MQ=255)
                                                               |                                                                
GATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCATCGGCTTATTAAGCATCAGGTTGAATACTAAAAAGGCAAAAATCACCTTT  >  minE/1138384‑1138511

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: