Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,664 A→G 5.8% N396N (AAT→AAC)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,6640AG 5.8% 245.0 / 3.6 86N396N (AAT→AAC) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base G (2/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCT  >  minE/2093594‑2093731
                                                                      |                                                                   
aaaCACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCa                                                                     >  1:5084048/1‑71 (MQ=255)
 aaCACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCAt                                                                    >  1:4501693/1‑71 (MQ=255)
  acacGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCAt                                                                    >  1:2520186/1‑70 (MQ=255)
    acGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTAc                                                                 >  1:1049574/1‑71 (MQ=255)
    acGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTAc                                                                 >  1:4586301/1‑71 (MQ=255)
        ggTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:4699382/71‑1 (MQ=255)
        ggTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:1307244/71‑1 (MQ=255)
         gTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             >  1:1051328/1‑70 (MQ=255)
         gTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:109478/70‑1 (MQ=255)
          tGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCaa                                                           <  1:4429569/71‑1 (MQ=255)
            ccGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATg                                                         >  1:2188227/1‑71 (MQ=255)
               aGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCAc                                                      >  1:3772774/1‑71 (MQ=255)
               aGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCAc                                                      >  1:3369357/1‑71 (MQ=255)
                   gcggcgCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGa                                                  >  1:3882536/1‑71 (MQ=255)
                   gcggcgCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGa                                                  >  1:2599924/1‑71 (MQ=255)
                         cATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGt                                             <  1:1185664/70‑1 (MQ=255)
                          aTCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:1292180/53‑1 (MQ=255)
                          aTCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAc                                           >  1:5358735/1‑71 (MQ=255)
                           tCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCaa                                                           >  1:5364519/1‑54 (MQ=255)
                           tCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCaa                                                           >  1:759590/1‑54 (MQ=255)
                             gCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAAt                                                          >  1:3372029/1‑53 (MQ=255)
                             gCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTg                                        <  1:1356241/71‑1 (MQ=255)
                             gCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTg                                        >  1:374560/1‑71 (MQ=255)
                              ccATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAc                                           >  1:1770903/1‑67 (MQ=255)
                                 tCCTCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGAtt                                                <  1:5050246/59‑1 (MQ=255)
                                 tCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCAtt                                                                   <  1:1651732/40‑1 (MQ=255)
                                 tCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGAt                                                 <  1:5245108/58‑1 (MQ=255)
                                 tCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTATCTGTCCg                                    <  1:4469926/71‑1 (MQ=255)
                                 tCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGt                                       <  1:4249346/68‑1 (MQ=255)
                                  ccGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGAttt                                               >  1:3383283/1‑59 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGtt                                                              <  1:836306/43‑1 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGtt                                                              <  1:1549572/43‑1 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGtt                                                              <  1:2241053/43‑1 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                   >  1:1152682/1‑70 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                   >  1:1237253/1‑70 (MQ=255)
                                   cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                   >  1:602812/1‑70 (MQ=255)
                                    gtgcggcgTTCACCGCTTTTACTGCCGCAATGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGcg                                 <  1:140230/69‑2 (MQ=255)
                                    gtgcggcgTTCACCGCTTTTACTGCCGCAATGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGcg                                 <  1:2659374/69‑2 (MQ=255)
                                    gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:5290970/43‑1 (MQ=255)
                                    gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:2725174/43‑1 (MQ=255)
                                    gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:3806820/43‑1 (MQ=255)
                                    gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGc                                  >  1:367462/1‑70 (MQ=255)
                                    gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGc                                  >  1:4318269/1‑70 (MQ=255)
                                     cgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCa                                                            >  1:5461145/1‑43 (MQ=255)
                                        ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:37973/39‑1 (MQ=255)
                                        ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:4747398/39‑1 (MQ=255)
                                        ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCg                             <  1:2229390/71‑1 (MQ=255)
                                        ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCg                             <  1:3983882/71‑1 (MQ=255)
                                         gcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCgg                            <  1:47870/71‑1 (MQ=255)
                                          cgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCaa                                                           >  1:5091175/1‑39 (MQ=255)
                                          cgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCaa                                                           >  1:5543898/1‑39 (MQ=255)
                                          cgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGa                           >  1:4124898/1‑71 (MQ=255)
                                           gTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGAt                          >  1:2398074/1‑71 (MQ=255)
                                           gTTTACTGCTTACACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGAt                          >  1:4688571/1‑71 (MQ=255)
                                            tttACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:718600/35‑1 (MQ=255)
                                            tttACTGCTTTCACGGCATTAATGGCATTACGTTc                                                             <  1:4501103/35‑1 (MQ=255)
                                            tttACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCa                                 <  1:1033904/63‑1 (MQ=255)
                                             ttACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAAc                               <  1:2483797/64‑1 (MQ=255)
                                              tACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                        >  1:3354777/1‑70 (MQ=255)
                                              tACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                        >  1:1807330/1‑70 (MQ=255)
                                              tACTGCTTTCACGGCATTAATGGAATTACGTTCAATGCACCGGATTTGTACCTg                                        >  1:360010/1‑54 (MQ=38)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACggg                                                   >  1:2728630/1‑42 (MQ=255)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAcc                                          >  1:384749/1‑51 (MQ=255)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAcc                                          >  1:4997200/1‑51 (MQ=255)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                        >  1:428430/1‑69 (MQ=255)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                        >  1:3965734/1‑69 (MQ=255)
                                               aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGc                       <  1:1246911/70‑1 (MQ=255)
                                               aCTGCTTTCACGGCATTAAAGGCATTACGTTCAATGCACGGGATTTGTAcc                                          >  1:224919/1‑51 (MQ=255)
                                               aCTGCTTTCACGGAATTAATGGCATTACGTTCAATGCACggg                                                   >  1:1392820/1‑42 (MQ=255)
                                                 tGCTTTCACGGCATTAATGGCATTACGTTCAATGCACgg                                                    <  1:327343/39‑1 (MQ=255)
                                                   cTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGt                   <  1:870489/70‑1 (MQ=255)
                                                     ttCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAg                >  1:3275880/1‑71 (MQ=255)
                                                      tCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGc               >  1:5304007/1‑71 (MQ=255)
                                                         cGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCa                                 <  1:4757448/50‑1 (MQ=255)
                                                          ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa           <  1:4417386/71‑1 (MQ=255)
                                                          ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa           <  1:4238672/71‑1 (MQ=255)
                                                          ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa           <  1:1316666/71‑1 (MQ=255)
                                                            cATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg         >  1:1290078/1‑71 (MQ=255)
                                                             aTTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt        <  1:3518293/71‑1 (MQ=255)
                                                                aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt     >  1:62536/1‑71 (MQ=255)
                                                                 aTGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGaccg            <  1:4700496/63‑2 (MQ=255)
                                                                  tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa            <  1:4082327/62‑1 (MQ=255)
                                                                  tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa            <  1:2850388/62‑1 (MQ=255)
                                                                   ggCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt  >  1:5369122/1‑71 (MQ=255)
                                                                   ggCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt  >  1:72115/1‑71 (MQ=255)
                                                                   ggCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCt  >  1:2510150/1‑71 (MQ=255)
                                                                    gCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAg                     >  1:1661883/1‑51 (MQ=255)
                                                                      |                                                                   
AAACACGCGGTGCCGAGGTGCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCT  >  minE/2093594‑2093731

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: ‑