Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,802,519 T→A 56.7% D196V (GAT→GTT)  qor ← quinone oxidoreductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,802,5190TA56.7% 8.1 / 38.0 30D196V (GAT→GTT) qorquinone oxidoreductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base A (12/5);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

CACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTT  >  minE/2802455‑2802577
                                                                |                                                          
cACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCctc                                                        >  1:1298263/1‑69 (MQ=255)
cACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCctc                                                        >  1:1146251/1‑69 (MQ=255)
      aTCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGa                                                 >  1:1042626/1‑70 (MQ=255)
          tACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCAc                                                   <  1:796923/64‑1 (MQ=255)
           tcaccacGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGATAGtt                                            <  1:1852394/69‑1 (MQ=255)
              ccacGCTCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATc                                                           >  1:284699/1‑52 (MQ=255)
              ccacGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATc                                                           >  1:1373204/1‑52 (MQ=255)
                 cgcgCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAAc                                      >  1:1286584/1‑70 (MQ=255)
                        tttCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCAcg                              >  1:1443872/1‑71 (MQ=255)
                          tCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGATAGTtaat                                         >  1:1278118/1‑58 (MQ=255)
                                cgcccGTGATCTCTTTTAACCGCTCGACCAGAACCTCTTCACGATAGtt                                            <  1:1848358/49‑1 (MQ=38)
                                    ggTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAAcc                                     <  1:11163/52‑1 (MQ=255)
                                     gTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCAc                                                   <  1:1299137/37‑1 (MQ=255)
                                     gTGATCTCTTTTAACCGCTCGACCAGAACCCCTTCAc                                                   <  1:1827696/37‑1 (MQ=37)
                                      tGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTg                                   >  1:677449/1‑52 (MQ=255)
                                       gATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgcg               >  1:1854718/1‑71 (MQ=255)
                                          ctctTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaataa                                       >  1:2114780/1‑44 (MQ=255)
                                           tctTTTAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGct           >  1:1106453/1‑71 (MQ=255)
                                            ctTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctc          >  1:1483354/1‑71 (MQ=255)
                                            ctTTTAACCGCTCGACCAGAACCTCTTCACGATAGTtaat                                         >  1:1045529/1‑40 (MQ=255)
                                             ttttAACCGCTCGACCAGAACCTCTTCACGATAGTtaata                                        <  1:1254863/40‑1 (MQ=255)
                                             ttttAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctct         >  1:1274986/1‑71 (MQ=255)
                                             ttttAACCCCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctct         >  1:2113440/1‑71 (MQ=255)
                                               ttAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAAc                                      >  1:384969/1‑40 (MQ=255)
                                               ttAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTgc       >  1:2162805/1‑71 (MQ=255)
                                               ttAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTgc       >  1:55439/1‑71 (MQ=255)
                                               ttAACCGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTgc       >  1:828125/1‑71 (MQ=255)
                                                   ccGCTCGACCAGAACCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgcg               <  1:1258873/59‑1 (MQ=255)
                                                    cGCTCGACCAGCTCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCttt  <  1:2206300/71‑1 (MQ=255)
                                                               aTCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTa                  <  1:1425983/44‑1 (MQ=255)
                                                                |                                                          
CACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTT  >  minE/2802455‑2802577

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: