Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,036,505 G→A A632A (GCG→GCA rsxC → fused predicted 4Fe‑4S ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,036,5050GA98.5% 189.9 / NA 65A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/32);  minor base T (0/1);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036437‑1036572
                                                                    |                                                                   
gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:381999/1‑71 (MQ=21)
gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:1149904/1‑71 (MQ=21)
gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:1297638/1‑71 (MQ=21)
    aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                               >  1:2036628/1‑71 (MQ=21)
    aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                               <  1:1480103/71‑1 (MQ=21)
    aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                               >  1:1440995/1‑71 (MQ=21)
     aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                                                               >  1:899521/1‑70 (MQ=25)
     aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGc                                                              <  1:163705/71‑1 (MQ=21)
      ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                                             >  1:929001/1‑71 (MQ=21)
         gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGa                                                          <  1:2603528/71‑1 (MQ=21)
           agaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:494092/1‑60 (MQ=16)
                            gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       >  1:1040414/1‑71 (MQ=14)
                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       >  1:21547/1‑69 (MQ=14)
                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       >  1:2086049/1‑69 (MQ=14)
                              aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       >  1:25349/1‑69 (MQ=14)
                               aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:526550/1‑71 (MQ=25)
                               aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:2964950/1‑71 (MQ=25)
                                 gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                     >  1:2769306/1‑68 (MQ=25)
                                 gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc                                  <  1:5462/71‑1 (MQ=255)
                                    gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               <  1:2097257/71‑1 (MQ=255)
                                     ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                     <  1:3011547/64‑1 (MQ=25)
                                        tCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:2562654/69‑1 (MQ=255)
                                           aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                        <  1:813691/71‑1 (MQ=255)
                                           aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                        <  1:2028167/71‑1 (MQ=255)
                                            aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                       <  1:2271504/71‑1 (MQ=255)
                                            aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                       >  1:2542300/1‑71 (MQ=255)
                                            aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                       >  1:1939619/1‑71 (MQ=255)
                                            aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                       >  1:143746/1‑71 (MQ=255)
                                              ccGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                     >  1:652466/1‑71 (MQ=255)
                                              ccGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTc                     >  1:2123425/1‑71 (MQ=255)
                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:230157/62‑1 (MQ=255)
                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:2054068/62‑1 (MQ=255)
                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:2919543/62‑1 (MQ=255)
                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                            <  1:1631831/62‑1 (MQ=255)
                                                gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   >  1:3114743/1‑71 (MQ=255)
                                                 cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                  <  1:1603692/71‑1 (MQ=255)
                                                 cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                  <  1:1378975/71‑1 (MQ=255)
                                                 cTATTGCCCGTGACAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                  <  1:239609/71‑1 (MQ=255)
                                                  tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc                 >  1:1788272/1‑71 (MQ=255)
                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:2159730/70‑1 (MQ=255)
                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:2414377/70‑1 (MQ=255)
                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:1007476/70‑1 (MQ=255)
                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:182489/70‑1 (MQ=255)
                                                    ttGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcg               <  1:1005453/71‑1 (MQ=255)
                                                        ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCgcgc            >  1:1321535/1‑70 (MQ=255)
                                                        ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCgcgc            >  1:1882593/1‑70 (MQ=255)
                                                         cGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcg               <  1:1348854/66‑1 (MQ=255)
                                                         cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    >  1:2234690/1‑61 (MQ=255)
                                                         cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:680107/71‑1 (MQ=255)
                                                         cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:3026142/71‑1 (MQ=255)
                                                         cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa           >  1:1505786/1‑70 (MQ=255)
                                                         cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGACg                    >  1:1384817/1‑61 (MQ=255)
                                                           tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:475788/69‑1 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:2008137/58‑1 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:3011038/58‑1 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       >  1:3053220/1‑71 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       >  1:1910659/1‑71 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       >  1:2668298/1‑71 (MQ=255)
                                                            gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgc       >  1:2021390/1‑71 (MQ=255)
                                                             ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:2897750/57‑1 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:472945/60‑1 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:1838161/60‑1 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:1430169/60‑1 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATcc                <  1:1307023/60‑1 (MQ=255)
                                                              cAAAGCACGCAAGCTGGAACAGCAAAAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                  >  1:2205714/1‑58 (MQ=14)
                                                                aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgcc   >  1:2307078/1‑71 (MQ=255)
                                                                aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg        >  1:239391/1‑66 (MQ=255)
                                                                aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg        >  1:1737909/1‑66 (MQ=255)
                                                                 aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  >  1:878660/1‑71 (MQ=255)
                                                                    |                                                                   
GCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036437‑1036572

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: