Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,599,135 A→T 18.8% I33I (ATT→ATA gor ← glutathione oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,599,1350AT18.8% 202.7 / 39.9 85I33I (ATT→ATAgorglutathione oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/45);  new base T (4/12);  total (28/57)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CCGCGTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCTG  >  minE/2599065‑2599195
                                                                      |                                                            
ccgcGTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCa                                                              <  1:2738527/71‑1 (MQ=255)
ccgcGTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCa                                                              <  1:1395052/71‑1 (MQ=255)
ccgcGTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCa                                                              <  1:2314039/71‑1 (MQ=255)
  gcgTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAAt                                                            <  1:901072/71‑1 (MQ=255)
       ccACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAgcg                                                       <  1:4095330/71‑1 (MQ=255)
       ccACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAgcg                                                       <  1:3543614/71‑1 (MQ=255)
       ccACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAgcg                                                       <  1:2652588/71‑1 (MQ=255)
       ccACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAgcg                                                       <  1:2080420/71‑1 (MQ=255)
       ccACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAgcg                                                       >  1:1727680/1‑71 (MQ=255)
           aTCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCaa                                                             <  1:96681/61‑1 (MQ=255)
           aTCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCaa                                                             <  1:1079573/61‑1 (MQ=255)
            tCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACAt                                                  <  1:3685821/71‑1 (MQ=255)
            tCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACAt                                                  >  1:705794/1‑71 (MQ=255)
            tCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACAt                                                  >  1:1403718/1‑71 (MQ=255)
              aCTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACAttt                                                >  1:470157/1‑71 (MQ=255)
              aCTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACAttt                                                <  1:3650821/71‑1 (MQ=255)
              aCTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACAtt                                                 <  1:2180158/70‑1 (MQ=255)
               cTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCTCACATTTc                                               <  1:1954759/71‑1 (MQ=255)
               cTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTc                                               <  1:2181304/71‑1 (MQ=255)
               cTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTc                                               <  1:2876506/71‑1 (MQ=255)
               cTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTc                                               <  1:3188744/71‑1 (MQ=255)
                   tttCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTCTgg                                            <  1:3977316/70‑1 (MQ=255)
                       ggCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTa                                       <  1:334251/71‑1 (MQ=255)
                          acacaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTCTg                                             >  1:3311400/1‑62 (MQ=255)
                          acacaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTCTg                                             >  1:2849943/1‑62 (MQ=255)
                          acacaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTAc                                      >  1:1731373/1‑69 (MQ=255)
                          acacaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACa                                     >  1:2703767/1‑70 (MQ=255)
                            acaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTg                                             <  1:2570656/60‑1 (MQ=255)
                            acaGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTg                                             <  1:3009812/60‑1 (MQ=255)
                                ccAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAg                                                         >  1:2157777/1‑44 (MQ=255)
                                ccAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTCTGGCCGTACATAGCcgc                              >  1:3144049/1‑71 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:114861/45‑1 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:3584464/45‑1 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:2165391/45‑1 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:2121764/45‑1 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:142370/45‑1 (MQ=255)
                                   aCATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGca                                                     <  1:1946107/45‑1 (MQ=255)
                                    cATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTa                                       <  1:4031214/58‑1 (MQ=255)
                                    cATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTa                                       <  1:1850779/58‑1 (MQ=255)
                                    cATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTa                                       <  1:3524188/58‑1 (MQ=255)
                                    cATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTa                                       <  1:1823428/58‑1 (MQ=255)
                                     aTTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGt                         >  1:3908558/1‑71 (MQ=255)
                                     aTTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGt                         <  1:206565/71‑1 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCAGATTTCTGGCCGTACATAGCCGCGCGGtt                        <  1:381711/71‑1 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGc                                           >  1:3623706/1‑52 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGc                                           >  1:1315538/1‑52 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGtt                        <  1:234713/71‑1 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGtt                        <  1:2773174/71‑1 (MQ=255)
                                      tttACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGtt                        <  1:1735692/71‑1 (MQ=255)
                                        tACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCa                                                          >  1:253875/1‑35 (MQ=255)
                                        tACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCa                                                          >  1:1089223/1‑35 (MQ=255)
                                         aCGCAGGTGCCGCCCAGCTCTTTGGCTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGa                      <  1:1561394/70‑1 (MQ=255)
                                          cGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATg                    <  1:819263/71‑1 (MQ=255)
                                               gTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGcc                                <  1:2667428/54‑1 (MQ=255)
                                                tGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTgg                                            >  1:2098388/1‑41 (MQ=255)
                                                tGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGcc                                <  1:2672285/53‑1 (MQ=255)
                                                tGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGc               >  1:2309523/1‑70 (MQ=255)
                                                tGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGc               >  1:887465/1‑70 (MQ=255)
                                                       cAGCTCTTTGGCTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:1656001/71‑1 (MQ=255)
                                                       cAGCTCTTTGGCTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:3226304/71‑1 (MQ=255)
                                                       cAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCcgcgc                            <  1:2323125/50‑1 (MQ=255)
                                                       cAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGa                      >  1:1186543/1‑56 (MQ=255)
                                                       cAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       >  1:3532713/1‑71 (MQ=255)
                                                       cAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccg        >  1:1644795/1‑70 (MQ=255)
                                                       cAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccg        >  1:2914401/1‑70 (MQ=255)
                                                       aaGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACa                                     >  1:1781103/2‑41 (MQ=255)
                                                        aGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgcc      <  1:2788737/71‑1 (MQ=255)
                                                         tctctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGtt                        <  1:2800536/51‑1 (MQ=255)
                                                           tctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCt   >  1:3783434/1‑71 (MQ=255)
                                                           tctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCt   >  1:1813287/1‑71 (MQ=255)
                                                           tctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCt   >  1:2903417/1‑71 (MQ=255)
                                                           tctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCt   >  1:3551364/1‑71 (MQ=255)
                                                           tctTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCt   >  1:2668018/1‑71 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:2280714/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:3844730/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:3877556/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:4115902/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:842220/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATa                                   <  1:632931/37‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:1551895/65‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:2873848/65‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:3283320/65‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:1509760/65‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:4121025/65‑1 (MQ=255)
                                                             tttGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgc       <  1:3493507/65‑1 (MQ=255)
                                                               tGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCTg  <  1:2505202/68‑1 (MQ=255)
                                                                  cTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgcc      <  1:719078/61‑1 (MQ=255)
                                                                  cTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgcc      <  1:1293314/61‑1 (MQ=255)
                                                                  cTTCTATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgcc      <  1:1275494/61‑1 (MQ=255)
                                                                    tCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATAccgcc      <  1:3345376/59‑1 (MQ=255)
                                                                      |                                                            
CCGCGTGCCACATCACTTTTTTCGGCACACAGCCAACATTTACGCAGGTGCCGCCCAGCTCTTTGGCTTCAATCAGCGCACATTTCTGGCCGTACATAGCCGCGCGGTTGATGGAGGCGATACCGCCGCTG  >  minE/2599065‑2599195

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: