Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,080,002 G→A 14.1% P507L (CCA→CTA)  ydhV ← predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,080,0020GA14.1% 179.2 / 18.2 71P507L (CCA→CTA) ydhVpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (7/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GTTTTTGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGTTAACA  >  minE/1079932‑1080070
                                                                       |                                                                    
gTTTTTGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATg                                                                      <  1:2834698/71‑1 (MQ=255)
gTTTTTGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATg                                                                      >  1:2698887/1‑71 (MQ=255)
   tttGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGc                                                                    <  1:236892/70‑1 (MQ=255)
    ttGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCa                                                                   <  1:534044/70‑1 (MQ=255)
          tCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCgg                                                             >  1:291926/1‑70 (MQ=255)
          tCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCgg                                                             <  1:2506437/70‑1 (MQ=255)
          tCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCgg                                                             >  1:3844727/1‑70 (MQ=255)
             tCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGa                                                            <  1:3450454/68‑1 (MQ=255)
             tCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAAc                                                          <  1:1249713/70‑1 (MQ=255)
              cGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAAcc                                                         >  1:3604048/1‑70 (MQ=255)
                 aaGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTCAGTTTCAATGGCAAGCCGGAACCAATa                                                     >  1:2079257/1‑70 (MQ=255)
                  aGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATaa                                                    <  1:1274652/71‑1 (MQ=255)
                   gCATCTTCAGAACCAAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATAGCAAGCCGGAACCAATa                                                     >  1:2832763/1‑70 (MQ=255)
                     aTCTTCAGAACCAAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATAGCa                                                                   >  1:868232/1‑54 (MQ=38)
                     aTCTTCAGAACCAAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATAGCAAGCCGGAACCAATAAAg                                                  >  1:3587773/1‑71 (MQ=255)
                     aTCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGt                                                 <  1:3510901/71‑1 (MQ=255)
                     aTCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGt                                                 >  1:163169/1‑71 (MQ=255)
                        ttCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGa                                                            <  1:3523241/57‑1 (MQ=255)
                        ttCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGa                                              >  1:2049911/1‑71 (MQ=255)
                        ttCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGa                                              >  1:415837/1‑71 (MQ=255)
                          cAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATg                                                                      >  1:63639/1‑45 (MQ=255)
                           aGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGatat                                           >  1:2274679/1‑71 (MQ=255)
                             aaCCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATgg                                         >  1:1487156/1‑71 (MQ=255)
                             aaCCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATgg                                         >  1:3908691/1‑71 (MQ=255)
                               ccAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGt                                       <  1:1325850/71‑1 (MQ=255)
                                cAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCAGGAACCAATAAAGTTGATATGGGTa                                      <  1:3009436/71‑1 (MQ=255)
                                  aaaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTAt                                     <  1:67687/71‑1 (MQ=255)
                                   aataGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCTAGGCGGAACCAATAAAGTTGat                                             <  1:704822/59‑1 (MQ=255)
                                   aaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAg                                                                 >  1:474967/1‑42 (MQ=255)
                                   aaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAg                                                                 >  1:1802934/1‑42 (MQ=255)
                                   aaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGaa                                                           >  1:909804/1‑48 (MQ=255)
                                   aaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGaa                                                           >  1:3284312/1‑48 (MQ=255)
                                   aaaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGat                                             <  1:908320/62‑1 (MQ=255)
                                    aaaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGc                                   >  1:458639/1‑71 (MQ=255)
                                     aaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAg                                                  >  1:2944297/1‑55 (MQ=255)
                                     aaGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAg                                                  >  1:768478/1‑55 (MQ=255)
                                      aGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGt                                 >  1:2891842/1‑71 (MQ=255)
                                      aGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGt                                 >  1:1196908/1‑71 (MQ=255)
                                      aGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGt                                 >  1:2257651/1‑71 (MQ=255)
                                      aGTTCTTTCGCCACTTCACGTTGCAGTTTCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGt                                       <  1:3657943/65‑1 (MQ=255)
                                        ttCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAg                                                                 <  1:2577760/37‑1 (MQ=255)
                                        ttCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGatat                                           <  1:2149765/59‑1 (MQ=255)
                                         tCTTTCGCCACTTCACGTTGCAGTTTCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCAt                              >  1:2499112/1‑71 (MQ=255)
                                          cTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGtt                                                >  1:2848415/1‑52 (MQ=255)
                                           tttCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCAt                              >  1:3771057/1‑69 (MQ=255)
                                            ttCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATg                             >  1:2823591/1‑69 (MQ=255)
                                               gCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGat                                             <  1:3869691/50‑1 (MQ=255)
                                               gCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGat                                             <  1:276392/50‑1 (MQ=255)
                                                ccACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAAt                                                      <  1:489815/40‑1 (MQ=255)
                                                 cACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGatat                                           >  1:318398/1‑50 (MQ=255)
                                                 cACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCAc                      <  1:3427797/71‑1 (MQ=255)
                                                 cACTTCACGTTGCAGTTTCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCa                       <  1:2254107/70‑1 (MQ=255)
                                                  aCTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCGATAAAGTTGATATGGGTATGCGTCATGCAGTCAc                      <  1:418498/70‑1 (MQ=255)
                                                   cTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGata                                            >  1:3684652/1‑47 (MQ=255)
                                                   cTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGata                                            >  1:697270/1‑47 (MQ=255)
                                                   cTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGata                                            >  1:1009015/1‑47 (MQ=255)
                                                   cTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATg                             <  1:1832776/62‑1 (MQ=255)
                                                      cACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGtt                  <  1:393856/70‑1 (MQ=255)
                                                      cACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGt                   <  1:937934/69‑1 (MQ=255)
                                                       aCGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTg                 >  1:36015/1‑70 (MQ=255)
                                                           tGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGt                                 >  1:2524893/1‑50 (MQ=255)
                                                               gTTTCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACa             <  1:3517316/66‑1 (MQ=255)
                                                                tttCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCg                                  <  1:479349/44‑1 (MQ=255)
                                                                tttCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGt       >  1:48569/1‑71 (MQ=255)
                                                                tttCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGt       >  1:3579864/1‑71 (MQ=255)
                                                                  tCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTa                                      >  1:3489822/1‑38 (MQ=255)
                                                                  tCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTa                                      >  1:2470612/1‑38 (MQ=255)
                                                                  tCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCa                           >  1:1598487/1‑49 (MQ=255)
                                                                  tCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATAc          >  1:3020371/1‑66 (MQ=255)
                                                                  tCAATAGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATAc          >  1:1831863/1‑66 (MQ=255)
                                                                   cAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCa                               <  1:119304/44‑1 (MQ=255)
                                                                   cAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGTTa     <  1:1938689/70‑1 (MQ=255)
                                                                    aaTGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGTTAAc   <  1:51123/71‑1 (MQ=255)
                                                                     aTGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGTTAACa  <  1:905302/71‑1 (MQ=255)
                                                                       |                                                                    
GTTTTTGGTTTCATCGTAAGCATCTTCAGAACCAAAAAGTTCTTTCGCCACTTCACGTTGCAGTTTCAATGGCAAGCCGGAACCAATAAAGTTGATATGGGTATGCGTCATGCAGTCACGGTTGAACATACAGTTAACA  >  minE/1079932‑1080070

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: