Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,682 G→C 12.8% I390M (ATC→ATG) ‡ tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,6820GC12.8% 239.9 / 21.5 93I390M (ATC→ATG) ‡tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base C (7/5);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCG  >  minE/2093613‑2093749
                                                                     |                                                                   
gcggcgCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGa                                                                    >  1:3882536/1‑71 (MQ=255)
gcggcgCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGa                                                                    >  1:2599924/1‑71 (MQ=255)
      cATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGt                                                               <  1:1185664/70‑1 (MQ=255)
       aTCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAc                                                             >  1:5358735/1‑71 (MQ=255)
          gCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTg                                                          >  1:374560/1‑71 (MQ=255)
          gCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTg                                                          <  1:1356241/71‑1 (MQ=255)
           ccATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTAc                                                             >  1:1770903/1‑67 (MQ=255)
              tCCTCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGAtt                                                                  <  1:5050246/59‑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)
                cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                                     >  1:1237253/1‑70 (MQ=255)
                cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                                     >  1:1152682/1‑70 (MQ=255)
                cgcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCgg                                                     >  1:602812/1‑70 (MQ=255)
                 gtgcggcgTTCACCGCTTTTACTGCCGCAATGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGcg                                                   <  1:2659374/69‑2 (MQ=255)
                 gtgcggcgTTCACCGCTTTTACTGCCGCAATGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGcg                                                   <  1:140230/69‑2 (MQ=255)
                 gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGc                                                    >  1:4318269/1‑70 (MQ=255)
                 gcgcgGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGc                                                    >  1:367462/1‑70 (MQ=255)
                     ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCg                                               <  1:3983882/71‑1 (MQ=255)
                     ggcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCg                                               <  1:2229390/71‑1 (MQ=255)
                      gcgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCgg                                              <  1:47870/71‑1 (MQ=255)
                       cgTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGa                                             >  1:4124898/1‑71 (MQ=255)
                        gTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGAt                                            >  1:2398074/1‑71 (MQ=255)
                        gTTTACTGCTTACACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGAt                                            >  1:4688571/1‑71 (MQ=255)
                         tttACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCa                                                   <  1:1033904/63‑1 (MQ=255)
                          ttACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAAc                                                 <  1:2483797/64‑1 (MQ=255)
                           tACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                                          >  1:1807330/1‑70 (MQ=255)
                           tACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                                          >  1:3354777/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:3965734/1‑69 (MQ=255)
                            aCTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCg                                          >  1:428430/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)
                                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:1316666/71‑1 (MQ=255)
                                       ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa                             <  1:4417386/71‑1 (MQ=255)
                                       ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa                             <  1:4238672/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:2850388/62‑1 (MQ=255)
                                               tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa                              <  1:4082327/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)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCa                                                   <  1:4158560/36‑1 (MQ=255)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                           <  1:1899907/60‑1 (MQ=255)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                           <  1:3464394/60‑1 (MQ=255)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                           <  1:5595479/60‑1 (MQ=255)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                           <  1:1350301/60‑1 (MQ=255)
                                                    ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                       <  1:4430465/64‑1 (MQ=255)
                                                    ttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAgg                           <  1:5253238/60‑1 (MQ=255)
                                                    ttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAt                <  1:2779624/71‑1 (MQ=255)
                                                     tACGTTCAATGCACGGGATTTGTACCTGTCCGGCa                                                   <  1:3204794/35‑1 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGt                                     >  1:2178079/1‑48 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTc                                    >  1:511074/1‑49 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTc                                    >  1:121262/1‑49 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGccc                               >  1:4994782/1‑54 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg              >  1:114394/1‑71 (MQ=255)
                                                      aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg              >  1:5653179/1‑71 (MQ=255)
                                                      aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAgac                                >  1:5432081/1‑51 (MQ=255)
                                                      aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGtt                         >  1:605005/1‑60 (MQ=255)
                                                      aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt                    >  1:346/1‑65 (MQ=255)
                                                      aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAtgg              >  1:1677202/1‑69 (MQ=255)
                                                      aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAtgg              >  1:944373/1‑69 (MQ=255)
                                                       cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc             <  1:1106021/71‑1 (MQ=255)
                                                       cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc             <  1:3410395/71‑1 (MQ=255)
                                                        gTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGt                                     <  1:2960225/46‑1 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                        >  1:4516193/1‑56 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                        >  1:4494426/1‑56 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                        >  1:3693864/1‑56 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt         >  1:3584056/1‑71 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt         >  1:1451747/1‑71 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt         >  1:3590660/1‑71 (MQ=255)
                                                           cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAt          <  1:4537438/70‑1 (MQ=255)
                                                            aaTGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt        <  1:1529390/71‑1 (MQ=255)
                                                              tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                          <  1:4015058/51‑1 (MQ=255)
                                                              tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                          <  1:226556/51‑1 (MQ=255)
                                                              tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                          <  1:4681619/51‑1 (MQ=255)
                                                               gCACGGGATTTGTACCTGTCCGGCAACCGGATCGc                                         <  1:1708365/35‑1 (MQ=255)
                                                               gCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                       <  1:236864/53‑1 (MQ=255)
                                                                cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc    >  1:1744041/1‑71 (MQ=255)
                                                                 aCGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGc                     >  1:216264/1‑53 (MQ=255)
                                                                  cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg  <  1:381/71‑1 (MQ=255)
                                                                  cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg  <  1:5343954/71‑1 (MQ=255)
                                                                   gggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc      >  1:1219149/1‑66 (MQ=255)
                                                                     gATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                       <  1:4799612/47‑1 (MQ=255)
                                                                     |                                                                   
GCGGCGCATCGCCATCCGCGCGGCGTTTACTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCG  >  minE/2093613‑2093749

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