Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,355,004 A→G 32.8% N64N (AAT→AAC dsbA ← periplasmic protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,355,0040AG32.8% 104.7 / 71.9 132N64N (AAT→AACdsbAperiplasmic protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base A (35/53);  new base G (7/36);  total (42/90)
Fisher's exact test for biased strand distribution p-value = 9.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

CACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGAAA  >  minE/2354940‑2355064
                                                                |                                                            
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:1880703/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:2892560/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:1897659/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:1913320/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:2494500/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:864150/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACAtt                                                            >  1:943559/1‑67 (MQ=255)
caccacCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTAACAtt                                                            >  1:2471693/1‑67 (MQ=255)
             gTTGACGTTGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:128768/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCa                                              >  1:78522/1‑68 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAg                                             >  1:2552337/1‑69 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAg                                             >  1:3338475/1‑69 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAg                                             >  1:233387/1‑69 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1460274/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1649440/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:3391385/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1839665/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1842479/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1589400/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:3333868/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:2397002/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:3301894/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:1019712/69‑1 (MQ=255)
             gTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAg                                             <  1:3059994/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACTTTATCAGAAATATGCAGAACttcttc                                    <  1:1654369/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:2186863/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:1008339/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:558745/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:2818165/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:2839914/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:501433/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:911244/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:1485681/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:1187207/69‑1 (MQ=255)
                      gTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:408384/69‑1 (MQ=255)
                       tATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:3023795/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:3211633/67‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:1842470/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:312978/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:630554/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:1963623/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:1550829/68‑1 (MQ=255)
                        aTTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   <  1:1324469/68‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:368026/1‑65 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:2907717/69‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               >  1:1059162/1‑69 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:1377342/69‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:3166756/69‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               >  1:2534015/1‑69 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:3386006/69‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:193459/69‑1 (MQ=255)
                           tAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               >  1:2030234/1‑69 (MQ=255)
                            aGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:340220/68‑1 (MQ=255)
                            aGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACt                               <  1:1496972/68‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:537202/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:1626957/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:519911/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:1447454/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:2287927/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:3254870/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:273701/69‑1 (MQ=255)
                                aTCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAg                          <  1:3387410/69‑1 (MQ=255)
                                  cTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2169569/68‑1 (MQ=255)
                                  cTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2395682/68‑1 (MQ=255)
                                  cTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:1334913/68‑1 (MQ=255)
                                  cTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2545771/68‑1 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATGAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:219449/1‑69 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:488971/1‑57 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:499455/1‑57 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:2894024/1‑57 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:2482710/1‑57 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCa                                   >  1:694012/1‑57 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:2212786/1‑69 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:3356334/1‑69 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:139052/1‑69 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:45033/1‑69 (MQ=255)
                                   ttCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAg                       >  1:943684/1‑69 (MQ=255)
                                          cTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTGAAACTGATAGCAGTGCg                   >  1:686654/1‑66 (MQ=255)
                                          cTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCg                   >  1:2756128/1‑66 (MQ=255)
                                              cGGCAGTTTTTTCTTCACGTTATCAGAAATATGCAGAACttcttc                                    <  1:1230644/45‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:657001/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:785771/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:910729/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:953957/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2945520/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2869843/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2768363/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:2600096/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:1957222/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:1624850/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:1430858/48‑1 (MQ=255)
                                                      ttttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:101273/48‑1 (MQ=255)
                                                       tttCTTCACGTTATCAGAAATATGCAGAACTTCTTCAAACTGATAGc                         <  1:1403051/47‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:2951521/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1890630/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1863574/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1859678/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1762941/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1759056/68‑1 (MQ=255)
                                                         tcttcACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGaaa  <  1:1154492/68‑1 (MQ=255)
                                                                |                                                            
CACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGAAATATGCAGAACTTCTTCAAACTGATAGCAGTGCGGGCAGAAGAAAGAGAAA  >  minE/2354940‑2355064

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: