Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 851,033 A→G 26.2% T56T (ACT→ACC yciA ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE851,0330AG26.2% 111.6 / 38.9 99T56T (ACT→ACCyciApredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (27/46);  new base G (23/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 102 total aligned reads displayed.

CGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACA  >  minE/850965‑851097
                                                                    |                                                                
cGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGcaca                                                                  <  1:2688002/69‑1 (MQ=255)
cGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGcaca                                                                  <  1:3099488/69‑1 (MQ=255)
cGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGcaca                                                                  <  1:1949463/69‑1 (MQ=255)
  tGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:2727414/69‑1 (MQ=255)
  tGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:498727/69‑1 (MQ=255)
  tGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:144221/69‑1 (MQ=255)
  tGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:1906074/69‑1 (MQ=255)
  tGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:1583689/69‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCTCACAGt                                                                <  1:441426/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:3302817/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:1503848/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:556885/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:3255261/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:1526783/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:731620/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:1180266/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:2251821/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:2431586/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:837487/63‑1 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTa                                                           >  1:1894939/1‑68 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:337072/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:917955/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:1715348/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:1906415/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:2149590/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:3012439/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:2965182/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:2917679/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:373536/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:1158978/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:1330396/1‑69 (MQ=255)
        gcagcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:1211419/1‑69 (MQ=255)
        gcaacaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAc                                                          >  1:3277004/1‑69 (MQ=255)
          agcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcgc                                                        >  1:2731749/1‑69 (MQ=255)
          agcaCACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTAcgc                                                        >  1:1151886/1‑69 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCAGAGTCAc                                                             <  1:442485/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                     <  1:1120890/65‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:2693020/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:2476086/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:1200088/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:3053991/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:3082247/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:2082228/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:2048645/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:181996/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:3160168/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:191539/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:623825/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:188471/57‑1 (MQ=255)
                 cacaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:1428199/57‑1 (MQ=255)
                  acaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                     <  1:1493995/64‑1 (MQ=255)
                  acaTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACGGTCACTACACGAc                                                     <  1:2471109/64‑1 (MQ=255)
                    aTCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:2787092/51‑1 (MQ=255)
                     tCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGt                                                                <  1:652975/50‑1 (MQ=255)
                         cGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCAc                                                             <  1:116752/49‑1 (MQ=255)
                                   ctgCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:3422524/67‑1 (MQ=255)
                                   cGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:1417919/69‑1 (MQ=255)
                                   cGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:3100567/69‑1 (MQ=255)
                                   cGGCCGTAAGAAAGTCATTCCTTCAACCCCCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:1065536/69‑1 (MQ=255)
                                   cGGCCGTAAGAAAGTCATTCCTACAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:1058667/69‑1 (MQ=255)
                                   cGGCCGTAAGAAAGTCATTCCATCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               <  1:2242693/69‑1 (MQ=255)
                                    ggCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               >  1:1779479/1‑68 (MQ=255)
                                    ggCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAg                               >  1:1065227/1‑68 (MQ=255)
                                      ccGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAAt                            <  1:196863/69‑1 (MQ=255)
                                                 tCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCaa                 >  1:2589372/1‑69 (MQ=255)
                                                 tCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCaa                 >  1:2163290/1‑69 (MQ=255)
                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                >  1:787755/1‑69 (MQ=255)
                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                >  1:2445577/1‑69 (MQ=255)
                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                >  1:415248/1‑69 (MQ=255)
                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                >  1:3182290/1‑69 (MQ=255)
                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAat                >  1:2019077/1‑69 (MQ=255)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATgggg                                              >  1:721540/1‑35 (MQ=25)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:553606/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1633796/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1605331/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1535841/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1507194/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:3269759/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1434385/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:136036/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1266828/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1249013/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:778052/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1008285/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1793633/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1904338/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:1976684/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:3016844/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:3006899/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:3000926/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:2379315/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:2922481/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:247747/1‑36 (MQ=37)
                                                     tCCTTCAACCCGCACGGTCACTACACGACCATGGGc                                              >  1:253272/1‑36 (MQ=37)
                                                      ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATcc            >  1:1897177/1‑69 (MQ=255)
                                                      ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATcc            >  1:2849274/1‑69 (MQ=255)
                                                                gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa  <  1:30582/69‑1 (MQ=255)
                                                                gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa  <  1:1490792/69‑1 (MQ=255)
                                                                gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa  <  1:218347/69‑1 (MQ=255)
                                                                gCACAGTCACAACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACa  <  1:2878724/69‑1 (MQ=255)
                                                                    |                                                                
CGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACA  >  minE/850965‑851097

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