Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,314,252 T→C 42.5% A169A (GCT→GCC hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,314,2520TC42.5% 26.4 / 72.2 94A169A (GCT→GCChisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base C (21/19);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCT  >  minE/1314185‑1314302
                                                                   |                                                  
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:968678/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:235518/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1773522/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2647601/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2693717/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2760401/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2820928/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2920568/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:827846/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:837632/68‑1 (MQ=255)
gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:731789/68‑1 (MQ=255)
      aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2015430/67‑1 (MQ=255)
      aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:85617/67‑1 (MQ=255)
      aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2393262/67‑1 (MQ=255)
      aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:863382/67‑1 (MQ=255)
                     aCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1487137/47‑1 (MQ=255)
                      cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1383668/46‑1 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:25853/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2200629/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1004565/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2372258/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2428536/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2457868/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2531587/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:842111/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:3071039/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:320169/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:675360/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2022661/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1351958/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1294335/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1767672/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1808150/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1850906/3‑68 (MQ=255)
                           cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1176008/3‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTATGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:2625292/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:1161370/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:512660/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:1001104/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:802545/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2811485/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:3340604/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2143869/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:443951/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:1656977/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2680487/42‑1 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:1546485/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:506822/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:467611/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:3278024/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:1717406/1‑68 (MQ=255)
                               gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:844075/1‑68 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGTTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2157673/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:463524/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:392819/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:3119884/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:1514384/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:384040/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:36162/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:3427525/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2890087/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:1401765/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:285864/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2827078/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2783626/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:836983/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:1695059/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:843167/1‑69 (MQ=255)
                                aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2248628/1‑69 (MQ=255)
                                  tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:645457/1‑69 (MQ=255)
                                  tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:1550899/1‑69 (MQ=255)
                                  tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:1570672/1‑69 (MQ=255)
                                  tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:2177280/1‑69 (MQ=255)
                                       aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:899019/69‑1 (MQ=255)
                                       aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:2142160/69‑1 (MQ=255)
                                       aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:220111/69‑1 (MQ=255)
                                       aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:3198566/69‑1 (MQ=255)
                                       aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTCACGTTc            <  1:435799/69‑1 (MQ=255)
                                              cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATGCCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:1304822/69‑1 (MQ=255)
                                              cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:2246130/69‑1 (MQ=255)
                                              cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:2858445/69‑1 (MQ=255)
                                              cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:2886370/69‑1 (MQ=255)
                                              cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:1397864/69‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCTCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1715190/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAATAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:978431/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2865435/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:596153/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2791658/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1575179/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2946927/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:3038616/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2000525/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2255749/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:3122365/66‑1 (MQ=255)
                                                 taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2014797/66‑1 (MQ=255)
                                                  aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1460387/66‑1 (MQ=255)
                                                  aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:914324/66‑1 (MQ=255)
                                                  aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:429962/66‑1 (MQ=255)
                                                            gTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2562875/58‑1 (MQ=255)
                                                               aCGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:808087/55‑1 (MQ=255)
                                                                   |                                                  
GCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCT  >  minE/1314185‑1314302

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: