Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 329,715 T→C 40.4% H374H (CAT→CAC gsk → inosine/guanosine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE329,7150TC40.4% 41.1 / 72.6 103H374H (CAT→CACgskinosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base T (36/25);  new base C (20/21);  total (57/46)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 105 total aligned reads displayed.

AATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGT  >  minE/329653‑329783
                                                              |                                                                    
aaTGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGc                                                                >  1:1631544/1‑69 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:2938083/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:1814986/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:1829086/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:2142397/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:237247/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:245314/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:2710234/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:2805826/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:1068654/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:2953494/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:415654/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:509792/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:574844/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:605378/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:609783/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:756467/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:762512/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:1584407/58‑1 (MQ=255)
      gcaggggATGGCGCACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:1061583/58‑1 (MQ=255)
      gcaggggATGGCGAACTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTa                                                                  <  1:3041251/58‑1 (MQ=255)
         ggggATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCa                                                       >  1:3395284/1‑69 (MQ=255)
         ggggATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCa                                                       >  1:2315230/1‑69 (MQ=255)
         ggggATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCa                                                       >  1:2973181/1‑69 (MQ=255)
                gcgcATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTcc                                                 >  1:152450/1‑68 (MQ=255)
                       ggCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACg                                                            >  1:3223150/1‑50 (MQ=255)
                       ggCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAAc                                          >  1:3199641/1‑68 (MQ=255)
                       ggCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAAc                                          >  1:51937/1‑68 (MQ=255)
                       ggCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAAc                                          >  1:3383996/1‑68 (MQ=255)
                       ggCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAAc                                          >  1:1681886/1‑68 (MQ=255)
                          aGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATa                                       <  1:325095/68‑1 (MQ=255)
                          aGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATa                                       <  1:2305950/68‑1 (MQ=255)
                           gCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATaaa                                     <  1:2742866/69‑1 (MQ=255)
                              ttGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTc                                  <  1:139090/69‑1 (MQ=255)
                                gctgcATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAc                                >  1:2118252/1‑69 (MQ=255)
                                          aTTACCGCCAACATCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:1778320/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:215981/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:484616/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:3269915/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:196466/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:2547315/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:1105437/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:1880043/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:3150773/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:846969/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:1733783/69‑1 (MQ=255)
                                          aTTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:2757949/69‑1 (MQ=255)
                                           ttACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:625561/68‑1 (MQ=255)
                                            tACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    <  1:845930/69‑1 (MQ=255)
                                             aCCGCCAACAGCTACGATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:385189/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTa                         <  1:879383/63‑1 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTa                         <  1:1416779/63‑1 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTa                         <  1:1837578/63‑1 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTa                         <  1:2742006/63‑1 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1229869/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1551513/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:3272102/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2333526/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:3421966/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2248863/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:464731/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1592045/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1341216/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1493175/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2024936/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:654362/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:681726/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:707382/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1836629/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:3218628/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:1530101/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2947763/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:3181027/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2961932/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:3154440/1‑68 (MQ=255)
                                             aCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTa                    >  1:2678708/1‑68 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCGCCt                              >  1:2556324/1‑54 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:189019/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:3011834/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:637050/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:2968159/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:2902910/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:52431/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:1873003/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:787607/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:2936020/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:847972/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:1606134/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:2062074/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:487347/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:26777/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:3395051/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:3204792/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTaa                        >  1:2435388/1‑60 (MQ=255)
                                                 ccAACAGCTACCACCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTa                         >  1:2159611/1‑59 (MQ=255)
                                                 ccAACAGCTACCAACGTAGCAACGTACCAAACTCCAGCAAACATAAATACACCTGGTTaa                        >  1:2256996/1‑60 (MQ=255)
                                                    aCAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACt                      <  1:1165992/59‑1 (MQ=255)
                                                             aTCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAgg   >  1:42365/1‑69 (MQ=255)
                                                              tCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgt  <  1:1118548/69‑1 (MQ=255)
                                                              tCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgt  <  1:3077679/69‑1 (MQ=255)
                                                              |                                                                    
AATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGT  >  minE/329653‑329783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: