Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,314,240 T→A 28.7% R165R (CGT→CGA 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,2400TA28.7% 175.3 / 81.0 132R165R (CGT→CGAhisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (45/49);  new base A (21/17);  total (66/66)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

ACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCT  >  minE/1314175‑1314302
                                                                 |                                                              
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:1943170/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:1115276/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:1670300/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:2726652/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:3372931/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:3149141/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:1366870/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:2630526/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:87529/1‑69 (MQ=255)
aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAc                                                             >  1:2622425/1‑69 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:1811112/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:1693282/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:228502/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:2551961/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:1422218/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:974933/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:2751777/69‑1 (MQ=255)
   ggTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                          <  1:1057508/69‑1 (MQ=255)
       tACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACg                                                      >  1:679170/1‑69 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:235518/68‑1 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:2920568/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:968678/68‑1 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1773522/68‑1 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:731789/68‑1 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:823737/68‑1 (MQ=255)
          gCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:827846/68‑1 (MQ=255)
                aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:863382/67‑1 (MQ=255)
                aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2015430/67‑1 (MQ=255)
                aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2393262/67‑1 (MQ=255)
                               aCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1487137/47‑1 (MQ=255)
                                cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCt                                                    <  1:1383668/46‑1 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:3071039/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1351958/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2531587/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2457868/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2372258/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:842111/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1294335/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:320169/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:2200629/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1176008/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:405993/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1004565/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:675360/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1850906/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1767672/3‑68 (MQ=255)
                                     cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                         >  1:1808150/3‑68 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTATGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:2625292/1‑68 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2680487/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:3001556/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:3340604/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:443951/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:512660/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:1001104/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2143869/42‑1 (MQ=255)
                                         gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                     >  1:506822/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)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                               <  1:2295515/41‑1 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:836983/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:463524/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:392819/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:1695059/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:3119884/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:843167/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2890087/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:1514384/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:1401765/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2827078/1‑69 (MQ=255)
                                          aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                   >  1:2783626/1‑69 (MQ=255)
                                            tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:1570672/1‑69 (MQ=255)
                                            tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                 >  1:1550899/1‑69 (MQ=255)
                                                 aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:2142160/69‑1 (MQ=255)
                                                 aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc            <  1:899019/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:2858445/69‑1 (MQ=255)
                                                        cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:2886370/69‑1 (MQ=255)
                                                        cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg     <  1:2246130/69‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAATAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:978431/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1575179/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:3122365/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1542223/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2255749/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2946927/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2791658/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2865435/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:596153/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2000525/66‑1 (MQ=255)
                                                           taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:2014797/66‑1 (MQ=255)
                                                            aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt  <  1:1460387/66‑1 (MQ=255)
                                                                 |                                                              
ACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCT  >  minE/1314175‑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: