Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,942,643 T→C 14.4% K83K (AAA→AAG yqfB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,942,6430TC14.4% 179.8 / 15.7 103K83K (AAA→AAGyqfBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (49/39);  new base C (12/3);  total (61/42)
Fisher's exact test for biased strand distribution p-value = 9.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 104 total aligned reads displayed.

TTTAAAGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGTTA  >  minE/1942579‑1942710
                                                                |                                                                   
tttAAAGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACtttt                                                                  <  1:1249832/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:1792231/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:1541368/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:2795182/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:1882510/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:2387722/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:3298531/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:1506308/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:1479061/68‑1 (MQ=255)
     aGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:2175421/68‑1 (MQ=255)
      gACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTttcag                                                             <  1:3126542/67‑1 (MQ=255)
      gACATTTAAATTCAATCACATAAAACTGTGTCTGGTCGGGATAGATGTCGGCAATGACCTTTttcagtt                                                           <  1:3386627/69‑1 (MQ=255)
      gACATTTAAATTCAATCACATAAAACTGTGTCTGGTCGGGATAGATGTCGGCAATGACCTTTttcagtt                                                           <  1:1514975/69‑1 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTtcagtcag                                                    >  1:2193767/1‑68 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:1420045/1‑69 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:63528/1‑69 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:689309/1‑69 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:3394399/1‑69 (MQ=255)
              aaTTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:1927596/1‑69 (MQ=255)
                ttCAATCACATATAGTTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:3195485/1‑67 (MQ=255)
                ttCAATCACATATAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:958335/1‑67 (MQ=255)
                ttCAATCACATATAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:1549728/1‑67 (MQ=255)
                ttCAATCACATATAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAgg                                                   >  1:3301113/1‑67 (MQ=255)
                  cAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCa                                               >  1:2619222/1‑69 (MQ=255)
                      cacaTAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCatat                                            >  1:2628550/1‑68 (MQ=255)
                      cacaTAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATAtt                                           >  1:2435166/1‑69 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:1620968/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:457908/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:616570/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:3388366/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:2965099/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:1283367/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:1084370/68‑1 (MQ=255)
                                      tGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGt                           >  1:670436/1‑69 (MQ=255)
                                      tGACCGGGATAGATATCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATg                            <  1:1786458/68‑1 (MQ=255)
                                           gggATAGATGTCATCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:3324931/69‑1 (MQ=255)
                                           gggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:465570/69‑1 (MQ=255)
                                           gggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:2981984/69‑1 (MQ=255)
                                           gggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:1508691/69‑1 (MQ=255)
                                           gggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:1278576/69‑1 (MQ=255)
                                            ggATAGATGTCAGCAATGTCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     <  1:1528070/69‑1 (MQ=255)
                                            ggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     <  1:461091/69‑1 (MQ=255)
                                            ggATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtc                      <  1:1265022/68‑1 (MQ=255)
                                              aTAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     >  1:3423588/1‑67 (MQ=255)
                                              aTAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     >  1:1717803/1‑67 (MQ=255)
                                              aTAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     >  1:2943087/1‑67 (MQ=255)
                                              aTAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtca                     >  1:710716/1‑67 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:108521/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:565135/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:62469/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:2875010/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:910385/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:2823341/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:1481381/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:1511566/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:2524396/1‑69 (MQ=255)
                                                 gATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTtcagtcag                >  1:1594811/1‑69 (MQ=255)
                                                   tGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              <  1:1281423/69‑1 (MQ=255)
                                                   tGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              <  1:2795776/69‑1 (MQ=255)
                                                   tGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              <  1:667059/69‑1 (MQ=255)
                                                   tGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              <  1:56876/69‑1 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCt                                      >  1:2513718/1‑37 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCt                                      >  1:2996609/1‑37 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:1162889/1‑38 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:67667/1‑38 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:457319/1‑38 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:2474063/1‑38 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:3115314/1‑38 (MQ=255)
                                                           aTGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTg                                     >  1:649817/1‑38 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:583850/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:587639/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:602923/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:666528/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:1357624/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:335440/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:3158263/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:2584427/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:2905159/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGTATCCAGc      >  1:924570/1‑69 (MQ=255)
                                                           aTGACCTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCCGTCAGCGAATCCAg       >  1:862669/1‑68 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:1972966/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:1379041/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:869539/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:2459049/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:1229624/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:2901579/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:144539/1‑58 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:674325/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:1846951/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:389700/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:2080481/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:2084510/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:2419181/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:2495608/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:2758832/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCAGCATGTTTTTCAGTCAGCGTATCCAGCGt    <  1:1534088/68‑1 (MQ=255)
                                                              aCTTTTTTCAGTTCAGTCAGGGTCATATTTACCTGCTCTGCATGTTTTTCAGTCAGCg              >  1:2870990/1‑58 (MQ=255)
                                                               cTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGTTa  >  1:1758065/1‑69 (MQ=255)
                                                               cTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGTTa  >  1:1887669/1‑69 (MQ=255)
                                                               cTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGTTa  >  1:2161481/1‑69 (MQ=255)
                                                                |                                                                   
TTTAAAGACATTTAAATTCAATCACATAAAATTGTGTCTGACCGGGATAGATGTCAGCAATGACTTTTTTCAGTTCAGTCAGGGTCATATTTTCCTGCTCTGCATGTTTTTCAGTCAGCGTATCCAGCGTTA  >  minE/1942579‑1942710

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: