Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,356,630 C→T 22.3% T6M (ACG→ATG)  mobA → molybdopterin‑guanine dinucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,356,6300CT22.3% 123.0 / 25.4 89T6M (ACG→ATG) mobAmolybdopterin‑guanine dinucleotide synthase
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base T (13/7);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGGGCGG  >  minE/2356564‑2356678
                                                                  |                                                
gagaCAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATgacgac                                                  >  1:1489653/1‑67 (MQ=255)
gagaCAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGacgacg                                                 >  1:1060862/1‑68 (MQ=255)
gagaCAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGacgacg                                                 >  1:2995597/1‑68 (MQ=255)
    cagacaCGTTAGCAGGGTCAATCCCACAATATAAGAGGCGATATCGGTGAATCTGATGACGACGATaaa                                            <  1:1486393/69‑2 (MQ=255)
    cagacaCGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATaaa                                            <  1:2221325/69‑2 (MQ=255)
    cagacaCGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATaaa                                            <  1:380789/69‑2 (MQ=255)
    cagacaCGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATaaa                                            <  1:3375455/69‑2 (MQ=255)
             tAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTtgt                                   <  1:2063517/69‑1 (MQ=255)
             tAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTtgt                                   <  1:3214397/69‑1 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATa                                              >  1:1224033/1‑53 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACATGCGTtgtg                                  >  1:1695304/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:1494543/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:2250958/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:2712657/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:151077/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:3280216/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:3425059/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:328378/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:1329738/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:1101013/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:621797/1‑65 (MQ=255)
                  gggTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTtgtg                                  >  1:53312/1‑65 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:948062/53‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:1137360/53‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:1828131/53‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:345811/53‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGt                                      <  1:309932/59‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGt                                      <  1:1225802/59‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGt                                      <  1:2987595/59‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGt                                      <  1:123824/59‑1 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCa                            >  1:249561/1‑69 (MQ=255)
                    gTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCa                            >  1:848357/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:1889655/52‑1 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAAc                                            <  1:2882150/52‑1 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGt                                      <  1:2412435/58‑1 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:3021685/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:3242880/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:991854/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:2280581/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:2276634/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:2105295/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:1513072/1‑69 (MQ=255)
                     tCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAg                           >  1:1485855/1‑69 (MQ=255)
                       aaTCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATaaa                                            <  1:3129400/50‑2 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:1905353/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:2814775/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:2974061/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:578118/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:1144121/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               >  1:2390313/1‑38 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGAt                                               <  1:621842/38‑1 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:1447772/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:1569996/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:3091247/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:1771690/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:931386/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:1051250/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:102370/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:2881339/1‑66 (MQ=255)
                                aTAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTaaa                   >  1:2279656/1‑66 (MQ=255)
                                     agagGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:2171052/69‑1 (MQ=255)
                                     agagGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:563708/69‑1 (MQ=255)
                                     agagGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:1800481/69‑1 (MQ=255)
                                     agagGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:2167132/69‑1 (MQ=255)
                                     agagGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:2551020/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:2541498/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:819440/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:844523/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:2592005/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:980045/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:1838447/69‑1 (MQ=255)
                                      gagGCGATATCGGTGAATCTGATGACGATGATAACAGCCGTTGTGCTGGCGGGCGGTAAAGCCAGACGa          <  1:370230/69‑1 (MQ=255)
                                       agGCGATATCGGTGATTCTGATGACGACGATAACAGGGGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:3223502/69‑1 (MQ=255)
                                       agGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:190977/69‑1 (MQ=255)
                                       agGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:1208326/69‑1 (MQ=255)
                                       agGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:1392450/69‑1 (MQ=255)
                                       agGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:3422925/69‑1 (MQ=255)
                                       agGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGaa         <  1:2695929/69‑1 (MQ=255)
                                           gATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:140503/63‑1 (MQ=255)
                                           gATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:2538487/63‑1 (MQ=255)
                                             tatCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:2431060/61‑1 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:2976374/1‑69 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:1643310/1‑69 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:2342489/1‑69 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:1204960/1‑69 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:1435461/1‑69 (MQ=255)
                                              atCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:3384491/1‑69 (MQ=255)
                                              atCGCTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:2088541/1‑69 (MQ=255)
                                               tCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:241732/1‑68 (MQ=255)
                                               tCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:2049384/1‑68 (MQ=255)
                                               tCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGggcgg  >  1:3345156/1‑68 (MQ=255)
                                                    gAAACTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:325913/54‑1 (MQ=255)
                                                              acgacgATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACg           <  1:1053985/44‑1 (MQ=255)
                                                                  |                                                
GAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGCCAGACGAATGGGCGG  >  minE/2356564‑2356678

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: