Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,240 T→C 29.1% E28E (GAA→GAG hydN ← formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,826,2400TC29.1% 171.2 / 77.3 110E28E (GAA→GAGhydNformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (45/33);  new base C (17/15);  total (62/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAA  >  minE/1826170‑1826307
                                                                      |                                                                   
tGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGAttt                                                                     >  1:2348571/1‑71 (MQ=255)
     acacCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCt                                                                 >  1:3577443/1‑70 (MQ=255)
     acacCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCt                                                                 >  1:4775503/1‑70 (MQ=255)
      cacCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCtga                                                               <  1:5383607/71‑1 (MQ=255)
          tttAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCtgat                                                              >  1:5689004/1‑68 (MQ=255)
                   aTGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCt                                                                 >  1:5190936/1‑56 (MQ=255)
                      gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcacgcacg                                               >  1:3602715/1‑71 (MQ=255)
                      gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACcacgcacg                                               >  1:1488008/1‑71 (MQ=255)
                      gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACg                                               >  1:3988531/1‑71 (MQ=255)
                          cgcgGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGaaa                                                        >  1:4173205/1‑58 (MQ=255)
                              gTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATgaga                                                          <  1:5510429/52‑1 (MQ=255)
                               tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTc                                                                   <  1:4209967/42‑1 (MQ=255)
                               tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                      <  1:612773/71‑1 (MQ=255)
                               tAAAAAAGTTTCCGGGGGCAGCGATGCACAGTCCTGATTTTc                                                                   <  1:4191548/42‑1 (MQ=255)
                                aaaaaaGTTTCCTGGGTCAGCGATGCACAGTCCTGATTTTca                                                                  >  1:2592121/1‑41 (MQ=39)
                                 aaaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCtgat                                                              <  1:5298269/45‑1 (MQ=255)
                                  aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:4945897/71‑1 (MQ=255)
                                   aaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGt                                  <  1:3883524/71‑1 (MQ=255)
                                    aaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTa                                 >  1:2947631/1‑71 (MQ=255)
                                       tttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCt                                                                 <  1:5361557/36‑1 (MQ=255)
                                       tttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCt                                                                 <  1:4469653/36‑1 (MQ=255)
                                        ttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACgg                              >  1:737801/1‑70 (MQ=255)
                                         tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAAcca                                                     >  1:96347/1‑45 (MQ=255)
                                         tCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac                                      >  1:1112419/1‑61 (MQ=255)
                                            ggggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAAc                          >  1:5412158/1‑70 (MQ=255)
                                              ggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCaa                       <  1:2852499/71‑1 (MQ=255)
                                                 cAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCa                                                 >  1:2704404/1‑42 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAAcca                                                     >  1:4041080/1‑36 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACTCACGCTACTACACAGGTACGGCAACCAATACa                   >  1:508002/1‑71 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                 >  1:3008819/1‑57 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                       >  1:4262611/1‑67 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCaa                       >  1:4462017/1‑67 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATAc                    <  1:3326430/70‑1 (MQ=255)
                                                  aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa                   >  1:1007127/1‑71 (MQ=255)
                                                   gCGATGCACAGTCCTGATTTTACTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa                   <  1:765359/70‑1 (MQ=255)
                                                    cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                 >  1:3295559/1‑71 (MQ=255)
                                                    cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                 <  1:4151053/71‑1 (MQ=255)
                                                    cGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt                 <  1:1007139/71‑1 (MQ=255)
                                                      aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTa               >  1:1894282/1‑71 (MQ=255)
                                                      aTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:4132488/51‑1 (MQ=255)
                                                      aTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:966424/51‑1 (MQ=255)
                                                      aTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTa               <  1:1973093/71‑1 (MQ=255)
                                                       tGCACAGTCCTGATTTTCCTGATGAGAAACCacgcacgc                                              <  1:1415566/39‑1 (MQ=255)
                                                        gCACAGTCCTGATTTTGCTGATGAGAAACCACGCACGCTACTTCACAGGTa                                 >  1:3382647/1‑51 (MQ=38)
                                                        gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                    >  1:1436344/1‑48 (MQ=255)
                                                        gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                    >  1:4692973/1‑48 (MQ=255)
                                                        gCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAg                                    >  1:4513815/1‑48 (MQ=255)
                                                         cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTAc                                >  1:2238228/1‑51 (MQ=255)
                                                         cacaGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTAc                                >  1:2244745/1‑51 (MQ=255)
                                                           caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc          <  1:3192130/68‑1 (MQ=255)
                                                           caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc          <  1:1102946/68‑1 (MQ=255)
                                                           caatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgc          <  1:4541114/68‑1 (MQ=255)
                                                            aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAgg                                   <  1:1078367/43‑1 (MQ=12)
                                                            aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAgg                                   <  1:2021130/43‑1 (MQ=12)
                                                            aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgca         >  1:649410/3‑70 (MQ=255)
                                                            aGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:181518/45‑1 (MQ=255)
                                                            aGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:1912060/45‑1 (MQ=255)
                                                             gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg                                   <  1:1449875/44‑1 (MQ=255)
                                                             gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt             >  1:4059290/1‑66 (MQ=255)
                                                             gTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt             >  1:133643/1‑66 (MQ=255)
                                                             atCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat        >  1:4393609/2‑69 (MQ=255)
                                                             atCCTGAATCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgcat        <  1:1165970/70‑3 (MQ=255)
                                                              tCCTGATTTTCCTTATGAGAAACCACGCACGCTACTTCACAGGTACGGc                             <  1:3634267/49‑1 (MQ=255)
                                                              tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                             <  1:5149652/49‑1 (MQ=255)
                                                              tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                             <  1:1982038/49‑1 (MQ=255)
                                                              tCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACg                               >  1:2740211/1‑47 (MQ=16)
                                                                cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg     >  1:880577/1‑71 (MQ=255)
                                                                cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg     >  1:3502090/1‑71 (MQ=255)
                                                                cTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAg     >  1:4335482/1‑71 (MQ=255)
                                                                cTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAt                  >  1:1394768/1‑58 (MQ=255)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:4574773/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:4594554/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:5273999/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:5233976/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:3370833/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:386318/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTa                                 >  1:168476/1‑43 (MQ=12)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAAccga                       <  1:1103046/53‑3 (MQ=255)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                  >  1:4947987/1‑58 (MQ=255)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                  >  1:1124460/1‑58 (MQ=255)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACAt                  >  1:2844205/1‑58 (MQ=255)
                                                                cTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg     >  1:1523518/1‑69 (MQ=255)
                                                                  gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt             >  1:3681554/1‑61 (MQ=255)
                                                                  gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACt             >  1:3799546/1‑61 (MQ=255)
                                                                  gATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc    >  1:5298515/1‑70 (MQ=255)
                                                                  gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                            <  1:4620319/46‑1 (MQ=16)
                                                                  gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                            <  1:2131427/46‑1 (MQ=16)
                                                                  gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCa                            <  1:4354397/46‑1 (MQ=16)
                                                                  gATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc    >  1:4850430/1‑70 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCa                            >  1:1634268/1‑45 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCaa  >  1:5036326/1‑71 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa   >  1:3737640/1‑70 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCa   <  1:2285167/70‑1 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACAt                  >  1:2524375/1‑55 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg           >  1:515577/1‑62 (MQ=255)
                                                                   aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc    <  1:1264433/69‑1 (MQ=255)
                                                                   aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc    <  1:3267245/69‑1 (MQ=255)
                                                                   aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc    <  1:4796910/69‑1 (MQ=255)
                                                                   aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCa   <  1:3064599/70‑1 (MQ=255)
                                                                      ttCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGc                             <  1:2303942/41‑1 (MQ=255)
                                                                      |                                                                   
TGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAA  >  minE/1826170‑1826307

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: