Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 387,114 A→C 7.4% L12F (TTA→TTC) ‡ ybdD → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE387,1140AC 7.4% 236.7 / 12.2 94L12F (TTA→TTC) ‡ybdDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (52/35);  new base C (5/2);  total (57/37)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAAC  >  minE/387050‑387183
                                                                |                                                                     
tttCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAg                                                                      >  1:510136/1‑66 (MQ=255)
  tCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAgg                                                                <  1:436353/70‑1 (MQ=255)
  tCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAgg                                                                <  1:4838640/70‑1 (MQ=255)
   cAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcg                                                              >  1:5618781/1‑71 (MQ=255)
       aaCATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGgcggcg                                                           >  1:4153458/1‑70 (MQ=255)
           tGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                                                      <  1:2848884/71‑1 (MQ=255)
                aGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatg                                                  <  1:4917597/70‑1 (MQ=255)
                 gCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTgatgat                                                <  1:1957616/71‑1 (MQ=255)
                    gAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCa                                                                  >  1:215504/1‑50 (MQ=255)
                    gAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                                 >  1:5376051/1‑67 (MQ=255)
                    gAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:1631052/71‑1 (MQ=255)
                    gAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:2177515/71‑1 (MQ=255)
                    gAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:1790062/71‑1 (MQ=255)
                     aGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtg                                                     <  1:1440691/62‑1 (MQ=255)
                      ggCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTgatgat                                                >  1:4865496/1‑66 (MQ=255)
                       gCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcggc                                                            <  1:1628028/53‑1 (MQ=255)
                       gCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTg                                              <  1:1261203/67‑1 (MQ=255)
                         cAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGc                                        >  1:3990063/1‑71 (MQ=255)
                          aaGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                                                      <  1:4857845/56‑1 (MQ=255)
                           aGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTc                                                                   <  1:2949382/42‑1 (MQ=255)
                           aGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTc                                                                   <  1:675230/42‑1 (MQ=255)
                           aGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAg                                                        >  1:3806213/1‑53 (MQ=255)
                           aGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAg                                                        >  1:1387564/1‑53 (MQ=255)
                           aGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:172929/64‑1 (MQ=255)
                            gATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAg                                                                      <  1:3068182/38‑1 (MQ=255)
                             aTGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcgg                                                             >  1:5330562/1‑46 (MQ=255)
                             aTGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    >  1:1601391/1‑71 (MQ=255)
                             aTGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    >  1:4886846/1‑71 (MQ=255)
                              tGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAt                                   >  1:2949322/1‑71 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGa                                                          <  1:2905471/46‑1 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                                                      <  1:1866695/50‑1 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                                                      <  1:4828258/50‑1 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt                                            >  1:1954260/1‑60 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt                                            >  1:5027236/1‑60 (MQ=255)
                                tttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTa                                 <  1:1122434/71‑1 (MQ=255)
                                 ttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    >  1:2650542/1‑67 (MQ=255)
                                 ttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                                >  1:2445598/1‑71 (MQ=255)
                                 ttGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                                >  1:1661989/1‑71 (MQ=255)
                                  tGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                                >  1:1729895/1‑70 (MQ=255)
                                  tGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                                >  1:2071149/1‑70 (MQ=255)
                                  tGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc                                >  1:2450142/1‑70 (MQ=255)
                                      tCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTAt                                          <  1:5027329/56‑1 (MQ=255)
                                       cACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACaa                           >  1:2493443/1‑70 (MQ=255)
                                       cACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAAc                          >  1:4476256/1‑71 (MQ=255)
                                        aCTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                                                      >  1:2551315/1‑42 (MQ=255)
                                        aCTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAAc                          >  1:305563/1‑70 (MQ=255)
                                         cTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAAc                          <  1:598127/69‑1 (MQ=255)
                                           ggCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                                 >  1:1795185/1‑44 (MQ=255)
                                           ggCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTa                                           >  1:4115185/1‑50 (MQ=255)
                                           ggCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTAt                                          >  1:596595/1‑51 (MQ=255)
                                            gCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTg                                              <  1:27387/46‑1 (MQ=255)
                                             cAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGt                     >  1:169441/1‑70 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                                 >  1:2336443/1‑41 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:981279/45‑1 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:2746943/45‑1 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTgg                                             <  1:2263016/45‑1 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAt                                   >  1:5572510/1‑55 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAt                                   >  1:5565807/1‑55 (MQ=255)
                                              aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTc                    >  1:4569649/1‑70 (MQ=255)
                                                aaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtg                                                     <  1:3806854/35‑1 (MQ=255)
                                                 aGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGaa                 >  1:3932333/1‑70 (MQ=255)
                                                 aGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAAc                <  1:2296528/71‑1 (MQ=255)
                                                   ccGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                                                 >  1:4265816/1‑36 (MQ=255)
                                                   ccGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACa               >  1:4169278/1‑70 (MQ=255)
                                                    cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    <  1:3742088/48‑1 (MQ=255)
                                                    cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    <  1:4007465/48‑1 (MQ=255)
                                                    cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGa                                    <  1:3234295/48‑1 (MQ=255)
                                                    cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACaa                           <  1:4582039/57‑1 (MQ=255)
                                                    aaaaaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCata             >  1:1368742/4‑71 (MQ=255)
                                                    aaaaaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCata             >  1:2503899/4‑71 (MQ=255)
                                                    aaaaaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCata             >  1:1375512/4‑71 (MQ=255)
                                                      gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACt                         >  1:1563660/1‑57 (MQ=255)
                                                      gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACt                         >  1:4746026/1‑57 (MQ=255)
                                                      gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACt                         >  1:3960817/1‑57 (MQ=255)
                                                      gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACt                         >  1:2710601/1‑57 (MQ=255)
                                                      aaaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCATATg           <  1:2994722/70‑1 (MQ=255)
                                                      aaaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCATATg           <  1:784947/70‑1 (MQ=255)
                                                       aaaaTATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAt                                   >  1:1798256/1‑46 (MQ=255)
                                                       aaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTAc                                >  1:307246/1‑49 (MQ=255)
                                                       aaaaTATCTCGGCCAGGCGGCGAAGATGTTGATTGGCATTCCAGACTACGACAACTACGTCGAGCATATg           >  1:2617577/1‑70 (MQ=255)
                                                        aaaTATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCg         >  1:331928/1‑71 (MQ=255)
                                                        aaaTATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCg         >  1:405943/1‑71 (MQ=255)
                                                        aaaTATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCg         <  1:3826447/71‑1 (MQ=255)
                                                          atatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCggg       >  1:5142009/1‑71 (MQ=255)
                                                            atTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGt      >  1:2606716/1‑70 (MQ=255)
                                                            atTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGt      >  1:3267710/1‑70 (MQ=255)
                                                            atTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGt      >  1:4408503/1‑70 (MQ=255)
                                                            atTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGt      >  1:1579984/1‑70 (MQ=255)
                                                            atTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGt      >  1:3907523/1‑70 (MQ=255)
                                                             tttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTa    <  1:4862363/71‑1 (MQ=255)
                                                             tttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTa    <  1:3716080/71‑1 (MQ=255)
                                                              ttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTaa   >  1:2608251/1‑71 (MQ=255)
                                                              ttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTaa   >  1:967722/1‑71 (MQ=255)
                                                               tAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAAc  >  1:4330098/1‑71 (MQ=255)
                                                                aGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAAc  <  1:342600/70‑1 (MQ=255)
                                                                |                                                                     
TTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAAC  >  minE/387050‑387183

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: