Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,884,215 T→C 25.6% A131A (GCA→GCG cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,884,2150TC25.6% 131.4 / 43.2 113A131A (GCA→GCGcpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base C (16/13);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

AACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGTA  >  minE/2884153‑2884281
                                                              |                                                                   
aaCATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAg                                                                >  1:3430864/1‑68 (MQ=255)
aaCATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAg                                                                >  1:2697997/1‑68 (MQ=255)
aaCATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAg                                                                >  1:2270418/1‑68 (MQ=255)
aaCATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAg                                                                >  1:143650/1‑68 (MQ=255)
aaCATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAg                                                                >  1:1210456/1‑68 (MQ=255)
         tGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTc                                                     <  1:3002649/69‑1 (MQ=255)
                          aaTGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTa                                    >  1:691740/1‑69 (MQ=255)
                              aCGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:2699372/69‑1 (MQ=255)
                              aCGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:1449549/69‑1 (MQ=255)
                              aCGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:2355411/69‑1 (MQ=255)
                              aCGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:1583995/69‑1 (MQ=255)
                              aCGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:579496/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCTCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:2248782/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:889106/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:2030050/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:3100550/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:3403236/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:1874056/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:45748/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:1741061/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:1581303/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:2453136/69‑1 (MQ=255)
                                 ttGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:2482035/69‑1 (MQ=255)
                                  tGGCATTTACATAAGGGAATTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAAc                             <  1:1327225/68‑1 (MQ=255)
                                   ggCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGt                                     >  1:2158928/1‑59 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGt                                                 >  1:2923329/1‑41 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGt                                                 >  1:807538/1‑41 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGt                                                 >  1:319060/1‑41 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2818508/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2804854/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2923648/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2734686/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2634381/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2157527/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:2560452/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:651350/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:687345/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:1481373/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:713073/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:1415743/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:1388624/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:130886/1‑69 (MQ=255)
                                         tACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGtt                     >  1:804895/1‑69 (MQ=255)
                                               aGGGAATTTCGCTCCTGCCAGCGCATTTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCaa                >  1:2503444/1‑69 (MQ=255)
                                               aGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:520419/52‑1 (MQ=255)
                                               aGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTg                        >  1:2086780/1‑60 (MQ=255)
                                                 ggAATTTCGCTCCTGCCAGCTCATTTTTCAGGTAATCCAGACCGTAGTTa                                <  1:790752/50‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:2628293/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            >  1:2611849/1‑69 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:3079099/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:3418441/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:813984/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:1021171/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:1671954/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:1937481/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGt            <  1:1956058/69‑1 (MQ=255)
                                                  gAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGAGGTTGCCAAGCGt            <  1:757372/69‑1 (MQ=255)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACtt                           >  1:2077093/1‑51 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACt                            >  1:186667/1‑51 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:2865100/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:1034355/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:850205/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:1260814/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:1305044/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:671686/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:548395/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:523283/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:385552/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:377486/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:1905284/1‑52 (MQ=39)
                                                    aTTTCGCTCCCGCCAGCGCATTTTTCAGATAATCCAGACCGTAGTTAAACTc                           >  1:2492044/1‑52 (MQ=39)
                                                      ttCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCg        >  1:3261968/1‑69 (MQ=255)
                                                      ttCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCg        >  1:3374398/1‑69 (MQ=255)
                                                      ttCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCg        >  1:309382/1‑69 (MQ=255)
                                                      ttCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCg        >  1:2669869/1‑69 (MQ=255)
                                                      ttCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCg        >  1:2030957/1‑69 (MQ=255)
                                                           tCCTGCGAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:3070933/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:2556839/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:814191/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:1046551/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:2225422/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:1347055/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:1619617/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:2429053/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:1894224/68‑1 (MQ=255)
                                                           tCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:1820988/68‑1 (MQ=255)
                                                           tCCTGACAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:161665/68‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:1459735/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:2345146/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:2906317/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:2623682/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:2205241/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:308268/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:438519/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:1576202/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:2412847/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:1836710/69‑1 (MQ=255)
                                                            ccTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:3132839/69‑1 (MQ=255)
                                                             cTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGta  <  1:673243/68‑1 (MQ=255)
                                                              tGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACgg    <  1:497550/65‑1 (MQ=255)
                                                              |                                                                   
AACATTGGCTGTTTGGTTCTGGCGTCAATGACGTTGGCATTTACATAAGGGAATTTCGCTCCTGCCAGCGCATTTTTCAGGTAATCCAGACCGTAGTTAAACTCGTGGTTGCCAAGCGTTCCGACGGTA  >  minE/2884153‑2884281

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: