Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 197,781 T→C 31.4% L213L (CTA→CTG yaeB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE197,7810TC31.4% 96.7 / 56.2 111L213L (CTA→CTGyaeBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/60);  new base C (26/9);  total (42/69)
Fisher's exact test for biased strand distribution p-value = 1.41e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTG  >  minE/197713‑197843
                                                                    |                                                              
tAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCTAACGTTAAAATCATGt                                                                <  1:2525548/69‑1 (MQ=255)
tAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGt                                                                <  1:1987976/69‑1 (MQ=255)
tAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGt                                                                <  1:498527/69‑1 (MQ=255)
  aCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTa                                                               <  1:3192485/68‑1 (MQ=255)
  aCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTa                                                               <  1:802249/68‑1 (MQ=255)
  aCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTa                                                               <  1:1292755/68‑1 (MQ=255)
    gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGc                                                             <  1:3120160/68‑1 (MQ=255)
    gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGc                                                             <  1:2028597/68‑1 (MQ=255)
    gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGc                                                             <  1:1160068/68‑1 (MQ=255)
    gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGc                                                             <  1:1178177/68‑1 (MQ=255)
    gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGc                                                             <  1:1287406/68‑1 (MQ=255)
       gTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGa                                                         >  1:3084774/1‑69 (MQ=255)
          ccAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa                                                           >  1:51386/1‑64 (MQ=255)
          ccAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACgg                                                      >  1:2335505/1‑69 (MQ=255)
          ccAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACgg                                                      >  1:526262/1‑69 (MQ=255)
          ccAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACgg                                                      >  1:1101660/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:38741/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:2143731/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:3216185/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:530204/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:2025784/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:1063684/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:537904/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:607216/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:1846094/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:1290502/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:9362/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:1135056/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGACAGACGGCATAGGttt                                             >  1:2206837/1‑69 (MQ=255)
                   aGACTTCAAAACCTGCGTCGGTGACGCGCCAGAGAACGTTAAAATCATGCAGCCAGACGGCATAGGttt                                             >  1:1458729/1‑69 (MQ=255)
                   aGACTTAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGtttt                                            >  1:2275266/1‑69 (MQ=255)
                       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:100808/69‑1 (MQ=255)
                       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:409108/69‑1 (MQ=255)
                       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:10689/69‑1 (MQ=255)
                       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:552313/69‑1 (MQ=255)
                       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:568055/69‑1 (MQ=255)
                        tCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:2836772/68‑1 (MQ=255)
                        tCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:1256891/68‑1 (MQ=255)
                         cAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa                                                           >  1:2667504/1‑49 (MQ=255)
                         cAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa                                                           >  1:2212887/1‑49 (MQ=255)
                         cAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa                                                           >  1:2663428/1‑49 (MQ=255)
                         cAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCa                                                           >  1:2499487/1‑49 (MQ=255)
                            aaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGtt                                     <  1:1663552/68‑1 (MQ=255)
                            aaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGtt                                     <  1:3212755/68‑1 (MQ=255)
                            aaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGtt                                     <  1:1132623/68‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:2632406/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:256240/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:3075612/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:2239556/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:630124/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:614788/69‑1 (MQ=255)
                               cTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:1844073/69‑1 (MQ=255)
                                tGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTc                                <  1:2428841/69‑1 (MQ=255)
                                 gCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTc                                <  1:2513220/68‑1 (MQ=255)
                                 gCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTc                                <  1:2464559/68‑1 (MQ=255)
                                 gCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTc                                <  1:2059048/68‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGGCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:607588/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTTCCCGTTTCCTCACCTTTAc                        <  1:2912851/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:867869/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:3303016/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:774595/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:951402/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:556109/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:3032893/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:1819877/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:2253271/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:2375047/69‑1 (MQ=255)
                                        tGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:2591097/69‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:2810107/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:2814051/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:1934767/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:3192275/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:2903600/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:2136920/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:806332/51‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:909396/68‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAc                        <  1:1942780/68‑1 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTAACCTTTAc                        >  1:1746541/1‑68 (MQ=255)
                                         gACGCGCCAGCGAACGTTAAAATCATGTAACCAGACGGCATAGGTTTTGcc                                         <  1:2653384/51‑1 (MQ=255)
                                            gcgcCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGcc                                         <  1:731675/48‑1 (MQ=255)
                                            gcgcCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATa                    >  1:918881/1‑69 (MQ=255)
                                            gcgcCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATa                    >  1:1392475/1‑69 (MQ=255)
                                            gcgcCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATa                    >  1:2482588/1‑69 (MQ=255)
                                                  gCGAACGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCt                                 <  1:460744/50‑1 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTAGGATAGGCCGGGcgcg         >  1:350308/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:1568453/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:1689340/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:1276323/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:1882343/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:1910092/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:2086863/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:512903/1‑69 (MQ=255)
                                                       cGTTAAAATCATGCAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGcgcg         >  1:2707021/1‑69 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:689142/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:556258/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:1681537/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:1936996/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:3266615/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:1306292/69‑1 (MQ=255)
                                                              aTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg  <  1:291260/69‑1 (MQ=255)
                                                                    |                                                              
TAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTG  >  minE/197713‑197843

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: