Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,002,764 C→T 27.0% L105L (CTG→CTA hybF ← protein involved with the maturation of hydrogenases 1 and 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,002,7640CT27.0% 119.1 / 40.2 103L105L (CTG→CTAhybFprotein involved with the maturation of hydrogenases 1 and 2
Reads supporting (aligned to +/- strand):  ref base C (47/28);  new base T (15/13);  total (62/41)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

ACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTGCGCA  >  minE/2002711‑2002821
                                                     |                                                         
aCGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTa                                            >  1:913565/1‑69 (MQ=255)
aCGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTa                                            >  1:631180/1‑69 (MQ=255)
 cGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGt                                             <  1:2823506/67‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:152550/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:925139/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:155327/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:159514/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:79707/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:1722441/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:1928253/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       <  1:2234414/49‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:2526491/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:2738869/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:3001179/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:303044/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:3181246/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       <  1:3235483/49‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       <  1:558287/49‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       <  1:554380/49‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:1332596/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:1384605/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:1154913/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCg                                                       >  1:848901/1‑49 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACg                                    >  1:3369920/1‑68 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:3253638/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:3287564/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2535056/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2569595/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2502626/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2582977/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:868584/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:1102398/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2707871/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:1217114/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:3348788/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:1310029/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:2773340/69‑1 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:1897620/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:66705/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:1114050/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:2927644/1‑69 (MQ=255)
         cacaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   >  1:1564286/1‑69 (MQ=255)
          acaTTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:1719602/68‑1 (MQ=255)
            attattCACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:2892/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCTGGTATc                                          >  1:3075269/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:2921727/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:713346/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:22028/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:1446533/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:1373094/1‑59 (MQ=255)
            attattAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATc                                          >  1:3430585/1‑59 (MQ=255)
               attaACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGc                             >  1:1091945/1‑69 (MQ=255)
               attaACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGc                             >  1:833382/1‑69 (MQ=255)
                  aaCTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGa                                                      >  1:2307184/1‑41 (MQ=255)
                  aaCTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGa                                                      >  1:100101/1‑41 (MQ=255)
                  aaCTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGa                                                      >  1:2776605/1‑41 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:2059084/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:2886337/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:112002/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:1607653/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:470714/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:750745/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:2608113/69‑1 (MQ=255)
                      ccGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGt                      <  1:2572441/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:366711/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:1084478/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:3287607/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:592671/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2676509/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2563282/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2539790/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2530891/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:1248302/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:796098/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2173091/69‑1 (MQ=255)
                           tATTCAACTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2135837/69‑1 (MQ=255)
                                aaCTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGc                                   <  1:2773342/46‑1 (MQ=255)
                                  cTTCAATACTTTTGACGATTAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCATGACag                 <  1:2869642/62‑1 (MQ=255)
                                  cTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGac          <  1:2882332/69‑1 (MQ=255)
                                  cTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGac          <  1:1679598/69‑1 (MQ=255)
                                  cTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGac          <  1:1687886/69‑1 (MQ=255)
                                   ttCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGaca         <  1:1350485/69‑1 (MQ=255)
                                   ttCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGaca         <  1:1150909/69‑1 (MQ=255)
                                   ttCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGaca         <  1:164303/69‑1 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:2583795/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:3407278/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:3403083/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:3126471/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:2840839/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:2791629/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:2185277/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:1682600/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:1612051/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:1482262/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:12744/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGacac                                     >  1:1061326/1‑39 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCCGGTATCGaca                                      >  1:3411605/1‑38 (MQ=255)
                                     cAATACTTTTGACGATTAGCGAATCGCAGGTATCGacac                                     >  1:2891127/1‑39 (MQ=38)
                                          cTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTGCGCa  >  1:2724612/1‑69 (MQ=255)
                                          cTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTGCGCa  >  1:1132576/1‑69 (MQ=255)
                                                     |                                                         
ACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTGCGCA  >  minE/2002711‑2002821

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: