Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,526,156 A→G 6.6% V150V (GTA→GTG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,526,1560AG 6.6% 362.2 / 11.7 121V150V (GTA→GTGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (52/61);  new base G (2/6);  total (54/67)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 121 total aligned reads displayed.

ACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAT  >  minE/2526088‑2526225
                                                                    |                                                                     
aCCGCCGGAACGGCGATTCGCGAGTCGATGGAAATTATTCAGGCCAATGGCGCGACGCTTGCAGGCGTGtt                                                                     <  1:4545109/71‑1 (MQ=255)
aCCGCCGGAACGGCGATTCGCGAGTCGATGGAAATTATTCAGGCCAATGGCGCGACGCTTGCAGGCGTGtt                                                                     <  1:4637294/71‑1 (MQ=255)
aCCGCCGGAACGGCGATTCGCGAGTCGATGGAAATTATTCAGGCCAATGGCGCGACGCTTGCAGGCGTGtt                                                                     <  1:1948679/71‑1 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                                   <  1:4552598/71‑1 (MQ=255)
   gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                                  >  1:1386170/1‑71 (MQ=255)
     cgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAttt                                                                <  1:2995334/71‑1 (MQ=255)
      ggAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAttt                                                                <  1:1167571/70‑1 (MQ=255)
      ggAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                               >  1:1567972/1‑71 (MQ=255)
       gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                                              <  1:3602513/71‑1 (MQ=255)
       gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                                              <  1:2540319/71‑1 (MQ=255)
       gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                                              <  1:509088/71‑1 (MQ=255)
         aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                                             <  1:3364586/70‑1 (MQ=255)
         aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                                            >  1:3189930/1‑71 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           >  1:2655230/1‑71 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           >  1:1717276/1‑71 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                                            >  1:3802852/1‑70 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                         <  1:4708359/71‑1 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                         <  1:2359869/71‑1 (MQ=255)
               aTTCGCGAGTCGATGGAAATTATTCAGGCCAATGGCGCGACGCTTGCAGGCGTGtt                                                                     <  1:931789/56‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     <  1:1813623/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     <  1:4430620/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     <  1:3203008/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     <  1:4433346/71‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                       <  1:930754/68‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                       <  1:3444885/68‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                       <  1:1388861/68‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                       <  1:3668100/68‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                       <  1:4612966/68‑1 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     >  1:3989636/1‑69 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     >  1:3101991/1‑69 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                     >  1:2176534/1‑69 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                    <  1:333901/70‑1 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                    <  1:3382243/70‑1 (MQ=255)
                      aaTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGaa                                               <  1:1515140/71‑1 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCTTATTGATTTCGCTTGAt                                                        >  1:2449045/1‑61 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAttt                                                                <  1:1466031/53‑1 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                        >  1:5095759/1‑61 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                        >  1:3872310/1‑61 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                        >  1:3819201/1‑61 (MQ=255)
                       aTCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                        >  1:4755970/1‑61 (MQ=255)
                       aTCGATGGAGATTAATCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                        >  1:418557/1‑61 (MQ=255)
                              gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcg                                       <  1:1439874/71‑1 (MQ=255)
                               aaattattCAGGCCAATGGCGCGACGCTTGCAGGCGTGtt                                                                     <  1:1927061/38‑1 (MQ=25)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:3490531/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:116802/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:4296429/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:4931546/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:3146338/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:267229/48‑1 (MQ=255)
                                 attattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                           <  1:2330446/48‑1 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCATGAACGGGGGCgcggc                                   >  1:2736331/1‑70 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                   >  1:835618/1‑70 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                   >  1:3737231/1‑70 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                   >  1:4687147/1‑70 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                   >  1:652822/1‑70 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:4071926/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:602590/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:756836/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:4924159/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:3393762/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:3763731/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                      <  1:332421/49‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                                 <  1:5112914/70‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                                 <  1:811052/70‑1 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:1879797/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:1949338/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:4883271/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:5148637/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:4701396/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:1540140/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:775077/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:4059623/1‑52 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                                  >  1:3924447/1‑52 (MQ=255)
                                          gCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                            <  1:327484/70‑1 (MQ=255)
                                             aaTGGCGCGACGCTTGCTGGCGTATTGATTTCTCTTGATCg                                                      <  1:2169293/41‑1 (MQ=255)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                        <  1:3679602/70‑1 (MQ=255)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGGCAGGAACGGGGGCGCGGCGAGATTTCGGc                        <  1:4174961/70‑1 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                 <  1:3025023/44‑1 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                 <  1:4618617/44‑1 (MQ=255)
                                                ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                 <  1:2956479/43‑1 (MQ=255)
                                                 gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                    >  1:2791372/1‑71 (MQ=255)
                                                 gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                    >  1:1451258/1‑71 (MQ=255)
                                                 gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCACGAACGGGGGCGCGGCGAGATTTCGGc                        <  1:3272224/67‑1 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                 >  1:1282537/1‑70 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                 >  1:3720776/1‑70 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                 >  1:5036252/1‑70 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                 >  1:3007146/1‑70 (MQ=255)
                                                       cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                >  1:5035118/1‑69 (MQ=255)
                                                       cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGaa              <  1:1916898/71‑1 (MQ=255)
                                                           tGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                 >  1:938234/1‑64 (MQ=255)
                                                           tGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt           >  1:2653922/1‑70 (MQ=255)
                                                           tGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt           >  1:3522613/1‑70 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATGGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt           <  1:2525949/69‑1 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt           <  1:902630/69‑1 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa         >  1:1889858/1‑69 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa         >  1:1597250/1‑69 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa         >  1:2677946/1‑69 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa         >  1:1338175/1‑69 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa         >  1:1208277/1‑69 (MQ=255)
                                                                 cGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt     <  1:3255443/70‑1 (MQ=255)
                                                                    |                                                                     
ACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAT  >  minE/2526088‑2526225

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: